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Section 1 


INTRODUCTION 


The National Aeronautics and Space Administration has undertaken 
a significant program of active control technology research and development. 

The feasibility of the use of this technology has been demonstrated in a num- 
ber of military research programs and by NASA (References 1, 2). These and 
other programs of research and study (References 3, 4) indicate that the prop- 
er application of Active Control Technology (ACT) can reduce the structural 
weight and drag and therefore bringing about substantial weight savings . In 
addition, active control technology can be expected to have direct benefits in 
terms of fatigue life of the vehicle, thereby making the technology addition- 
ally cost effective. The potential benefits to transport aircraft are obvi- 
ous in terms of fuel saving and extended vehicle service life. 

Under the sponsorship of NASA Langley,' the Calspan Corporation has 
undertaken a program of vehicle definition or mathematical modeling of the 
Total In-Flight Simulator (TIPS) aircraft in preparation for control system de- 
velopment and design of an ACT system for the TIPS. In addition, a' major part 
of the effort also involves design changes to the control system that would 
allow the ailerons to be operated collectively as well as differentially to en- 
hance the ability of the vehicle to perform the dual function of maneuver load 
control and gust alleviation. 


Ref. 1. NASA Symposium on Advanced' Control Technology and Its Potential for 
Future Transport Aircraft. NASA TM X-3409, August 1976. 

2. Flood, R. V.: "A Summary of the Application of Active Controls Techno- 

logy on the ATT System Studies," NASA Symposium on Advanced Control 
Technology and Its Potential for Future Transport Aircraft. 

NASA TM X-3409, August 1976. 

3. Rynaski, E. G. and Weingarten, N. C. : "Flight Control Principles for 

Control Configured Vehicles." Air Force Flight Dynamics Laboratory 
Report AFFDL-TR- 71-154, January 1972. 

4. Stockdale, C. K. and Poyncer, R. D. : "Control Configured Vehicle Ride 

Control System (CCV RCS)." Air Force Flight Dynamics Laboratory Report 

' AFFDL-TR-73-83, July 1973. 


The major emphasis of this study, however, was to mathematically predict the 
rigid body and flexible equations of longitudinal motion using the FLEXSTAB 
program . 


This report is organized in three main sections. Section 2 de- 
scribes the elastic and aerodynamic analysis of the TIFS airplane using the 
FLEXSTAB program. Included in this analysis is a definition of the vehicle 
geometry., the mass and stiffness distribution, the calculated mode frequencies 
and mode shapes and the resulting aerodynamic equations of motion of the flex- 
ible vehicle. Comparisons with data taken during the ground vibration and 
flight tests of the vehicle are made. Section 3 describes the control and 
instrumentation system that presently exists on TIPS., including performance 
and flight test data on the control and instrumentation system. Section 3 also 
contains an analysis of the existing elevator servo system and proposed mod- 
ifications to the system to improve its performance. 

Conclusions and recommendations are included at the ends of Sections 
2 and 3. Large parts of these sections are devoted to an evaluation of the 
accuracy and usefulness of existing data and recommendations for the accumula- 
tion of additional flight and ground data. 
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Section 2 


ELASTIC AND AERODYNAMIC ANALYSIS 


2.0 THE TIFS AIRPLANE 

The USAF/Calspan Total In-Flight Simulator (TIFS) airplane, as 
shown in Figure 1 is probably the most complete and versatile control con- 
figured and variable stability aircraft ever constructed. TIFS has been de- 
signed to have independent control of all six degrees of rotational and trans- 
lations of freedom of rigid body motion. The aircraft is equipped with large 
servo-driven direct lift flaps for independent control of lift forces. Side 
forces are independently controlled through large vertical surfaces built into 
the wings. The propellers are used for longitudinal force control and the con- 
ventional ailerons, elevator and rudder are used for roll, pitch and yaw con- 
trol. All of these control surfaces are driven by full authority electrohy- 
draulic servos. Aside from the side force control surfaces, TIFS' most dis- 
tinctive feature is the addition of a complete and removable second cockpit on 
the nose of the airplane. The capability of either masking the presently in- 
stalled simulation cockpit, or easily removing it and adding a completely new 
one is unique. The instrument displays are readily changed. The airplane has 
been flown with a wheel controller, a side stick controller and with a center 
stick. The wheel controller, center stick and rudder pedals are equipped with 
variable feel systems including nonlinear functions. A more detailed descrip- 
tion of this aircraft as well as its capability is given in Reference 5. 

A number of research projects and simulations of specific airplanes 
have been performed using the Air Force TIFS airplane. These include simula- 
tions of the space shuttle, the B-l bomber and the Concorde supersonic trans- 
port; a flight research program to explore use of side force control in both 
automatic and manual modes, to counter crosswinds during crosswind landings; 
and flight demonstrations of multiplex techniques for flight control. 

.Ref. 5. Reynolds, P. A., Wasserman, R., Fabian, G. J. and Motyka, P. R. : 
"Capability of the Total In-Flight Simulator (TIFS)." Air Force 
Flight Dynamics Laboratory, AFFDC TR-72-39, July 1972. 
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craft into linear equations of motion and to investigate control system- 
airframe interactions and b) the computer program was immediately available 
to CaTspan through the Air Force Flight Dynamics Laboratory. 

FLEXSTAB provided the following combined computational capability 
that made its application suitable for this program. 

1) Capability to accurately analyze a very large number of struc- 
tural modes. 

2) Capability to retain the static deflection of all modes deleted 
in the analysis (Residual Elastic Formulation of the Equations 
of Motion) . 

3) A first order unsteady aerodynamic effect. 

4) Capability to improve the completeness of the geometric model 
and therefore improve the aerodynamic model (e.g., engine na- 
celles and side force surfaces) . 

5) Ability to very accurately compute symmetric, antisymmetric 
and combined dynamic properties of an airframe. 

In the discussions to follow concerning the finite element struc- 
tural and aerodynamic models of the TIFS prepared for input to Flexstab, the 
detailed description of analytical procedures employed in FLEXSTAB will be 
omitted. Detailed documentation of the program is provided by Reference 7. 

A simplified block diagram showing the sequence of computations 
in the FLEXSTAB system, as used for this ACT TIFS study, is shorn in Figure 2. 


Ref. 7. *'A Method for Predicting the Stability Characteristics of Control 

Configured Vehicles". Air Force Flight Dynamics Laboratory Report 
AFFDL-TR- 74-91, November 1974. 
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t Figure 2 FLEXSTAB SYSTEM COMPUTATIONAL SEQUENCE 








Each block is a separate program or module and is run separately. The start- 
ing point is the Geometry Definition (GD) program which defines the geometry 
of the vehicle. Structural properties of the vehicle are input in terms of 
mass and inertia distributions into the Internal Structural Influence Coeffi- 
cient (ISIC) program. This program uses thin beam .theory to compute structur- 
al influence coefficients for the vehicle. Computation of the normal vibra- 
tion modes is performed in the Normal Modes (NM) program using the ISIC pro- 
gram outputs. The steady and unsteady aerodynamic properties of the vehicle 
are computed in terms of aerodynamic influence coefficients in the Aerodynamic 
Influence Coefficient (AIC) program. The aerodynamic theory employed is a lin- 
ear first order approximation to unsteady, dnviscid, unseparated subsonic or 
supersonic flow equations. Computation of stability and control derivatives, 
analysis of vehicle static and dynamic stability and formation of the vehicle 
equations of motion are performed in the Stability Derivative and Static Sta- 
bility (SD+SS) program. The final program employed in the TIFS analysis was the 
Linear Systems Analysis (LSA) program. The classical techniques of open or 
closed loop analysis are performed in this program, i.e. Bode plots, Nyquist 
plots , and root locus . 

2.2 TIFS GEOMETRY DEFINITION 


The baseline TIFS configuration consisted of the General Purpose 
canopy with side-force surfaces installed and gear up, i.e. the standard TIFS 
configuration. The resulting geometric representation used in FLEXSTAB 
is shown in Figure 3. 

The following geometric approximations were found to be either 
necessary or desirable in formulating the aircraft model: 

1) The point at which the wing dihedral changes from 4. 83° (inner 
wing) to 6.5° (outer wing) was moved from Convair Station 8 
(a distance of 4.2909 meters or 168.933 inches perpendicular 
to the aircraft plane of symmetry) to the engine nacelle cen- 
terline (3.81 meters or 150 inches from the plane of symmetry) . 

8 

} 



Y-fiXIS METE 



Figure 3 TIFS FLEXSTAB GEOMETRIC REPRESENTATION 




This was desirable to simplify the FLEXSTAB model but had the 
effect of raising the wing panels outboard of Convair Station 
8 by a vertical distance o_f ,.0142. .meters, or .558 inches-. 


2} In FLEXSTAB the geometry of the inner wing panel is described 
with respect to a local axis system (X N , Y N , Z^) whose X^ axis 
is parallel to the X axis shown on Figure 3 and whose X^Y N plane 
coincides with a mean surface for this thin body. Plane X^Yj^ 
is oriented with .respect to the Reference Axis System by a 
rotation about the X^ axis of magnitude equal to the dihedral 
angle. The manufacturer’s inner panel chord plane of the TIFS 
airplane does ( not lie in the X^Y^ plane because of a 4° angle 
of incidence. Consequently, the geometry for the FLEXSTAB 
model of TIFS was developed with the incidence removed by 
assuming the actual wing rotated by 4° about a line perpendicu- 
lar to the plane of symmetry and located at the intersection of 
the elastic axis of the wing with the plane of symmetry. Be- 
cause the wing has dihedral, this rotation also had the effect 
of moving the wing tip .115 meters or 4.528 inches forward. 

3) The side-force surface fairing was ignored. 

4) The engine nacelle was taken to be a closed body, i.e. no flow 
through the nacelle was assumed. The TIFS engine is a turbo- 
prop engine which derives only approximately 10% of its total 
thrust from flow through the engine. 

i 

5) The fuselage and engine nacelle were represented with circu- 
lar cross sections with. the correct cross sectional area and 
correct area centroid location in the Z direction. 

6) The only planar lifting surface with a thickness distribution 
was the wing. 
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7) The interference body cross section which is used in FLEXSTAB 
to determine the interference effects between the horizontal 
stabilizer, vertical tail and fuselage is required to be iden- 
tical to one used for the wing/ fuselage intersection. The 
interference body used in the TIFS model had a larger radius 
than the actual fuselage in the region of the horizontal tail 
and thus enclosed a portion of the horizontal tail. Consequent- 
ly, the geometric positions of the horizontal and vertical tails 
were changed in accordance with the recommendations in the FLEX- 
STAB User’s Manual. The thin bodies used to represent the - hor- 
izontal and vertical tails were moved outboard and upward re- 
spectively enough to expose the same surface areas as obtained 
with the actual fuselage. Correct sweep angles of leading and 
trailing edges and correct tip chords were maintained in this 
process . 

8) Exposed areas used in FLEXSTAB are listed in the table below: 

2 2 

Area m ft 


Horizontal Tail per side (with 
elevator) 

9.46 

101.815 

Vertical Tail (with rudder, 
without dorsal fin) 

13.0 

139.84 

Elevator per side (with tabs) 

3.34 

35.87 

Direct Lift Flap per side (aft 
of hinge line) 

4.19 

45.0 

Aileron per side (aft of hinge 
line) 

2.09 

22.6 


9) The TIFS was assumed to have a vertical plane of symmetry. 

Actual TIFS modifications of the forward fuselage are not sym- 
metric due to the crew access tunnel to the evaluation cock- 
pit. 
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It is recognized that the horizontal tail aspect ratio is effec- 
tively increased when the thin bodies representing the left hand and right hand 
sides of the tail are moved outboard as discussed above in Item 7. This geo- 
metric change would have some effect on rolling moments developed by the hor- 
izontal tail but is believed to be of secondary importance in considering long- 
itudinal stability and control. .Of the geometric approximations made, the shape 
of the fuselage cross section and the treatment of the engines and nacelles are 
likely to most seriously affect the accuracy of the analytical results. The for- 
ward fuselage of the TIFS is roughly the shape of a rounded triangle. As a re- 
sult, the lifting effectiveness and therefore the destabilizing effect on long- 
itudinal static stability may be greater than predicted by FLEXSTAB. To correct- 
ly model the engines both the propeller and direct thrust slipstreams and their 
corresponding effects on wing, fuselage and tail would have to be incorporated 
in some way into the FLEXSTAB model. A study of techniques to accomplish this 
task was beyond the scope of this contract and therefore not attempted. 

2.3 DESCRIPTION OF MASS AND STIFFNESS PROPERTIES OF USAF/TIFS 

In this section a description’ is given of the mass and stiffness 
properties of TIFS. The basic TIFS airplane is the Convair 580 which is the 
version of the Convair powered by two Allison S01-D13 turboprop engines. 

Mass and stiffness distributions for TIFS were derived by review and analysis 
of data for the reciprocating engine Convair 340 airplane, data describing the 
Allison prop-jet engine conversion, and data concerning structural modifica- 
tions made in the development of TIFS. 

Mass and stiffness parameters given in this section were the basis 
for inputs prepared for the Structural Influence Coefficient Program (ISIC) of 
FLEXSTAB which was used in computing structural vibration modes. A printout 
of the computer inputs to the ISIC program is reproduced in Appendix A. In 
general, the basic parameters which are presented are limited to those which 
could be utilized within the modeling constraints of the ISIC program. How-* 
ever, the fore and aft bending stiffness of the wing and some engine inertia 
and stiffness parameters which cannot be handled directly by ISIC are also 
included. 
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2.3.1 Wing Parameters 


The TIFS airplane was received by Calspan prior to modification 
as a reciprocating-engine C-131B aircraft. TIFS direct-lift wing flaps and 
side-force surface actuators were installed on this airplane before power- 
plant conversion. Inertia parameters for making a wing flutter analysis of 
this interim configuration were derived by correcting data for the Convair 
340 given in Reference 8. Inertia parameters for making a wing flutter 
analysis of this interim configuration were derived by correcting data for 
the Convair 340 given in Reference 8. Intertia data for 30 sections of the 
TIFS modification. Contributions due to the TIFS direct-lift flaps and side- 
force surface modifications were distributed to the same stations and com- 
bined with the distributions for the 340 wing. 

Subsequently, these data were updated for the completely modified 
TIFS wing including the engine change and were used in the wing flutter analy- 
sis for the Air Transport TIFS (AT/TIFS) . These distributions of wing mass, 
mass unbalance about the elastic axis, and moment of inertia about the elastic 
axis are also applicable for USAF/TIFS and are used in this report. Table 1 
lists the distributions for the zero fuel condition. The table also shows the 
weight distribution for a fuel load of 1451. 5 Kg/side (3200 Ibs/side) 
obtained by interpolation of data for other fuel conditions. The parameters 
given are in English units as obtained from the original data. They are 
converted to metric units in Table 2 which were used in the FLEXSTAB solution. 


Ref. 8. Johnson, N. S., and' Krueger, N. A.: "Wing Stress Analysis - Model 

340." Report ZS-340-006, Part I, Consolidated Vultee Aircraft Corp., 
31 October 1951, NASA CR-158891. 


13 



The wing elastic .axis is assumed to be at 36.5% of thejving chord 
which is the position used in the Convair 340 flutter analysis (Reference 9} . 
This elastic axis is midway between the front and rear spars. The c.g. of 
the fuel in the tanks between the spars is assumed to be on the wing elastic 
axis. The effective moment inertia of the fuel was assumed to be 2/5 of the 
values computed treating the fuel as a rigid body. Table 2 presents the total 
inertia distributions for the wing including both structural and fuel weights. 
An equivalent dumbbell mass representation is also given for the distributed 
weight at each listed station except Sta. 3.825 m (Sta. 150.53 in.). The 
equivalent properties are converted to metric units which were used in the 
FLEXSTAB representation of TIFS. 

The treatment of the engine sprung weight and nacelle weight which 
are assumed to act at Sta. 3.825 m (150.53 in.) are taken up in a later sec- 
tion. The properties of the side-force surface vanes are not included in 
Table 1 and are also discussed separately in a later section. However, the 
side-force surface actuators, fairing, etc. which are integral parts of the 
wing are included with the wing weight. 


Ref. 9 Pancu, C. D., and Hiroshige, K. : "Model 340 Main Fixed and Movable 

Control Surface Flutter Analysis." Report ZU-340-003, Consolidated 
Vultee Aircraft Corp., 28 August 1951. NASA CR- 158893. 
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TABLE 1 


DISTRIBUTIONS OF WING WEIGHT, STATIC MOMENT ABOUT ELASTIC 
AXIS AND MOMENTS OF INERTIA ABOUT ELASTIC AXIS 

NOTES: 1) Elastic axis at 36.5% chord. . _ . . . 

2) Values of parameters are per side and are in English units as obtained from original 
data. These parameters are converted to metric units in Table 2. 


Sta. -in. 

Weight (lbs) 1 

Static Moment 

(in. lbs) 

Moment of Ineri 

tia (lbs. in^) 

(Chord 

Zero Fuel* 


Zero Fuel* 


Zero Fuel* 


plane) 

Cond. 

Fuel 

Cond. 

Fuel 

Cond. 

Fuel 

28.38 

502 


5,036 


480,243 


48.38 

566 


5,941 


534,359 


150.53 


3,319 

270 

[ -146,977 

0 

J 10,600,042 

23,316 





\ 



(Engine sprung. 

150.53 


3,122 


[-400,552 


[ 58 3 663 , 798 

weight) 

205.75 

501 

1,358 

7,831 

0 

422,119 

110,382 

259.53 

337 

1,009 

4,431 

0 

271,366 

68,028 

308.53 

372 

454 

6,596 

0 

249,886 

24,950 



516 

109 

- 591 

0 

332,523 

5,165 

363.53 










26 


- 734 


20 , 935 


401.25 


129 


1,470 


76,908 


436.00 

99 


490 


35,353 


478.07 

131 


874 


42,796 


531.16 

94 


776 


35,074 


594.37 

91 


754 


25,311 


Total 

9,805 

3,200 






* Values as determined for AT/TIFS less effect of wing tip booms not on USAF/TIFS. 




TABLE 2 


REPRESENTATION _OF WING INERTIA PROPERTIES WITH DUMBBELL MASSES 
[F uel Mass V 1451 . .5 Kg/Si de_ (Fuel Weight = 3200 lbs /Side)] 


E.A. (36_. 5-fi_ chord) representation of inertia 



Parameters from Table 1 in English Units | 

Equivalent 

Properties- 

Metric Units 

Sta. 

(in) 

Sta. 

Cm) 

Weight 
(lbs) . 

Unbalance 
about EA 

(.in. lbs) . 

Moment of 
Inertia 
about EA 
fibs in 2 ) 

* Lm Z 

(m) 

X, 

L Mi 

On) 

I 

0.0 

0.0 








28.38 

0.721 

502 

5,036 

480,243 

-.488 

.998 

113.8 

88.38 

2,245 

566 

5,941 

534,359 

-.467 

1.000 

128.4 

150.53 

3.823 

6,711 

-547,529 

69,287,156 




205.75 

5.226 

1,859 

7,831 

532,501 

-.309 

.523 

421.6 

259.53 

6.592 

1,346 

4,431 

339,394 

-.316 

.478 

305.3 

308.53 

7.837 

826 

6,596 

274,836 

-.214 

-.619 

187.3 

363.53 

9.234 

651* 

-1,325 

358,623 

-.646 

.542 

197.6 

401.25 

10.192 

129 

1,470 

76,908 

-.259 

.838 

29.3 

436.00 

11.074 

99 

490 

35,353 

-.338 

.589 

22.5 

478.07 

12.143 

131 

874 

42,796 

-.257 

.596 

29.7 

531.16 

13,491 

94 

776 

35,074 

-.234 

.653 

21.3 

594.37 

15.097 

91 

754 

25,311 

-.157 

.578 

20.6 

Total 


13,005 







*Does not include SFS vanes. 
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Figure 4 AREA MOMENT OF INERTIA I vv VS. SPANWISE POSITION ON WING 

A A 


Comparatively few modifications were made to the basic wing 

structure in the development of TIFS and the wing stiffness of TIPS was 

assumed to be the same as. for the reciprocating-engine version - of the Convair. 

Figure 4 compares structural moments of inertia used in the Convair 340 

flutter analysis (Reference 9) to the average section properties computed 

with compression on the upper or lower surfaces as given in Reference 10 . 

The out of plane bending stiffness used in the TIFS analysis is based on the 

I curve from Reference 9. Figure 5 provides the I„- moment of inertia 
xx ia 

from Reference 10 which is the basis for the assumed fore and aft bending 
stiffness used for ACT TIFS. Figure 6 presents a plot of the torsional 
stiffness factor J from Reference 9 which was used to obtain the torsional 
stiffness of TIFS. 

In analyzing TIFS, the stiffness curves or Figures 4, 5 and 6 were 
approximated by step curves which were constant between concentrated mass 
stations. Table 3 lists the constant values used between stations in this 
approximation. 

The wing vertical shears on the TIFS airplane are primarily car- 
ried into the fuselage by wing fittings located approximately at body line 
(1.397 m) . Figure 7 taken from the Convair 440 maintenance manual presents 
sketches of the wing to fuselage attachment fittings. In developing a model 
for the wing to fuselage attachment, the fuselage vertical deflection was as- 
sumed equal to the wing vertical deflection at the attachment fittings rather 
than at the centerline of the airplane. However, the wing cannot be assumed 
to be cantilevered from the fuselage at the attachment fittings because there 
is a change of wing slope at the fittings due to bending deflections in the 
wing carry-through structure. 

Figure 8 indicates how the effect of wing bending inboard of the 
fuselage attachment fittings was simulated in the structural model. As shown 

Ref. 10. Krueger, J. A., and Johnson, N. S. : "Wing Stress Analysis — Model 

340." Report ZS-340-006, Part II, Consolidated Vultee Aircraft 
Corp., 18 February 1952. NASA CR-158892. 


18 



SCiXMD 'Hboa 3Q 

m. mvz BKNioreo 


to 



) I 1 I ! I 1 i I ! I 1 I I 1 1 

o 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 

Wing Station, meters 


Figure 5 STRUCTURAL MOMENT OF INERTIA 1^ jVS. SPANWISE POSITION 
(UPPER SURFACE COMPRESSION) (Reference ,10) 
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Wing Station, meters 


Figure 6 TORSIONAL STIFFNESS FACTOR J VS. SPANWISE POSITION ON WING (Reference 9 ) 


TABLE 3 


WING STIFFNESS 


Station j 

El , Beam Bending 

EL.,, Fore $ Aft Bend. 


JG, Torsional 

Stiffness 

in. 

m 

1010 lbs in? 

10 8 N m 2 

10?-0 lbs in? 

10 8 N m 2 

10^0 lbs in 1 . 2 

107 N m 2 

0.0 

0.0 

6.747* 

1.936* 





28.38 

0.721 

6.747* 

1.936* 





55.00 

1.397 

5.614 ■ 

1.611 

10.043 

2.882 

2.420 

6.945 

88.38 

2.245 

4.429 

1.271 

7.952 

2.282 

2.002 

5.745 

150.53 

3.823 

3.142 

.902 

6.386 

1.833 

1.348 

3.869 

205.75 

5.226 

2.189 

.628 

4.491 

1.289 

.9548 

2.740 

259.53 

6.592 

1.494 

.429 

3.203 

.919 

.6853 

1.967 

308.53 

7.837 

.9476 

.272 

2.112 

. 606 

.4810 

1.380 

363.53 

9.234 

.6386 

.1833 

1.494 

.429 

.3385 

.9714 

401.25 

10.192 

.4841 

.1389 

1.082 

.311 

.2503 

.7183 

436.00 

11.074 

.3502 

.1005 

.876 

.251 

.1887 

.5415 

478.07 

12.143 

.2575 

.0739 

.464 

.133 

.1309 

.3757 

531.16 

13.491 

.1802 

.0517 

.278 

.0798 

.0770 

.2210 

594.37 

15.097 

.144 ' 

.041 

.206 

.0591 

.0470 

.0135 

635.64 

16.145 



• 





*Represented 
by an equiv- 
alent sec- 
tion in com- 
puter model. 


,, I B.L. 55.0 



Figure 7 ’WING TO FUSELAGE ATTACHMENT 


(a) 


SYMMETRIC BENDING 
OF TIFS WING 


£ B.L. (55 in. = 1.347 m) 



\ 


^ • LM SLOPE DUE 

^TIFS (EI) tif$ TO BENDING 

(El } TT rc BENDING STIFFNESS OF 

CARRY-THROUGH STRUCTURE 


(b) WING MODEL IN 
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on diagram (a) , symmetric bending of TiF^-wgs analyzed assuming that all wing 
shear is taken into the fuselage at the attachment (BL 55 in = 1.397 m) . 
However., it is assumed that no wing trending moments are taken out at these 
fittings and that the bending moments for the left and right wings balance 
each other at the centerline. 

The beam model representing the fuselage is at the centerline in 
diagram (b) of Figure 8. The inbound portion of the wing is simulated by two 
beams designated A and B. Beam A which extends from the centerline almost to 
B.L. (1.397 m) is assumed very stiff so that its deflections are negligible. 

A low bending stiffness is selected for beam B in order that the change in 
slope of the wing model due to bending will be the same as that for the 
actual TIFS wing at the spanwise location of the attachment fittings. How- 
ever, beam B is so short that its deflection does not result in appreciable 
difference between the fuselage vertical deflection and the vertical deflec- 
tion of the wing at the attachment. 

2.3.2 Engine Sprung Weight, Engine Mounts Installation and Nacelle 

The FLEXSTAB structural model developed for TIFS included modeling 
of motion of the engine sprung mass due to vibration isolator deflections. 
Results of vibration mode computation given in this report are limited to the 
case where the nacelle >and engine isolators were 'assumed rigid. This was ac- 
complished by assuming beams used to represent vibration isolators to be very 
stiff. Consequently, some of the parameters given below were not required 
but have been included to indicate how engine motions could be simulated in 
future computations with FLEXSTAB. 

Figure 9 is a drawing of the engine mounts installation taken 
from the Allison/Convair maintenance manual. It shows the engine sprung 
weight system consisting of the propeller, transmission and engine which are 
supported by three Lord vibration isolators, two f orward and one aft. Perti- 
nent parameters are given below for the inertia properties, mount geometry, 
and isolator stiffnesses of this system. 
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1 . FORWARD ENGINE MOUNT SUPPORT-UPPER 

2. FORWARD ENGINE MOUNT SUPPORT-LOW ER 

3. FORWARD ENGINE MOUNT FRAME ASSEMBLY -STA. 24.08 

4. FRAME INSTALLATION ~$TA. 33.10 

5. FRAME ASSEMBLY *• STA. 94.4S8 

6. REAR ENGINE MOUNT ASSEMBLY 

7. FORWARD ENGINE MOUNT ASSEMBLY 


Figure 9 ENGINE MOUNTS INSTALLATION 
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Inertia Properties of Sprung Weight 


M 

"'pitch 

I 

yaw 

Z roll 

Ref. c.g.= 


1416.1 Kg (W = 3122 lbs) 
2128 Kg m 2 
2127 Kg ra 2 
295.9 Kg m 2 


(0.483 m) aft of reduction gear mount pad 
(0.178 m) above power section centerline 


i 


Mount Geometry 

F.S. 239.67 in. (6.088 m) F.S. 386.98 in. 

(9.829 m) 
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Vibration Mount Stiffnesses 


Right front mt. (15.375 in. = 0.391m right of thrust line): 

K vert (static) = 13,000 lbs/in. = 2.28 x 10 6 N/m 

K ^ (.dynamic) = 17,000 lbs/in. = 2.98 x 10 6 

vert 

Kfo re g a £t (static) = K^ atera j = 16,000 lbs/in. = 2.80 x 10° 

<W « aft Cd >' namic: ' - ^lateral = 21 ’ 000 lbs/in ' " 3 ' 68 * 1° 6 


Left front mt. (15.375 in. = 0.391m left of thrust line): 
Same as for right front mount. 


Aft mt . (On thrust line) : 

K ^ (static) 

■ vert 

K vert ( dynamic) 


lateral (static3 


= 10,000 lbs/in. = 1.75 x 10 6 N/m 
= 13,000 lbs/in. = 2.28 x 10 6 
= 12,000 lbs/in. = 2.10 x 10 6 
= 15,000 lbs/in. = 2.63 x 10 6 


Since only simple beam finite elements are used in the ISIC pro- 
gram of FLEXSTAB, it was not possible to model all six rigid body degrees of 
freedom of the sprung weight as restrained by the vibration isolators. How- 
ever, the vertical and side translation degrees of freedom could be repre- 
sented by the model indicated on the sketch below where the sprung weight 
system is treated as a concentrated mass. 


SPRUNG MASS 
(1416.1 kg) 


WING ELASTIC AXIS 




_E 


■A 


A 


(814 kg) 


29.82 in. 
(0.757 m) 


r 


B 


98.48 in. (2.502 m) 


-128.3 in. (3.259 m) 



•(814 kg) 


16.575 in. 
(0.421 m) 


A — Beam representing vertical and lateral 
stiffness of vibration isolators 

B — Rigid beam connecting dumbbell masses 
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The effect of the concentrated sprung mass of 1416.1 Kg (W = 3122 
lbs) was subtracted from the total inertia properties distributed to wing Sta. 
150.53 in. (3.823m)- as- follows and the remaining properties were represented 
by the dumbbell masses indicated in the sketch. 



Weight 

[Mass] 

Unbalance 
About EA 

Moment of 
Inertia 
About EA 

Total at Sta. 150.53 
(from Table 2) 

6711 (lbs) 
[3044 (Kg) ] 

-547529 (in. lbs) 
Q- 6308 (m Kg)J 

69287156 (lbs in 2 ) 
[20276 (Kg m 2 )] 

Effect of concentrat- 
ed engine mass 

3122 (lbs) 
[l416.1 (Kg)] 

-400553 (in. lbs) 
[- 4615 (m Kg)] 

51390898 (lbs in 2 ) 
[l5039 (Kg m 2 )] 

Remaining properties 
which were repre- 
sented by dumbbell 
masses 

3589 (lbs) 
[l627.9 (Kg)] 

-146976 (in. lbs) 
Q 1693 (m Kg)] 

17896257 (lbs in 2 ) 
[ 5237 (Kg m 2 )] 


The required vertical bending stiffness of beam "A" in the sketch 
to represent the dynamic properties of the vibration isolators is as follows: 


K vert ( Total ) = (2 x 17,000) + 13,000 = 47,000 lbs/in. 
(El ) = vertical bending stiffness of beam "A" 

= i 3 K/3 = (47,000) (29. 82) 3 /3 = 4.15 x 10 8 lbs. in 2 

= 1.194 x 10 6 N m 2 


The uncoupled vertical vibration frequency of the sprung weight on beam A is 
then. 


vert 


1 / (47000) (386) 

27 r H 3122 


12.13 Hz. 


Similar computations for the lateral stiffness give. 


K lat (Total) = (2 X 21000) + 15000 = 57000 lbs/in. 
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(EI^) = Required lateral bending stiffness of beam A 

= 5.03 x 10 8 lbs in 2 (1.448 x 10 6 N m 2 ) 

f = 13.36 Hz = uncoupled lateral vibration frequency of 

1 cl L 

sprung weight on vibration isolators 


2.3.3 Fuselage Parameters 

The development of TIFS involved extensive structural modifica- 
tions affecting both the stiffness and weight distributions of the fuselage. 
Installation of special electronic equipment made further changes in weight 
distribution. 

Table 4 presents the fuselage mass distribution used in the pre- 
sent report in preparing inputs for FLEXSTAB. It was obtained by redistri- 
buting the panel point masses to fuselage stations located at 1.5 m intervals 
along the fuselage axis. As noted, the distribution does not include the 
weights of the horizontal and vertical tails. The last column gives the 
distributed masses for one side of the fuselage as required for the ISIC 
inputs. 


The representation of the fuselage by concentrated masses along 
its axis is an approximation since there are mass moments of inertia for each 
fuselage section associated with each concentrated mass. However, these sec- 
tion moments of inertia are not given because they are neglected in the ISIC 
program. 


Figure 10 shows plots of the effective structure moments of inerti 
used in computing fuselage vertical and lateral bending stiffnesses while 
Figure 11 gives a plot of the factor J used in computing fuselage torsional 
stiffness. These curves were approximated by a series of steps in treating 
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TABLE 4 

FUSELAGE MASS DISTRIBUTION* 


Sta. (in) 

X* 

(in) 

X* 

(m) 

Weight (lbs) 
(Both Sides) 

M(kg) 
(One Side) 

-117.82 

19.69 

.5 

432.1 

98.0 

- 78.44 

59.06 

1.5 

604.1 

137.0 

- 19.38 

118.12 

3.0 

1708.9 

387.6 

+ 39.68 

177.18 

4.5 

1866.2 

423.2 

+ 98.74 

236.24 

6.0 

2922.1 

662.7 

+157.80 

295.30 

7.5 

1824.4 

413.8 

+216.86 

354.36 

9.0 

1324 

300.3 

275.92 

413.42 

10.5 

1116.9 

253.3 

334.98 

472.48 

12 

1422.5 

322.6 

394.04 

531.54 

13.5 

1278.5 

290.0 

453. 10 

590.60 

15 

2022.2 

458.6 

512.16 

649.66 

16.5 

1197.5 

271.6 

571.22 

708.72 

18 

1744.6 

395.7 

630.28 

767.78 

19.5 

1826.3 

414.2 

689.34 

826.84 

21 

1619.9 

367.4 

748.40 

885.96 

22.5 

1469.7 

333.3 

807.46 

944.96 

24 

1136.4 

257.7 

866.52 

1004.02 

25.5 

994.2 

225.5 

949.11 

1086.61 

27.6 

155.8 

35.3 


26666 . 0 

6047.8 ' 


*Does not include horizontal or vertical tail 
% - Pus. Sta. = 137.5 in. = 3,493 m 
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Notes: 

(1) Values are totals for both sides 
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-2 0 2 4 6 8 10 12 14 16 18 20 22 24 

Fuselage Station, meters 


Figure 10 TIFS EFFECTIVE AREA MOMENTS OF INERTIA USED IN COMPUTING FUSELAGE 
VERTICAL AND LATERAL BENDING STIFFNESSES 
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‘Figure 11 TIFS FUSELAGE TORSIONAL STIFFNESS FACTOR, J {(Reference 11) 




the fuselage as a series of finite element beams each one having constant stiff 
ness. Table 5 indicates the stiffnesses used in this approximation. 

It will be noted that the vertical, lateral, and torsional stiff- 
nesses given in the last three columns of Table 5 are for both sides of the 
fuselage. The ISIC program requires fuselage stiffnesses for both sides as 
inputs while the mass inputs are for only one side. 

2.3.4 Parameters for Side Force Surfaces (SFS) 

Table 6 presents SFS inertia distributions which were used for 
both the upper and lower SFS vanes. Both the unbalance and moments of inertia 
are given about the SFS pivot or hinge line which was taken to be the SFS elas- 
tic axis. A dumbbell representation of this data is also presented for use in 
FLEXSTAB. 


The structural model for the SFS's used finite elements which 
permitted torsion about their elastic axes and transverse bending. However, 
the following stiffness data are only given for future reference and were not 
used for the structural mode computations discussed in this report. The com- 
putations were carried out assuming the side force surfaces to be effectively 
rigid and substituting very high stiffnesses for the finite elements. 

Bending and torsional stiffnesses of the SFS's were approximated 
by a series of steps giving constant stiffness parameters between the stations 
to which the inertia parameters were distributed. These results are given in 
Table 7. 


The computed cantilever bending frequency for the SFS was 52.68 Hz 
while vibration test results (Reference 12) gave a SFS bending or rocking 


Ref. 12. Balcerack, J. C. and White, Jr., R. P.: "Ground Vibration Tests, 

AF/TIFS Convair 580 Airplane", RASA Report 71-06, May 1971. 

NASA CR- 15 8895. 
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TABLE 5 - FUSELAGE STIFFNESSES (TOTALS FOR BOTH SIDES) 


X* 

Vertical Bending 
El yy x 10" 8 

Lateral Bending 
El,, x 10 -8 

Torsion 
JG x 10“ 8 

M 

(Nm 2 ) 

(Nm 2 ) 

(Nm 2 ) 

0.0 

2.953 

2.923 

2.798 

1.5 

3.254 

3.978 

3.626 

3.0 

3.616 

4.882 

3.637 

4.5 

4.309 

5.273 

3.637 

6.0 

6.268 

5.665 

3.660 

7.5 

7.051 

6.027 

3.716 

9. 

5.966 

5.544 

4.052 

10.5 

6.509 

4.82i 

5.451 

12 

7.352 

4.761 

7.163 

13.5 

7.443 

5.484 

7.275 

15 

6.629 

5.454 

5.540 

16.5 

5.635 

5.743 

4.533 

18. 

5. 123 

5.002 

4.205 

19.5 

4.821 

4.761 

2.451 

21. 

4.550 

4.520 

1.242 

22.5 

4.369 

4.369 

.8394 

24 

4.249 

4.249 

.5260 

25.5 

4.188 

4.188 

.3917 

27.623 





*X - Fus. Sta. (m) = 3.493 m 
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TABLE 6 


SFS MASS DISTRIBUTION (Upper or Lower Vane) 


. SFS E.A. (Assumed at the centerline of the SFS pivot. The 
centerline of the pivot intersects the manufacturing chord 
plane at wing Sta. 363.534 in (9.234m), This intersection 
is 28.53 in (0.724m) forward of the wing elastic axis.) 



Chord Plane 


/ 


Wing E.A. (FLEXSTAB 
Model of TIFS) 


TIFS Wing E.A. 


«« 0.719 m = B**| 

■«* 0.725 m £» 


SFS 

STA. 

■ 

1 

Mom. About H.L. 
(Tail Heavy) 
(in. lbs.) 

Mom. of Inertia 
About H.L. 
(lbs. in ) 

FLEXSTAB 

PARAMETERS 

X L 

"a 

(m) 

L "l 

(m) 

gjgg 

1 1 

(in) 

(M) 

7.75 

0.197 

26.41 

63.9 

3157 

-0.209 

0.332 

5.989 

15.00 

0.381 

15.73 

138.7 

3883 

-0.106 

0.554 

3.568 

25.00 

0.635 

34.56 

291.1 

3490 

0.075 

0.353 

7.838 

42.50 

1.080 

0.77 

81.5 

1800 

-0.074 

0.458 

2.443 

55.00 

1.397 

11.79 

82.0 

1500 

-0.049 

0.402 

2.674 

70. 00 

1.778 

6.26 

39.. 5 

1020 

-0.122 

0.442 

1.420 

82.50 

2.096 

1.12 

6.3 

180 

-0.146 

0.431 

0.259 


24.186 
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TABLE 7 


SFS STIFFNESS PROPERTIES 
(UPPER OR LOWER VANES) 


SFS 

Sta 

El, 

Bending Stiffness 

GJ., Torsional Stiffness 

(in.) 

Cm) 

(lbs 

in. ) 

(N m 2 ) 

(lbs in 2 ) 

(N m 2 ) 

0 


.294 

.if 

x 10 

4 + 

8.43 x 10 

4.7 x 10 8 

1.35 x 10 6 

7.75 

0.197 

5.4 

X 

1 — 1 
o 

CO 

1.55 x 10 6 

4.05 x 10 8 

1.16 x 10 6 

15.00 

0.331 

4.7 

x 10 8 

1.35 x 10 6 

3.30 x 10 8 

9.47 x 10 5 

25.0 

0.635 

3.8 

x 10 8 

1.09 x 10 6 

2.30 x 10 8 

6.6 x 10 5 

42.5 

1.080 

2.0 

x 10 8 

5.7 x 10 5 

1.0 x 10 8 

2.9 x 10 5 

55.0 

1.397 

1-2 

x 10 8 

3.4 x 10 5 

.45 x 10 8 

1.29 x 10 5 

70.0 

1.778 

.05 

x 10 8 

' 1.4 x 10 4 

.08 x 10 8 

4 

2.3 x 10 

82.5 

2.096 

.04 

x 10 8 

1.1 x 10 4 

.07 x 10 8 

2.0 x 10 4 

85.0 

2.159 







^ Effective value accounting for flexibility of bearing housing attachment. 
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frequency of only 16.7 Hz. It is believed that this discrepancy was caused by 
flexibility in the attachment of the SFS bearing housing to the wing and was 
accounted for by reducing the effective stiffness of the inboard section of 
the SFS model. 

No correction was made to the computed SFS torsional stiffness be- 
cause the measured torsional frequency of 52.68 Hz is above the frequency range 
of interest in the current program. 

2.3.5 Horizontal and Vertical Tail Parameters 

Horizontal and vertical tail parameters are based on data given in 
Reference'll. Inertia parameters from this report were represented by dumb- 
bells for convenience in preparing ISIC inputs to the FLEXSTAB program. 

The elastic axis of the horizontal tail is modeled -by two segments 
as shown in the sketch on Table 8. A rigid segment perpendicular to the plane 
of symmetry extends out to the stabilizer-fuselage attachment where it joins 
a second segment at Fuselage Sta. 836.532 in. (21.248 m) and Body line 
38.938 in. (0.989 m) . The second segment is inclined aft at an angle of 7.07° 
and is at approximately 49% of the portion of the chord forward of the eleva- 
tor hinge line. 

Table 8 indicates the positions and properties of the dumbbells 
used to represent the inertia parameters of the horizontal tail. The most 
inboard dumbbell is located 3.5 in (.089 m) inboard of the connection between 
the elastic axis segments. The remaining dumbbells are all on the outer seg- 
ment. 


Figure 12 presents plots of the structural moments of inertia, I, 
and torsional stiffness factor, J, used in computing the out-of-plane bending 
stiffness (El) and torsional stiffness (GJ) of the horizontal tail. The 

Ref. 11. Rodriguez, M. B.: "Commercial Total In-Flight Simulator Empennage 

Flutter Analysis," General Dynamics Convair Division Report 
GDE-DDE-69-001 , 13 March 1969. NASA CR-15SS94. 


37 



TABLE 8 


HORIZONTAL TAIL INERTIA PROPERTIES* 



i 27.60 .701 -17.81 -.452 22.7 10.30 

11.22 0.285 38.13 .968 - 1.82 -.0462 40.1 18.18 

39.75 1..010 35.49 .901 - .61 -.0155 41.75 18.94 

68.22 1.733 30.28 .769 - 2.66 -.0676 32.25 14.62 

96.70 2.456 27.07 .688 - 1.47 -.0373 27.0 12.24 

125.19 3.180 22.54 .573 - 2.39 -.0607 19.7 8.93 

153.62 3.902 19.06 .484 - 3.54 -.0899 13.9 6.30 

182.14 4.626 15.64 .410 - 3.36 -.098 6.2 7.82 

192.87 4.899 17.60 .447 - 2.95 -.0749 2.16 ,98 

93.31 

* Based on data from Ref. 11..- 

^ Assumed to act as dumbbell parallel to and rigidly attached to fuselage 
3.5 in (.089m) inboard of indicated fuselage attachment. 

** In FLEXSTAB model of TIFS attachment point moved outboard to body line 
64.605 in (1.641m). 
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(Reference 11) 



Figure 12 ACT TIFS HORIZONTAL TAIL STIFFNESS 
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abcissa on this figure is the body line or distance outboard of the plane of 
symmetry as measured on the actual TIFS airplane. In the FLEXSTAB model of 
TIFS, the -root of the horizontal tail is .moved outboard- for aerodynamic reasons 
by extending the rigid segment of the elastic axis by 25.667 in. (0.652 m) . 

This extension acts as a portion of the fuselage and has no effect on the com- 
puted symmetric vibration modes. 

The elastic axis for the vertical tail is represented by two 
segments as shown in the sketch on Table 9. The inboard rigid segment runs 
up perpendicular to the fuselage axis to the vertical tail-fuselage attach- 
ment. It joins the second segment at Fuselage Sta. 827.4 in. (21.016 m) 
and waterline 94.788 in. (2.408 m) . The second segment which is inclined 
aft at 12.57° is at approximately 48% of the portion of the chord forward of 
the rudder hinge line. Table 9 gives the positions and properties of the 
dumbbells used to represent inertia parameters of the vertical tail. The 
structural moments of inertia, I, for out-of-plane bending and the torsion 
constant, J, for the vertical tail are plotted on -Figure 13 as a function 
of the waterline. The bending stiffness (El) and torsional stiffness (GJ) 
of the vertical tail were computed using these factors. 

Mass Summary 


A summary of the component masses for the TIFS airplane which 
were given previously is listed below: 


1/2 Wing (excluding engine 

and nacelle) = 2 x 1427.4 

One engine and nacelle (one side) 

1/2 Fuselage 

Upper SFS (one side) = 2 x 24.186 

Lower SFS (one side) = 2 x 24.186 

Horizontal Tail (one side) = 2 x 93.31 

Vertical Tail (one side) = 2 x 59.93 

Total (one side) 


= 2854.8 Kg 
= 3044.1 Kg 
= 6047.8 Kg 
= 48 . 4 Kg 
= 48.4 Kg 
= 186.6 Kg 
= 119.9 Kg 
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= 12,350.0 Kg 
=(27,227.1 lbs) 



TABLE 9 


VERTICAL TAIL INERTIA PROPERTIES* 


12.57 


Assumed Elastic Axis 


Fuselage Attachment 
F.S. = 827.4 in.= 21.016m 
W.L. = 94.788 in. = 2.408m 


Yl m z Y 


Dumbbell Masses 


0.973 m 


£ Fuselage (W.L. = 56.5 in. = 1.435m) 


y 

n 


yj. 



M Z 

M 1 

(in) 

(m) 

mssm 

Cm) 

(in) 

(m) 

msxsi 

9.630 

.245 

52.24 

1 

.327 

-30.30 

- .770 

30 

32.911 

.836 

56.56 

1 

.429 

-22.47 

-0.571 

24.4 

57.757 

1.467 

46.15 

1 

.172 

-19.86 

-0.428 

20.1 

83.582 

2.123 

38.31 

0 

.973 

-14.73 

-0.374 

16.35 

108.019 

2.744 

32.46 

0 

.824 

- 9.32 

-0.237 

13.75 

136.014 

3.455 

24.77 

0 

.629 

- 1.42 

-0.036 

15.95 

164.300 

4.173 

18.36 

0 

.466 

+ 1.40 

+0.025 

11.55 


13.61 

11.07 

9.12 

7.42 

6.24 

7.23 

5.24 


59.93 


*Based on data from Ref. 11 (Values are for one side) . 

In FL'EXSTAB model of TIFS attachment point moved up to water line 124.703 in. 
= 3.167m to put vertical tail above fuselage interference body. 
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2 3 4 5 6 7 


TIFS WATER LINE (m) 

Figure 13 ACT TIFS VERTICAL TAIL STIFFNESS 
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2.3.6 Results of Frequency and Mode Shape Computation Performed With the 
FLEXSTAB Program 

Table 10 summarizes natural frequencies for the first seven sym- 
metric vibration modes of the TIFS airplane which were computed with the FLEX- 
STAB program assuming high rigidity for the nacelles and SFS’s. These modes 
are believed sufficient to represent dynamic responses in the frequency range 
where control by collective ailerons, DLF's and/or elevators is feasible. The 
effect of deflections in higher frequency modes was treated quasi-statically 
by utilizing the Residual-Elastic option of FLEXSTAB. 

The computed modes are "free-free" modes for the entire airplane 
and involve motions of all structural components. Nevertheless, it is some- 
times helpful to characterize the motions in the various modes although admit- 
tedly a somewhat arbitrary procedure. This has been done on the table which 
indicates characteristic and auxiliary motions. 

Figures 13 to 19 present isometric sketches of the various modes 
which were drawn from modal deflections found in the FLEXSTAB computations. 
Some intuition concerning the relative effectiveness of the various control 
surfaces for controlling these seven modes can be obtained from an examination 
of" the drawings. In general, the vibration amplitudes of the TIFS airplane 
are comparatively small and can be treated adequately with linear analyses. 
However, the vibration amplitudes on the sketches have been greatly exaggerat- 
ed in order to indicate small motions which would otherwise not be visible. 
Motions of the SFS's and vertical tail have been omitted to simplify the 
drawings . 


The output of the Normal Modes program of FLEXSTAB gives the modal 
displacements and rotations at each aerocentroid used in the aerodynamic model. 
The motion at a few of these aerocentroids are indicated with arrows on the 
figures . 
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TABLE 10 


TIFS SYMMETRIC VIBRATION MODES COMPUTED WITH FLEXSTAB 



FREQUENCY (Hz) 



SYM 

RIGID 



MODE 

NACELLE 


SECONDARY MOTION 

NO. 

§ SFS'S 

CHARACTERISTIC MOTION 

1 

3.02 

1st Wing Bending 

Fus. § H. T. Bending 

2 

4.65 

1st Wing Torsion 

Wing, Fus., § H. T. Bending 

3 

6.09 

1st Fus. Bending in Phase 
With H. T. Bending 

Wing Bending 

4 

6.90 

H. T. Bending Out of Phase 
With 1st Fus. Bending 

Wing Bending 

5 

9.01 

2nd Wing Bending 

Fus. § H. T. Bending 

6 

15.10 

2nd Fus. § H. T. Bending 

Wing Bending and Torsion 

7 

17.25 

Outer Wing Panel Bending 
and Torsion 




















The lowest frequency mode (f^ = 3.02 Hz) which is shown on Figure 
14 is a typical 1st wing bending mode with the wing tips moving out-of-phase 
with the forward and -aft fuselage and horizontal tail motion. It would be 
expected that the collective ailerons and DLF's would be most effective in 
controlling this mode. Elevator control forces would also do appreciable work 
in a virtual deflection in this mode and might also be used for control. 

The second symmetric mode ( f ^ = 4.65 Hz) was characterized as 
being primarily first wing torsion because its frequency was close to the one 
computed for uncoupled torsion assuming no bending and the’ wing clamped at its 
root. However, it can be seen from. Figure is that considerable wing, fuselage 
and horizontal tail bending movements are involved in the coupled motion. 

Again it appears that collective aileron, DLF, and elevator inputs could all 
affect the modal response. 

Figure 16 shows the relative deflections for mode 3 (f^ = 6.09 Hz). 
Fuselage bending in this mode is such that the deflections of the nose and aft 
ends of the fuselage are in phase with wing tip deflections whereas they were 
out-of-phase for the first mode. Also the relative amplitude of the horizon- 
tal tail bending is very large and is in phase with deflections of the aft end 
of the fuselage. 

In contrast. Figure 17 shows the fourth mode (f^ = 6.90 Hz). The 
tip bending of the horizontal tail is out-of -phase with aft fuselage deflec- 
tion resulting in a nodal line running across the horizontal tail. However, 
this nodal line is quite' far inboard and it would be expected that the ele- 
vator would be effective in controlling this mode as well as mode 3 shown on 
Figure 16. 


Mode 5 (9.01 Hz) shown on Figure 18 is characterized as a typical 
wing second bending mode. Again appreciable fuselage and horizontal tail 
bending are present. It appears that the direct lift flaps would be the most 
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effective surfaces for controlling second wing bending. Nodal lines cross 
both the aileron and elevator which would tend to decrease their effectiveness 
for controlling this mode. 


The sixth mode (f = 15.10 Hz) which is shown on Figure 19 is 
characterized by 2nd fuselage bending and 2nd horizontal tail bending although 
considerable wing bending and torsion are also present. The seventh mode 
(f^ ~ 17.25 IIz) which is shown on Figure 20 is even more difficult to chai'ac- 
terize. However, it is noted that although significant bending and torsional 
deflections are found on the outer wing panel, deflections are small on other 
parts of the airplane. 

2.3.7 Comparison of TIFS Vibration Modes Computed Using FLEXSTAB with Previ- 
ous Results 

Table 11 presents a comparison of TIFS vibration frequencies com- 
puted by the FLEXSTAB program with previous results. The frequencies given 
for the preliminary ACT study were computed by coupling together modal results 
which had been determined for the wing, fuselage, and horizontal tail separate- 
ly. In contrast the FLEXSTAB frequencies were obtained from a solution for 
the eigenvalues of a structural model for the entire airplane which was based 
on somewhat different fundamental parameters than used in the earlier study. 

In both solutions, the nacelles and SFS's were assumed to have high rigidity. 

The listed ground vibration test modes correspond to the princi- 
pal peaks found in a symmetric wing survey up to 11.5 Hz carried out with 
vertical excitation applied by shakers located at the wing tips. The wing 
fore and aft bending and engine vibration modes were not excited in these 
tests or included in the FLEXSTAB results in Table 11. Vibration modes assoc- 
iated with these degrees of freedom are believed to have a comparatively small \ 
influence on the response of the TIFS airplane to atmospheric turbulence. 
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TABLE 11 


TIFS MEASURED AND COMPUTED STRUCTURAL FREQUENCIES 






FREQUENCY (Hz) 

, 

SYM 

MODE 

NO. 

MODE 

CHARACTERIZATION 

PRELIM. ACT 
STUDY 
ref. 6 

FLEXSTAB 

COMPUTATION 

GROUND 
VIBRATION 
TEST ref. 12 

FLIGHT TEST IN 
ATMOSPHERIC 
TURBULENCE* ** 

1 

1st Wing Bending 

2.80 

3.02 

3.3 

3.12 

2 

1st Wing Torsion 

4.58 

4.65 

4.9 

5.0 

3 

Fus. £j H. T. Bending 

5.56 

6.09 

6.6 

7.13 

4 

Fus. § H. T. Bending 

6.89 

6.90 

7.5 

■k* 

5 

2nd Wing Bending 

8.26 

9.01 

9.7 

9.99 

6 

2nd Fus. § H. T. 
Bending 

16.59 

15.10 

NA 

17.13 

7 

Wing Outer Panel 
Mode 


17.25 

NA 

NA 


* Peaks in power spectrum of normal acceleration at pilot station 
(insufficient data to characterize modes) 

** No peak observed 




The last column in the table presents results obtained from flight 
tests in atmospheric turbulence. The frequencies listed are the frequencies 
corresponding to the power spectrum peaks of normal acceleration at the pilot's 
station. Insufficient data were present to characterize the vibration modes 
associated with the peaks in the spectrum, but they have been associated some- 
what arbitrarily with the computed and ground vibration test modes. Further- 
more, it should be pointed out that the frequencies of the spectrum peaks can 
be influenced by zeroes in the gust input/normal acceleration transfer func- 
tion arising from unsteady aerodynamic effects at gust entry. Thus they would 
not be expected to agree exactly with frequencies obtained from the poles of 
the transfer function. 

Figures 21 to 25 present a comparison of vibration nodal lines 
determined from the mode computation by FLEXSTAB to those measured in the 
ground vibration tests (Reference 12 ) . Results for the first symmetric wing 
bending mode presented on Figure 21 show the measured and computed node lines 
to lie in the nacelle region and to be fairly close together. 

Computed and measured nodal lines for the first symmetric torsion 
mode which are presented on Figure 22 are of the same general form, but are 
widely separated in the region inboard of the nacelle. Nodal lines which are 
located off the wing surface are determined as if narrow rigid wing sections 
perpendicular to the wing elastic axis extended both forward and aft of the 
actual planform. 

The peculiar behavior of the nodal lines inboard of the nacelle 
can be explained by the spanwise variation of the pitching motion of the wing 
sections. This phenomenon will be discussed considering the situation depicted 
on Figure 5. Here there is a very small nose down pitch angle at the wing 
root, but torsional deformation of the wing structure causes nose-up sectional 
pitch angles a short distance outboard. 

The small nose-down pitch angle at the wing root in conjunction 
with an upward vertical deflection results in a section node which is too far 
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VERTICAL BENDING (FLEXSTAB) 
(PLOTTED IN HORIZONTAL PLANE) 


Figure 23 NODAL LINES FOR 1st SYMMETRIC FUSELAGE BENDING MODE 




Figure 24 NODAL LINES FOR 1st SYMMETRIC HORIZONTAL TAIL BENDING MODE 




AIRPLANE 


Figure 24 NODAL LINES FOR 2nd SYMMETRIC WING BENDING 


forward to plot on Figure 22. Moving slightly outboard, the section node goes 
forward toward infinity as the nose-down pitch angle approaches zero. Still 
further outboard where a nose-up pitch angle develops and where there is still 
an upward vertical deflection at the elastic axis, the section node is aft of 
the wing. Thus as the nose-up section pitch angle increases from zero with 
spanwise position, the nodal line comes in from infinity aft of the wing 
towards the trailing edge as indicated on Figure 22. It is evident that in 
the region where the section pitch angles are small that small differences 
between the computed and measured pitch angles results in large differences 
between the computed and measured nodal lines. 

Figure 23 shows nodal lines obtained for the mode characterized 
as first symmetric fuselage bending. The shape of the fuselage vertical bend- 
ing curve as computed by FLEXSTAB is plotted in the horizontal plane in order 
to indicate how the fuselage node was determined. The measured fuselage node 
was somewhat forward of the one obtained from the FLEXSTAB solution, but it 
should be noted that both theory and ground vibration test results indicated 
vertical motion of the aft fuselage and horizontal tail to be in phase. 

Theoretical and measured nodal lines on the wing are also compared 
on Figure 23 for the mode characterized as first symmetric fuselage bending. 
They are found to be in general agreement. The vibration test results as re- 
ported in Reference 12 did not define the position of the nodal line which is 
outside the planform of the outer wing panel, but it is believed to be approx- 
imately parallel to the one obtained from the FLEXSTAB results. 

Nodal lines are shown on Figure 24 for the mode characterized by 
first symmetric horizontal tail bending. Both FLEXSTAB and vibration test re- 
sults indicate nodal lines crossing the inboard portion of the horizontal tail. 
The line obtained in vibration tests is at a greater angle to chord line than 
the one predicted by the theoretical solution indicating greater horizontal 
tail torsion. This difference might be associated with the fact that the 
FLEXSTAB solution treated the stabilizer and elevator as an integral unit. 
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The FLEXSTAB solution for the first horizontal tail bending mode 
indicated a fuselage node approximately midway between the wing and horizon- 
tal tail. Presumably such a node was also present in the vibration tests, but 
its position was not recorded. 

Measured and computed nodal lines are shown on Figure 24 which .run 
diagonally across the outer wing panel and are in substantial agreement. How- 
ever, again the vibration test results do not define the nodal line position 
off the wing planform. 

Figure 24 indicates a nodal line observed at the wing root in the 
vibration test as reported in Reference 12. It is believed that the result is 
spurious, but insufficient original data were available to verify or disprove 
this conclusion. 

Results for the 2nd symmetric wing bending mode which are presented 
on Figure 25 .shows a good agreement between the computed inboard nodal line 
and the one found in the ground vibration tests. The measured and computed 
outboard nodal lines are found at approximately the same spanwise station. 
However, the nodal line obtained in the ground vibration tests crosses the wing 
at a greater angle relative to the centerline of the airplane indicating more 
torsion. Possibly this difference could be attributed to aileron or DLF de- 
flections which were not considered in the FLEXSTAB analysis, but might have 
occurred during the vibration tests Olthough not recorded in Reference 21) . 

2.3.8 Discussion of Results of Vibration Analysis 

In the preceding sections vibration mode computations have been 
discussed and compared with frequencies and mode shapes measured in ground 
vibration tests of the TIPS airplane. It is concluded that the analysis has 
provided a reasonable description of the principal structural vibration modes 
throughout the frequency range contemplated for active control of TIPS by 
collective aileron, direct lift flap and elevator inputs. All of the com- 
puted frequencies are within 8% of those measured in the ground vibration 
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tests . The analytical results have been obtained using the original estimates 
for the mass and stiffness distributions for the airplane and these parameters 
have not been adjusted or "tweaked" to improve agreement between test and 
■theory. 


It is recognized that several possible sources of error existed 
in obtaining these results. As discussed earlier, some known limitations were 
placed on the analysis by the use of the Internal Structural Influence Coeffi- 
cient Program (ISIC) of FLEXSTAB for the generation of coefficients used in 
the normal modes computations. ISIC places certain constraints on the repre- 
sentations of slender bodies (e.g. fuselages and nacelles) and on thin bodies 
(i.e. the wings, SFS's, tail surface). The effects of these limitations on 
the TIFS structural mode computations are summarized below. 

1. It was necessary to model fuselage inertia properties with 
concentrated masses spaced along the fuselage axis. Fuse- 
lage panel effects of the mass distribution away from the 
fuselage axis could not be represented with concentrated mo- 
ments of inertia in FLEXSTAB. (This approximation is not 
believed to have introduced appreciable errors.) 

2. Representation of wing inertia properties was by dumbbell 
masses in the plane of the wing. The ISIC routine does not 
permit modeling the effects of the mass distribution trans- 
verse to the chord plane. (These approximations are not 
believed significant.) 

3. The thin body representation of the wing did not permit fore 
and aft bending. On the other hand, the computed uncoupled 
fore and aft bending frequency of the TIFS wing is 4.83 Hz 
which is within the frequency range of interest for active 
control. It would be expected that a corresponding coupled 
mode of approximately the same frequency would be found by a 
more complete analysis including chordwise bending. However, 
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it is believed there should be little excitation in flight of 
the coupled mode primarily involving wing fore and aft bend- 
ing so its neglect should not have a large influence on the 
computed responses of the TIFS airplane. Some fore and aft 
wing bending exists in all actual coupled modes of the TIFS 
airplane. The importance of this motion could be evaluated 
by comparison of results obtained using the ISIC program of 
FLEXSTAB in computing normal modes with those obtained with 
the NASTRAN program including fore and aft bending. 

4. In the present FLEXSTAB solution, the engine and nacelle 
were treated as a slender body whose inertia properties 
could be represented by concentrated masses along its axis. 
The inertial forces due to the nacelles are not completely 
represented because the rotational inertias of the engine/ 
nacelle are neglected. The addition of engine degrees of 
freedom would have some influence on the vibration modes 
discussed in this report as well as adding modes primarily 
involving engine motion. It will be possible to evaluate 
these effects by comparison with vibration modes for TIFS 
being computed with NASTRAN. 

It is believed that further study will shoitf that the limitations 
of the ISIC program which have been .listed above are not serious in applying 
FLEXSTAB to ACT TIFS. However, the External Structural Influence Coefficient 
(ESIC) option of FLEXSTAB might be used if a more rigorous structural model- 
ing were desired; this program makes it possible to incorporate results from 
NASTRAN into the FLEXSTAB analysis. 

Uncertainties and shortcomings in the basic parameters provide 
another source of error in the vibration analysis results. Some of these 
are listed and discussed below: 
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1) Fineness of wing weight distribution 


Twelve spanwise stations were used in the wing mass 
distribution on Table 1. It is believed that this distribu- 
tion is fine enough for analyzing vibration modes in the fre- 
quency range treated in the present report. Additional panel 
points would be required to obtain a more accurate descrip- 
tion of the third and fourth wing bending modes. 

2 ) Wing elastic axis 

The wing elastic axis assumed in describing TIFS was 
based on an early analysis of the Convair airplane. No ex- 
perimental data is available to verify this assumption. 

The assumed elastic axis could easily be off by 5%. 

3) Engine sprung weight and nacelle inertia properties 

Only limited inertia data were available on the Alli- 
son Prop-Jet conversion of the Convair and the most recent 
data was obtained verbally. Consequently, the engine c.g. 
and moments of inertia are possible sources of error ._ 

4) Fuselage stiffness 

Estimates of fuselage stiffness properties could be 
in error for parts of the fuselage where there are cutouts 
and redundant structure. 

53 Control surface degrees of freedom 

The fact that control surface degrees of freedom are 
not included in the analysis may have resulted in some differ- 
ences between measured and computed nodal lines . However, 
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this is not strictly an error but the result of comparing 
modes' which should not be identical. 

2.3.9 Recommendations for Improving Structural Representation 

Several procedures might be used to refine the structural repre- 
sentation of TIFS developed in this report. Perhaps the simplest, which was 
suggested by NASA personnel at the final oral report, would be to assume that 
the computed mode shapes are correct and merely substitute measured for com- 
puted natural frequencies. 

A first step in a more extensive program to improve the structural 
representation would be an evaluation of the validity of the simplifications 
made in computing vibration modes as discussed previously.. This would' include 
an analytical evaluation of the effects of fore and aft bending of the wing 
and the engine degrees of freedom. 

However, it is more likely that remaining differences between 
measured and computed frequencies and' shapes of the principal structural' vi- 
bration modes are due to errors in basic parameters used in the analysis. 
Improved fundamental parameters might be obtained by more refined analyses or 
by appropriate tests while improvement in modal parameters could also be ob- 
tained directly by parameter identification techniques. 

A sophisticated modern finite element analysis might be made to 
improve structural parameters involving a detailed breakdown of the wing, tail, 
and side force surfaces. A structural model generated with an external pro- 
gram such as NASTRAN or Atlas could be incorporated into FLEXSTAB with the 
ESIC routine and used in subsequent computations of stability and control 
derivatives. This would be a costl-y procedure because of the difficulty of 
assembling the required data and determining required computer inputs. Such 
an approach is not believed justified because classical beam-type analyses 
should be adequate for the comparatively high aspect ratio TIFS configura- 
tion. However, some further analysis is believed worthwhile to check the 
validity of the assumed wing elastic axis. 



As mentioned previously, there are uncertainties concerning the 
stiffness of the TIFS fuselage because of cutout and redundant structure. 
Critical areas of the fuselage could be reanalyzed to obtain new effective 
stiffness estimates. Alternatively vibration analyses could be carried out 
to determine the effect of using different fuselage effective stiffnesses 
within the range of their uncertainty for critical fuselage stations. 

If engine degrees of freedom are found to influence the computed 
vibration modes, further verification should be made of engine sprung weight 
inertia properties. This might be accomplished by a detailed weight and bal- 
ance analysis of the engine-transmission system or by experimental determina- 
tion of the required properties. Furthermore, if it is established that en- 
gine degrees of freedom can make a significant contribution to the response 
of TIFS to turbulence, analyses should be carried out to determine the effect 
of gyroscopic moments on the responses. 

A possible explanation which has been suggested for some of the 
discrepancies between measured and computed mode shapes is that the control 
sui'face degrees of freedom were not included in the analytical solution, but 
may have influenced the vibration test results. This question could be inves- 
tigated by finding the effect of coupling control surface motions with some 
of the computed modes. However, it is believed that structural modes which 
are computed with locked control surfaces will be the most satisfactory for 
use as a basis for analyzing the flight responses of TIFS. Inertial forces 
due to control surface movement can then be treated as forcing functions 
exciting the structural modes. 

Although the currently proposed flight tests of the TIFS airplane 
are concerned with longitudinal response, it is believed important to obtain 
a model for antisymmetrical structural vibrations of TIFS. Atmospheric tur- 
bulence will excite both longitudinal and lateral responses and a knowledge 
of the antisymmetric structural modes will help in separating these effects. 
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Finally, it should be pointed out that the current vibration test 
data for the TIFS airplane does not provide a complete description of the 
structural vibration inodes to -compare with the analytical solutions . A more 
refined vibration test would be required to obtain this information. An 
alternative procedure would be to compare modal response data measured in 
flight with analytical predictions. This procedure is recommended since the 
instrumentation required could also be used in parameter identification to 
determine the actual parameters desired for active control system design. 

2.4 AERODYNAMICS CALCULATIONS AND RESULTS 


2.4.1 Selection of Flight Conditions 

In order to span an adequate flight range, three flight conditions 
were selected for analysis. These conditions roughly corresponded to landing, 
climb and cruise configurations for the TIFS airplane. The climb configura- 
tion was chosen because flight records of the TIFS flying in actual atmos- 
pheric turbulence were earlier obtained at this condition. 

Trim conditions are summarized in Table 12. Also shown in the 
table are the trim <x , and thrust required for trim as predicted by 

FLEXSTAB. A number of important properties of the TIFS used for all flight 
conditions are listed in Table 13. 

The relationship between the three flight conditions and the TIFS 
flight envelope is shown in Figure 26. Also shown in Figure 26 is the direct 
lift flap capability in terms of lift coefficient and acceleration and the 
various structural and aerodynamic limitations to the flight envelope. 

The aircraft weight of 2.42 x lO'’ N and the weight distribution was 
obtained based on the best information available for the TIFS as of 
May 1971, the date near which both ground vibration tests and flutter flight 
tests were conducted. Fuel load of 14,200 N/side corresponds to the fuel 
present during the ground vibration tests and is a typical full load for the TIFS. 
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TABLE 12 


TRIM FLIGHT CONDITIONS ANALYZED 



"LANDING" 

CLIMB 

CRUISE 


M 

.2 

.246 

.456 

FLEXSTAB 

H (meters) 

61. 

1402. 

3048 

INPUTS 

2 

q (Newtons/meter ) 

2819.68 

3629.23 

10155.37 


cc degrees 

7.0227 

4.641 

-.5188 


9 degrees 

7.0227 

4.641 

-.5188 

FLEXSTAB 

Thrust (N ewtons ) 

34,283. 

24,576. 

5973 

OUTPUTS 

Elevator angle (deg) 

-2.517 

-1.400 

+1.477 


Lift Coefficient 

1.005 

.781 

.279 



NOTES : 

1) For all flight conditions 

7=0, trim angular rates and bank angle are zero. 

2) "Landing" is actually a landing approach condition selected for analysis. 
No flap deflections were used for this condition which explains the high 
thrust required to obtain trim. 
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TABLE 13 


MISCELLANEOUS CONSTANTS 


Wing Mean Aerodynamic Chord 2.90 m 

2 

Wing Area 85.47 m 

Wing Span 32.11 m 

Location of LE of Reference Chord (from nose) 12.405 m 

Center of Gravity (Reference Axis System*) 

X = 13.2381 m 

Y = 0 m 

Z - -.3357 m 

C.G . Position 28.7% MAC 

2 

^pitch = 746,900 Kg m = Pitching moment of inertia of total 

airplane about C.G. 

M = 24700 Kg = Total airplane mass 

*The reference axis system is an orthogonal axis system 'with origin 
on the plane of symmetry at the nose of the TIFS. Positive X is 
aft, positive Y is out the right wing and positive Z is up. 
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All analysis in this report is for the TIFS in the clean configur- 
ation, i.e. all controls initially at zero deflection and gear up. Although 
the landing configuration would normally -require some flap deflection; this 
was not included in the analytical model because parameter identification re- 
sults were for essentially clean configurations and because of limitations 
in the possible trim inputs of the FLEXSTAB computer program. The effects of 
flap deflection on TIFS aerodynamics is included in Reference 5. 

2.4.2 Rigid Body Stability and Control Parameters 

Although the primary contribution of FLEXSTAB in aircraft analy- 
sis is in the formulation of the flexible equations of motion as opposed to 
computing rigid body stability and control derivatives, it is interesting to 
compare the FLEXSTAB rigid body or quasistatic results with results obtained 
using other techniques. Two sets of stability and control parameters are 
available. The first representing a combination of wind tunnel, analytical, 
and flight test derived results was published in Reference 5. This docu- 
ment contains the best description of the TIFS derivatives available at that 
time. The second source of stability and control parameters was obtained 
using the Calspan Beysian Maximum Liklihood Parameter Estimation program 
(BML) . This technique determines parameters of a given dynamic mathematical 
model which minimizes the mean square errors between the responses of the 
model and responses of the aircraft measured in flight. 

The results of the parameter estimation methods are subject to the 
following qualification. The elevator,, rudder and aileron control surface 
measurements were not available directly. Instead, the output of each control 
surface servo piston was measured which may differ from the control surface 
response due to effects such as linkage flexibility. Surface position is not a 
direct feedback quantity. The net effect of this is that available measure- 
ment of elevator motion (i.e. the servo output) for the elevator inputs avail- 
able may differ from the actual elevator motion. This may introduce errors in 
identification of the pitching moment derivatives and surface effectiveness. 
Flexibility in the aileron linkage is less significant. 


72 



The FLEXSTAB program yields stability derivatives computed in 
several ways. The first assumes that the aircraft is totally rigid. These 
parameters would he comparable to wind' tunnel test results using a rigid wind 
tunnel model. Another useful form of parameter estimates assumes that the 
aircraft is flexible but that the structural modes (generally at high frequen- 
cies) deflect instantaneously. Using this assumption, each rigid body sta- 
bility and control derivative will be modified by an aeroelastic correction 
term to account for elasticity. 

The nature of these quasistatic corrections are described more 
precisely below. The system equations are given by: 


X 1 - F, f x f *■ X'Z + G-fU. 

Z Z ~ F Z1 X ! + F 2Z Z Z + ^2 46 


where is a vector of rigid body response variables 

is a vector of structural mode response variables 
u is a vector of control inputs 
F a are matrices of stability derivatives 
& f are matrices of control effectiveness derivatives 

The quasistatic assumption requires that = 0. Therefore the following 
algebraic equation results: 


0 - F 2] X j + F 2 2 % 2 & 2 


or z z - F Z 2 F Z1 Zj - F 2 2 u 

Substituting this result into the differential equation for yields 

= \^11 ~ ^12 F ~22 + ~ F >Z ^ 2 ~\ U 

l 

%t ~ ^cp-S LL 
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, where and are the quasistatic stability and control parameter 
matrices . 


Note that the eigenvalues and eigenvectors of -F ^ are not the same as 


those for F 


qs 


Comparison of the various parameter estimates are shown in Table 
14 for the cruise condition. Lift and drag coefficients are the usual wind 
axis quantities and all derivatives are nondimens i onal . The column entitled 
"Best Estimates July 1972" . is from Reference 5 and the "Parameter Identifi- 
cation" column is from the BML results. Both contain quasistatic derivatives 
and can be compared directly only to each other and to the column "FLEXSTAB 
Quasistatic". The "FLEXSTAB, Rigid" column is included to show the contri- 
bution of the quasistatic correction. 


Drag parameters are difficult to predict using linear theories 
which ignore viscous and separated flow effects, are difficult to obtain 
from parameter identification methods using short record lengths, and rela- 
tively easy to obtain in wind tunnel models. As a result, the Cjj and 
from the "Best Estimates July 1972" are recommended for ACT modeling of the 
TIFS, Other drag parameters are not considered significant at this time. 


The estimates of lift curve slope, between the "Best" and para- 
meter identification results agree with each other to within three percent. 

As might be expected, the lift curve slope was overpredicted by FLEXSTAB. 

For the lift due to direct lift flap, , agreement is apparent between 

"Best, 1972" and FLEXSTAB Quasistatic but^not with the parameter identifica- 
tion results. Reference 5 has shown that the direct lift flap effectiveness 

is a nonlinear function of deflection whose most linear region is S„ - -3.5 

o ' 

-5 . The value shown in the "Best, 1972" column and for the FLEXSTAB result 

is linearized about = 0, while the parameter identification results define 
the best linear fit over the flap input range -8 < < 4. These re- 

sults are not directly comparable. To resolve this inconsistency it is 
recommended that inputs be applied to the TIFS direct lift flap in flight 
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TABLE 14 


Notes: 


COMPARISON OF PARAMETER ESTIMATES 
CRUISE CONDITION C.G. = .266 c 




BEST ESTIMATES 
JULY 1972 
(QUASISTATIC) 

BML 

PARAMETER 

IDENTIFICATION 

(QUASISTATIC) 

FLEXSTAB 

ESTIMATES 

(QUASISTATIC) 

c -z>& 

(1/rad) 

.124 

.09 

.185 

c °tf. 

(1/rad) 

0 

0 

.0278 

Cj >s e 

(1/rad) 

.0142 


-.00074 

C *8 % 

(1/rad) 

.026 

.004 

y— — ■ 

C» s 

°sa. 

(1/rad) 




C L 

(1/rad) 

6.02 

5.85 

6.38 


(1/rad) 

0 

0 ■ 

15.7 

< 

(1/rad) 

.63 


.727 


(1/rad) 

1.32 (Note 1) 

.96 

1.28 

C L r 

°$a. 

(1/rad) 

.53 (Note 2) 

.181 (Note 2) 

.595 

c 

m Qi 

(1/rad) 

-.43 

-.494 

-2.15 

C rrt ^ 

(1/rad) 

- 45 

-24.3 

-39.6 

Cm s e 

(1/rad) 

-2.18 

-2.2 

-2.71 

Cm So 

(1/rad) 

-.155 

-.117 

-.347 

r 

(1/rad) 



-.345 


(1/rad) 

-IS. 6 

-8.4 



1) This is a linearization about So - 0 of this nonlinear aero- 
dynamic effectiveness . 

2) Estimated from rolling moment due to asymmetric aileron deflection. 

3) The notation ~ denotes information not available. 
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which limit flap motion to the linear range S ^ = -3.5 i 5° and that the 
parameter estimate resulting from that input be utilized in the linear de- 
sign of the ACT control systems. 


In determining the lift due to symmetric aileron in the "Best" 

and parameter identification , it was assumed that the rolling moment due to 

asymmetric aileron deflection, , could be converted directly from a 

“a. 

moment to a force derivative by the appropriate nondimens ional length. The 
computation is as follows: 


Define the control deflections 


a. 




where and are the deflection of the right and left 

ailerons respectively, positive down. 


The total rolling moment due to asymmetric aileron deflection is given by: 

= c 2 « s? - c, 


C £ (& CL ) - $a 

a. 




R. r> P L 

L S, 


'a. 


The moment on each side of the aircraft is given by: 

jn 0 K 


= c 




a. 


<u = 

“a. 


If it can be assumed that the center of pressure of each aileron is a dis- 
tance X cp from the plane of symmetry, then the total lift force due to 
aileron deflection will be 


C / = C £ Cj_ 

TOTAL R./GHT LEFT 

AILERON A/CERON 



y- cj it 

* CP °a. 
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' CP 


C £ (&* + &$) = C £ - 

* *' X cp ^o.l * / 


2b r x 

xZ ^ 5 * 




a. 


If the ailerons are deflected asymmetrically, then cP, Z = - and =0 

^ ft/ y 

as expected. For the T1FS, the following properties apply: 


Tot a i 




Therefore, 


11.8 meters (38.7 ft) Reference 6 
32.2 meters (105.3 ft) 

-.0975 "Best, 1972" and Convair, Reference 24 
-.0333 Parameter Identification 
high speed cc =" 0° 



+.53 "Best 1972" and Convair 
+.181 Parameter Identification 


In the parameter identification effort, the actual aileron surface 
deflection was not measured. This may explain the discrepancy in rolling and 
lifting effectiveness of the ailerons. However, the measured aileron servo 
output is not believed to be substantially different from the aileron surface 
position. 


Additional experimental analysis performed in 1971 (Reference 11) 
using a combination of simple analytical techniques as well as analog match- 
ing indicate that = -.0418 at cc =5°. This compares favorable with 

the parameter identification result, C_g„ = -.039 at cc =5°, but not with 
the "Best Estimate July 1972", Cap = -.0975 for all CC . 

° CU _ 


Ref. 13. Riedler, A. D. et al: "Aerodynamic Data Flight Criteria and Unit 

Load Distribution for Structural Design of the Model 340 Airplane 
Turboprop Configuration," Consolidated Vultee Aircraft Corp. Report 
No. ZU-340-001, August 1951. NASA CR-158890. 
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This discrepancy in aileron lift effectiveness is considered to 
be perplexing in light of the fact that the TIFS ailerons are identical 
aerodynamically to the Conva-i-r 340’ ailerons. It is recommended that the 
actual aileron lifting and rolling effectiveness be measured in flight again 
with and without side force surfaces before ACT control policies are imple- 
mented. 


Few stability derivatives are more important than the pitching 
moment due to angle of attack, C m . Table 14 shows that the "Best Esti- 
mate July 1972" agrees well with the parameter identification result. How- 
ever, FLEXSTAB predicts a substantially more stable aircraft. The exact 
cause of this discrepancy is not known at this time. 

The pitch damping derivatives, C and C are in reason- 

able agreement between "Best" and FLEXSTAB. The parameter identification 
result indicates lower damping, but this is believed to be the result of 
phase errors in the measurements of the elevator due t'o flexibility in the 
actuating mechanism. It is recommended that these derivatives be verified 
by parameter identification using the elevator position sensor rather than 
the elevator servo output. 

The pitching effectiveness of the elevator, , was noted 

°e 

in the parameter identification effort to be a function of dynamic pressure. 
It is believed that this resulted due to aeroelastic effects on the elevator 
actuation mechanism. However, the extrapolation to zero dynamic pressure 
of the parameter identification results agreed reasonably well with the 
"Best' Estimate 1972". FLEXSTAB overpredicted the control effectiveness by 
23%. This overprediction is not unusual since flow over trailing edge de- 
vices is frequently partially separated, an effect' for which the theory has 
not accounted. 

It should be recalled that FLEXSTAB was not expected to produce 
as accurate rigid body stability and control derivatives as available from 
other sources. The proper use of the FLEXSTAB system would involve using the 
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AIC correction program to modify the lifting capability of panels whose ef- 
fectiveness was not correctly modeled. An example might be decreasing the 
lifting capability of the tail panels in order to more closely match and 

, Certainly all control surfaces should be so corrected. In this 
°e 

way the generalized forcing terms on the structural mode equations would be 
more accurate as well. 

Comparison of quasistatic derivatives at the landing and climb 
configuration are shown in Tables 15 and 16, respectively . 

2.4.3 Flexible Aircraft Equations of Motion 

The flexible equations of motion of the TIFS aircraft have been 
computed at two flight conditions, climb and cruise, for symmetric motion 
only. Time prohibited computation at the third flight condition. 

The form of the equations of motion is. as follows: 

Ax -h B x + C% ~ 0 

where X T = \_cl, -txr , , 9, ^ •, q z , q 3 7 , q s , , f 7 ] 

A, B, C are square matrices of dimension 11 x 11 
u = component of inertial velocity in X direction, m/sec 
w = component of inertial velocity in Z direction, m/sec 
q = pitch rate, rad/sec 
0 = pitch attitude, rad 

q ■ are structural mode deflections of the aircraft 

The seven structural modes are those modes discussed earlier in 
Section 2.3 of this report in order of increasing vibration frequency. Sen- 
sor equations have not yet been formulated. The equations of motion were 
computed using the Residual Elastic formulation in which the differential 
equations are corrected for the static effects of higher order structural 
modes above the seventh mode. The coefficients (A, B, C) to the above equa- 
tions are included in Appendix B. 
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TABLE 15 


COMPARISON OF PARAMETER ESTIMATES 
LAND-ING CONDITION, C.G. ='.266 c 




BML 



BEST ESTIMATE 

PARAMETER 

FLEXSTAB 


JULY 1972 

IDENTIFICATION 

ESTIMATES 


QUASISTATIC 

(QUASISTATIC) 

CQUASISTATIC) 

C »CL 

.447 

.69. 

1.15 

S 

0 

0 

1.33 

Cj>s e 

.0511 


.0832 

c L 

u oc 

6.02 

5.05 

5.78 

% 

0 

0 

15.20 

C **e 

.63 


.805 

C m 

m a. 

-1.33 

-1.05 

-1.97 

C/rip 

-45.0 

-24.3 

-39.1 

C^r. 

-2.18 

-2.2* 

-3.09 

m s e 




• 

cc 

-15.6 

-8.4 



Extrapolated to- q = 0. 
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TABLE 16 


COMPARISON OF PARAMETER ESTIMATES 
CLIMB CONDITION, C.G. = .266 c 



BEST ESTIMATE ' 
JULY 1972 
(QUASISTATIC) 

BML 

PARAMETER 
IDENTIFICATION 
(QUASISTATIC) 
(NOTE 1) 

FLEXSTAB 
ESTIMATES 
.(QUASI STATIC) 

Cj, 

i 

.348 

+ .43 

.854 

C H 

0 

o 

.919 


.0397 


.0559 

c L 

6.02 

5.3 

5.88 

u cc 




C, 

0 

0 

15.3 

9 




° L ^ 

.63 


.779 

c m 

-.93 

-.68 

-1.98 

cc 





-45.0 

-24.3 

-39.0 

0/n s e 

-2.18 

-2.2 (Note 2) 

-3.05 

c m . 

-15.6 

-8.4 



NOTES: 

1) These results are interpolated from data at other flight conditions. 
2} Extrapolated to q = 0. 
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The set of stability and control derivatives of the rigid body 
degrees of freedom corresponding to the above equations of motion are also 
shown in Appendix B. The column labeled "elastic increment" represents the 
static elastic correction for the structural modes not modeled dynamically 
in the equations of motion. 

For the cruise condition, the following rigid body stability der- 
ivatives were chosen for subsequent calculations. It is felt that these values 
are the most accurate for the TIFS airplane. 

C/ = 5.92/rad 

^ or _ 

C = -.26/rad about c.g.= .266 c 

m cc 

These changes are reflected in the cruise equations and stability derivatives 
in Appendix B. 

An assessment of the accuracy of these equations must be post- 
poned until sensor equations from FLEXSTAB can be computed. At that time 
comparisons of computed and observed responses of the TIFS in turbulence or 
responses to specific control inputs should be compared. These comparisons 
will be performed before control loops are closed around the TIFS aircraft. 

The flexible equations of motion computed to date have included 
symmetric degrees of freedom only. An analysis of the importance of anti- 
symmetric degrees of freedom on the overall objectives of this contract will 
be undertaken. 

2.4.4 Aerodynamic Hinge Moments 

Computation of hinge moments via theoretical means was not pro- 
ductive for the TIFS because of the presence of aerodynamic balance and tabs 
on the various control surfaces. The ailerons, for instance, are aerodyna- 
mically balanced using a Frise balance method. The resulting hinge moments 
are nonlinear functions of a. and . Estimated aileron hinge moments 
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obtained by Convair (Reference 13) are shown in Figure 27 and reveal that for 
small cc the hinge moments are not dependent on oc , for small 8 the hinge 
moments vary linearly with 8 ^ and for large 8 and cc the functional depen- 
dence is nonlinear. Trim tabs on both ailerons are actuated by two sources. 

A trim knob in the cockpit allows the pilot to change the trim tab position 
b y - 7°. In addition, the trim tab is slaved to the aileron position for hinge 
moment control. The combined relationship between trim tab deflection 8 ^ and 
aileron position is shown in Figure 28. Hinge moments about the aileron hinge 
axis can be computed as shown below. 


HM = 

where S ^ 


^ a. > <x ) 

% 




«-s 


tt 


// M 


9 S CL C *. 




= area of one aileron aft of hinge line (^) = 
= average aileron chord aft of - 


J21.91 ft 2 
[2.035 m 2 
[1.26 ft 
[.384 m 


= Cfj from Figure 4 
= dynamic pressure 

= -.0028 per degree = trim tab effectiveness 
= trim tab deflection, is a function of £ and pilot's 
trim command. Figure 5 

= hinge moment due to deflection of one aileron 


Elevator aerodynamic hinge moments are shaped by two servo tabs, 
one on each tailplane. Deflections of each tab are determined by elevator 
deflection in a nonlinear relationship (tab deflection measured with respect 
to the elevator is down for both positive and negative elevator deflections). 
Gearing on the left hand tab, called the control tab, is not the same as on 
the right hand tab, called the trim tab. The trim tab has the additional 
function of trimming control wheel forces, an action which is initiated via a 
trim wheel in the cockpit. Total control column forces experienced by the 
pilot are a function not only of the aerodynamic hinge moments, but also the 
actions of a bungee spring (down spring) and a bobweight (control yoke bal- 
ance weight) . 
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Figure 27 ESTIMATED AILERON HINGE MOMENTS 
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The aerodynamic hinge moment about the elevator hinge line can be 
computed as follows: 


where 


HM = c e 




Cm. 






tt 


e s. 




cc, + C, 


■O' 


C 




e S 


+ C H, 




G S, 


ct 


*«] 


S e = 


* 

e at 

4*., = 


. ct 
HM 


elevator area aft of tjf 
RMS chord of elevator aft of $ 
elevator hinge moment due to <x 
elevator hinge moment due to 6 e 
trim tab effectiveness 
'control tab effectiveness 
hinge moment due to combined action of both left and 
right elevator surfaces 


= 71.32 ft 
= 2.28 ft 
= -,002/deg 
= -.0033/deg 
= -.0051/deg 
* -.0037/deg 


Trim tab and control tab positions are not currently measured on 
the TIFS. However, for each control surface the hinge moment can be computed 
from the differential pressure at the surface actuator. 


Hinge moments for the direct lift flap which, like the aileron, 
has a nose seal, have been estimated from flight test measurements of differ- 
ential pressure at the DLF actuator. Data was obtained for various flap de- 
flections and angles of attack and was reported in Reference 14 , A cross 
plot of hinge moment versus DLF deflection for constant angle of attack, shown 
in Figure 29 , reveals that the hinge moment coefficient is nearly a linear 
function of DLF deflection. A linear approximation to this family of curves 
leads to the following representation of the DLF hinge moment. 

HM = 

where ^ = dynamic pressure 

S = area of one DLF, 4.161 m 2 (44.79 ft 2 ) ' 

0 

Ref. 14. Flight Research Staff: "Development, Design and Fabrication of the 

Total In-Flight Simulator (TIFS) . Air Force Flight Dynamics Labora? 
tory Report AFFDL-TR-71-77, August 1971. 


v/ & 


+ C 
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£ = average chord of DLF, .677 i, (2.22 ft) 

Cu = -.094/deg 

cc 

C H - -.0075/ deg 

O ^ 

- hinge moment due to deflection of one flap 

Aerodynamic hinge moment data for the aileron and elevator were 
obtained for the Convair 340 based on both wind tunnel testing and flight 
testing of the Convair 240. Because the 340 and TIFS are very similar air- 
craft, the hinge moment data so derived are considered to be representative 
of the actual TIFS hinge moment data. Direct assessment of accuracy is not 
possible without more thorough study of Convair 's techniques. More accurate 
hinge moment data might be obtained by direct measurement in flight using 
differential pressure across the surface actuator. 

Computations of the direct lift flap hinge moments were based on 
the following relationship: 

c _ *j) 


where AP^ 


S \ C Tr 


is the actuator pressure 

product of actuator piston. area and effective arm * 
and q are respectively flap area, flap mean chord and 
dynamic pressure. 


The primary sources of error in this computation will be in meas- 
uring dynamic pressure and actuator pressure. Estimates of the accuracy of q 
and AP^ are ± 5% and ± 8% respectively. For the worst case, the resulting 
accuracy for C ^ would be of the order of t 14%. Better accuracy could be 
obtained by more accurate calibration of both sensors. 
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2.4.S Pressure Distribution Data 

Two sources of pressure or load distribution data are available. 

The first was utilized by Convair in the basic design of the Convair 340 series 
aircraft. Reference 13. Except for the effects of flaps and side force sur- 
faces, this data is considered to be representative of the TIFS. The second 
source of data is FLEXSTAB which computes the pressure coefficient for numer- 
ous points on the aircraft for the trim configuration of the loaded, flexible 
aircraft. In producing the pressure distribution data, the pressures on the 
upper and lower lifting surfaces of the aircraft are differenced and nondi- 
mensionalized by the appropriate constant, i.e. 

p - P 

r Locozr u.pper 


If C-p.. and S:- are the pressure coefficient and area of the aero- 
dynamic panel in the ith spanwise row and jth chordwise row of the wing, then 


the lift load for the j chordwise row is given by 



A/C 

Z 

7T, H 



where is the lift on the ith chordwise row 
q is dynamic pressure 

AY. is the width of the ith chordwise row. 
1 


The lift load distribution obtained above can be compared ‘directly 
to the lift load distribution computed by Convair in the following manner. 

, c(y ) e £ (y) - c( S ) CjJy) C e * cCj k 

where dL'(tj) is the lift per unit span 

c(y) is the local wing chord 
C^(y) is the local lift coefficient 
C£ (y) is the additional lift distribution 
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Co is the basic lift distribution, the lift distribution 
-*-b 

when C, = 0 

C L i-s wing lift coefficient. 

A comparison of lift load distribution computed by FLEXSTAB and Convair is shown 
in Figure 3,0 . Very reasonable agreement is apparent except at the wing tip 
where an additional chordwise row in the FLEXSTAB model is needed. From this, 
it is reasonable to conclude that FLEXSTAB and Convair will also predict 
approximately equal wing root bending moments. 

The FLEXSTAB pressure coefficient data corresponding to this condi- 
tion is included in Appendix C. Convair data in the form of additional and 
basic lift load distributions are shown in Figures 31 .and 32 . The additional 
load distribution due to aileron deflection is not currently available from 
FLEXSTAB, however, the Convair estimate is shown in Figure 33. 

The moment distribution corresponding to the above-mentioned pres- 
sure distributions are not yet available from FLEXSTAB though they can be com- 
puted from the data of Appendix C. The Convair estimated moment distribution 
about the wing elastic axis (36.5% ) is shown in Figure 34,. 

Chordwise pressure distributions as computed by FLEXSTAB at three 
locations are shown in Figures 35 and' 36 . These figures represent pressure 
coefficient distributions at the trim configuration with flap and ailerons un- 
deflected. Pressure distribution plots provided by Convair based upon wind 
tunnel tests of a Convair 240 model corrected to the 340 configuration, are 
shown in Figures 37 and 38. Pressure distributions for two aileron deflec- 
tions are presented in Figures 39 and 40. The total pressure coefficient 
is computed from these figures from the following: 


AP I 


AP 


I C4+0 


\c_l 4o 


Wir 


^ / addit z-ona. 
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cCj or L-j/qAYi, meters 



! y, DISTANCE FROM PLANE OF SYMMETRY, meters 


j; Figure 30 TOTAL WING LIFT LOAD DISTRIBUTION 


ADDITIONAL LIFT, cZn , meters 



Figure 31 ADDITIONAL LIFT LOAD DISTRIBUTION 





BASIC LIFT, cC , , meters 


0.61k 2 






cC^ scale is adjusted to produce 
aCj = 1.00 with just one aileron 
deflected for positive C t 

1 • e ■ > Cl T0TAL = Cl LEFT AIL + Cl RIGHT AIL 

Load on right wing is equal to 

C? RIGHT x Th1s Distribut ' ion 

Load on left wing is equal to 
Cz , FFT x Negative of This 
LtM Distribution 


% SEMISPAN 


CP Position 


Figure 33 WING LOAD DISTRIBUTION DUE TO AILERON DEFLECTION 


CP POSITION, % CHORD 








.Moment Coefficient about 0.365 











Figure 35 PRESSURE COEFFICIENT AT TRIM 

cv FRL = 7.02°, LANDING, M = 0.2, H - 61 meters 
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ACp at HORIZONTAL TAIL ACp at DIRECT LIFT FLAP 4Cp at AILERON 




Figure 36 PRESSURE COEFFICIENT AT TRIM 

*TRL = -0.52°, CRUISE, M = 0.456, H = 3048 meters 


9.7 




Ap/q, FOR ADDITIONAL DISTRIBUTION 










AP/q s FOR. ADDITIONAL DISTRIBUTION 






















Ap/q 















where ~ section lift coefficient 

AP = local static pressure minus free stream static pressure 
^ = -dynamic pressure 

Direct comparison of FLEXSTAB and Convair pressure distributions 
have not as yet been completed. 

The accuracy of the pressure-load data presented in this section 
is not known at this time. A study of the accuracy of the Convair techniques 
was not made. FLEXSTAB pressure distributions, on the other hand, are subject 
to the following qualifications. Comparison between wind tunnel pressure data 
and FLEXSTAB predictions for the B-52 aircraft. Reference 7, Vol. IV, has 
indicated that in the region of the leading edge ( x/c <.05) FLEXSTAB substan- 
tially overpredicts the differential pressure. The excessive peak near the 
leading edge is caused by the leading-edge aerodynamic singularity inherent 
in the aerodynamic theory. In the TIFS model the leading edge panel was approx- 
1 imately 10% of • the local chord, much larger than for the B-52 model. Since 
one pressure is computed for each panel the effect of the larger panel was one 
of averaging the theoretical pressure distribution over the width of the panel 
thereby reducing the effect of the leading edge singularity. In addition, 
FLEXSTAB will not account for separated flow phenomena which may exist at cer- 
tain flight conditions (e.g., in landing where = 7°, Roor = H°)- 

In regions of separated flow FLEXSTAB will also overpredict the pressure coef- 
ficient. A quantitative assessment of the accuracy of pressure distribution 
data is- not available at this time. 
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Section 3 


SENSOR AND SERVO SYSTEM 


3.1 INTRODUCTION 


This section briefly describes the TIFS servo and sensor concepts, 
installation and general arrangement. Dynamic and static performance informa- 
tion for the servos and the instrumentation are described in more detail, because 
these are the characteristics of the system pertinent to an ACT control system 
configuration for the aircraft. 

3.2 GENERAL INSTRUMENTATION AND SERVO SUBSYSTEM DESCRIPTION 

The electronic system for TIFS is composed of the subsystems ar- 
ranged in Figure 41 and listed as follows: 

1. Sensors and Sensor Electronics — The sensors measure various 
angles, rates and accelerations of the TIFS aircraft. These 
sensor outputs are converted to properly scaled DC voltages 
representing parameters used in the control of the TIFS. 

2. Feel System — A feel system for each evaluation pilot provides 
variable force versus position gradients for the elevator and 
aileron control wheel and the rudder pedals. This control 
system can simulate a linear spring feel or nonlinear feel 
characteristics composed of variable amounts of deadband, 
breakout force and hysteresis. 

3. Surface and Throttle Servos — Electrohydraulic servos con- 
trol the position of the elevator, aileron, rudder, flaps and 
side force surfaces. Magnetic powder clutch servos are used 
for the throttle position servos. Rate and acceleration sta- 
bilization techniques are provided to shape the dynamic responses 
of these servos . 
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Figure 41 TIFS ELECTRONIC SYSTEM BLOCK DIAGRAM 











4. Model -Foil owing - Response-Feedback Patch Panel Electronics — 
The MF-RF patch panel ties together all of the other subsys- 
tems of the TIFS electronics. The most important function 

of the electronics is to process the signals between the TIFS 
input variables and airplane variables to provide the required 
surface and throttle position commands. 

5. Model Computer — The model computer provides command augmen- 
tation between the pilot control stick input commands and the 
input signals to the control surface servos. This computer 
can be programmed for as simple as a first order lag to a multi 
controller dynamic system representing the equations of motion 
of a complete aircraft for model -following operation. 

6. Pilot Instruments and Controls — Servo driven flight and 
engine instruments are provided for each evaluation pilot. 

The inputs to each instrument can 'be patched on the MF-RF panel 
from the sensors or the model computer. Throttle, wing-sweep, 
flap, and speed-brake levers are mounted on the center control 
console. Each provides a proportional electrical signal for 
simulated control inputs. 

7. Test Engineers' Control and Instrument Panels - Two test 
engineers control and monitor the operation and performance 
of the TIFS system on separate panels. One panel contains 
switches and status lights that control and indicate the en- 
gagement sequence of the various subsystems. The other panel 
contains instruments and a strip chart recorder to monitor 
the model-following performance. 

8. Digital Tape Recording System — Fifty-eight channels of digi- 
tal recording are provided with tape coding compatible with 
the IBM 360 computer system. A ground playback unit is pro- 
vided which can play back any 8 channels simultaneously on 

an analog recorder. 
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3 . 3 SURFACE SERVOS 


. _ A simplified -block -diagram orf ' a typical surface servo is presented 
in Figure 42. The surface position feedback signal is subtracted from the 
command signal, with the resultant error signal used to current drive the hy- 
draulic flow control valve. High frequency position feedback for the sys- 
tem is obtained from the strut position while low frequency and DC feedback 
are .obtained from the surface position. The surface rate and acceleration 
terms are available to shape the dynamic response of the servo. 

An additional feature of the servo loop is the auto balance system 
which electronically nulls surface command signals before system engagement 
to avoid undesirable transients. 

An electronic safety trip system is built into -the patch panel 
electronics. If a- predetermined maximum safe surface rate is exceeded or if 
a system failure causes an abrupt high level surface rate command, the system 
is automatically disconnected. Aiso available on the panel are the outputs 
of electrical differential pressure gages mounted across each actuator. Any 
combination of these signals may be monitored by the electronic safety trip 
system. 


Pressure limiting check valves are coupled around each hydraulic 
actuator to prevent damage to the surfaces due to excessive surface hinge 
moments. In' addition to these precautions, either safety pilot may disengage 
the system at any time by using a dump button if he feels that excessive con- 
trol rates or vibrations are occurring. 

To ensure that the flight control surfaces can be easily operated 
by the safety pilot when the TIFS system is "off", several safeguards are 
provided but are not discussed in this report. 
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Figure 42 SIMPLIFIED DIAGRAM OF TYPICAL SURFACE POSITION SERVO 



AIRCRAFT 
DYNAMIC 
RESPONSE . 


10.7 










3.3.1 Surface Servo Performance 


Intuitively, it can be reasoned that the best way to stabilize 
the highly resonant loads of the TIFS surface servos would be to use a notch 
or antiresonant network. At various times it has been proposed to insert 
passive networks inside a servo loop to accomplish this cancellation of sys- 
tem dynamics. The idea is fine on paper; however, in practice this approach 
presents problems. Once a configuration of passive circuit elements is tuned 
to give the required antiresonant characteristic, system operation relies on 
an unchanging load resonance and an unchanging compensation characteristic. 

Since both these effects are very powerful, any drift of load resonance or 
compensation circuitry will prevent cancellation of the resonance or compen- 
sation circuitry will prevent cancellation -of the resonance and may possibly re- 
sult in system instability. 

However, it can be shorn that pressure feedback in a resonant- 
loaded position servo is a form of antiresonant feedback which is continually 
slaved to the exact load resonance phenomena. For this reason, all of the 
surface servos incorporate pressure feedback. The block diagrams for the 
surface servos (Figure -42) show that the pressure signal is not connected 
directly to the servo summing junction but is first passed through a band- 
pass filter. -Without this filter static position errors would- be produced 
due to the force being generated -by the- actuator. 

A 

r ' t, ' 

The primary. position, feedback signal for the surface servos is 
shown in the block diagrams to come, from the actuator for the elevator, 

SFS and DLF servos, whereas the aileron and rudder servos use information 
near the surface being controlled. Surface position feedback was used in 
these servos to remove static position errors due to flexure in the links 
between the actuators and the surfaces. Linkage flexure is a more signifi- 
cant problem for the primary surfaces on the airplane because the actuators 
for these surfaces are connected to the original torque tubes and cables 
which are more flexible than the linkages in the side force surface and 
direct lift flap. 
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Surface position feedback was tried on the elevator servo as 
well as the aileron and rudder servos. However, the low natural frequency 
quency of the elevators produced significantly lower frequency response for 
this servo when surface position feedback was used. Consequently, it was 
decided to use strut position feedback for the elevator servo even though 
positional errors of 5 - .10% occur as a function of aerodynamic, loads on the 
surface . 


In addition to the primary feedback signals shown in Figure 42, 
the TIFS electronic system provides signals proportional to surface rate and 
surface acceleration which may be used to improve the response of the servo. 
Although rate feedback should theoretically provide damping, the amount of 
improvement which could be obtained when it was used in the TIFS servos was 
minimal . 


The block diagrams also show that when surface position feedback 
is used , it is low-pass filtered and the high frequency information is ob- 
tained from the strut position transducer. This technique makes it possible 
to obtain wide-bandwidth servos with good surface positioning accuracy. If 
surface position were used alone with the low-pass filter removed, stability 
criteria would necessitate lower gains with poorer dynamic performance. 

The frequency response of the TIFS surface servos is shown in 
Figures 43 through 50 . The in-flight servo surface frequency response is 
shown in Figure 51 through 55. On some of the surface servos, particularly 
the elevator servo, the effect of aerodynamic loading can be seen to have 
the effect of introducing a first order lag that will be a function of flight 
condition. At maximum dynamic pressure, this lag will reduce the bandwidth 
of the elevator surface servo down to about 1.5 Hz. This effect has been 
verified by in-flight experimentation. 


109 . 



A/R, Amplitude Ratio 

















Phase Angle, deg 



■ Figure 44 TIPS ELEVATOR FREQUENCY RESPONSE 



(A/R), Amplitude Ratio 


10.0 



Figure 45 TIFS AILERON FREQUENCY RESPONSE 



Phase Angle, deg 



Figure 46 TIFS AILERON FREQUENCY RESPONSE 



(A/R), Amplitude Ratio 










Phase Angle, deg 











(A/R), Amplitude Ratio 



Figure 49 TIFS DLF FREQUENCY RESPONSE 


Phase Angle, deg 



: Figure 50 TIFS DLF FREQUENCY RESPONSE 



AMPLITUDE RATIO 



AMPLITUDE RATIO 



Figure 52 AILERON ANGLE RESPONSE TO SERVO ELECTRICAL COMMAND 


119 


PHASE - DEGREES 



# FREQl 


Figure 53 RUDDER ANGLE 
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Figure 54 SIDE FORCE SURFACE ANGLE RESPONSE TO 
SERVO ELECTRICAL COMMAND 


PHASE - DEGREES 


AMPLITUDE RATIO 



FREQUENCY - HERTZ 


Figure 55 DIRECT LIFT FLAP ANGLE RESPONSE TO SERVO ELECTRICAL COMMAND 
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3.4 SENSORS AND AIR DATA SYSTEM 


This section presents pertinent details of the sensors, sensor 
signal conditioning and computed signals which are available in the TIFS 
aircraft . 


Scaled DC signals are trunked to the computer patch panel for 
use in the TIFS flight control system and for recording. These signals are 
arranged into three groups. Group one consists of directly measured parameters 
such as gyro and accelerometer signals. The signals are conditioned-demodulated 
if necessary, filtered, and scaled. Group two consists of air data type sig- 
nals. These signals are derived from air data sensors and the inertial measure-’ 
ments. In group three are the computed signals for inertial angle of attack, 
sideslip, rate of change of angle of attack, sideslip and gust. The features 
and capability of the onboard digital tape recording system are described. The 
physical locations of the basic sensors in the TIFS are identified in Figure 56. 
Distances of the C.G. and sensor locations from fixed aircraft reference points 
may be obtained from Table 17. 

3.4.1 Group One.. Directly Measured Parameters 

Attitude and Heading Channels 


The pitch, roll, and heading angle information is provided by syn- 
chro outputs from the C.G. located attitude and heading reference system (AHRS). 
The specifications for the AHRS are presented in Table 18. The AHRS signals 
are converted to sine and cosine terms, of pitch, roll and heading angle by con- 
trol transformers. The signals are scaled, amplified and then converted to 
DC by phase sensitive demodulators. The demodulated signals are filtered to 
remove the 800 Hz carrier frequency. The outputs of these channels are sine 
and cosine of pitch, roll, and heading angles, pitch angle change from pertur- 
bation engage and roll angle. 
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TABLE 17 


TIFS 

SENSOR LOCATIONS 




Vertical 

Longitudinal 

Lateral 

Sensor 

Location 

Location 

Location 

(In. From Waterline) (1) 

(in.) (2) 

(In. From Butt 

Radar altimeter antennas 

0.0 

654 (ave) 

0.0 

Attitude § Hdg. Gyro 

+35.5 

351 

4 LBL 

Angular Rate Gyro 

+44.5 

370 

3 LBL 

Angular Accelerometer 

+37 

369 

4 RBL 

Linear Accelerometer (eg) 

+37 

369 

1 LBL 

Linear Accel. (Sim. Cockpit) 

+ 3 

- 36 

- 0 

OC Vane 

+ 19 

- 16 

46 LBL 

Vane 

-11 

- 16 

0 

Static Pressure Source 

+ 10 

47 

43 LBL 

Dynamic Pressure Probe 

+65 

+ 78 

45 LBL 

TIFS c.g. Range 

« +36 

370-379 

0 


(1) Waterline — Horizontal plane tangent to lower fuselage skin. 

Measurement above waterline is positive. 

(2) Fuselage Ref. — Extreme forward point of original aircraft nose. 

Distances aft of point are positive. 

(3) Buttock Line — Vertical plane disecting aircraft fore-to-aft. 

, Distances are measured to left (LBL) or right (RBL) 
of buttock line. 

All distances measured in inches. 



TABLE 18 


ATTITUDE AND HEADING REFERENCE. SYSTEM 

Lear Siegler, Inc. 

AN/ASN-50 
Model 5010R 


Bench Vertical Error 
Dynamic Vertical Error 
Free Vertical Drift 
Erection Rate 
Erection Gradient 
Roll Follow-up Rate 


t .25 deg max 
1 1 deg max 
- 10 deg/hr max 
1.3 + .5 deg/min 
2.6 ± 1.0 deg/min/deg 
360°/sec 


Magnetic Heading (static) 
Heading Follow-up Rate 
Heading Slew Rate 
Heading Drift (corrected for 
Earth’s Rate and Gyro Bias) 
Magnetic Slaving Rate 

Magnetic Slaving Gradient 


+0.3 deg 
180 deg/sec 
50 deg/sec 
+ .3 deg/hr RMS 

1.0 + .5 deg/min 
- .3 

1.5 t .3 deg/min/deg 


Features 


Two-Gyro Platform with Servoed Outer Roll Gimbal . 
All-Attitude-No Gimbal Lock. 

Vertically Stabilized Directional Gyro. 

No Roll and Pitch Coupling. 


Functions 


Provides Gyro Roll and Pitch Data. 

Provides Three Modes of Heading 
Gyro Stabilized Compass. 

Free Gyro. 

Compass (no gyro). 

For additional information, reference; 

I 

NAVAIR 05-35-LAA-l 

Handbook, Operation and Maintenance 

Attitude and Heading Reference System AN/ASN-50 
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Three Axis Angular Rates 


The pitch, roll and yaw rate gyros sense body axis angular rates and 
are mounted mutually orthogonal in a single package located at the TIPS c.g. 

The performance specifications for the rate gyro package are given in Table 19. 
The AC outputs of the package are amplified, scaled and demodulated. The 
signals are then filtered to remove the 800 Hz carrier frequency. Additional 
first order filters are installed in each rate gyro channel. A notch filter 
is also provided in the roll rate channel. The outputs of these channels are 
pitch, roll and yaw rates. 


TABLE 19 

3 -AX IS RATE GYRO 


United States Time Corp. 
Part #400750 


jP Rated sensor output t 125 deg/sec = 4V DC 
Threshold ,005 deg/sec 
Resolution <.01 deg/sec 
Ave. Hysteresis .07 deg/sec 
Linearity .5% 

= 49 Hz 

r = .6 

o Rated sensor output + 40 deg/sec = 4V DC 
Threshold .005 deg/sec 
Resolution < ,01 deg/sec 
Ave. Hysteresis .080 deg/sec 
Linearity , 5% 
oO n = 28 Hz 
i = .7 

■f Rated sensor output ±.j40 deg/sec = 4V DC 
Threshold .006 deg/sec 
Resolution < .01 deg/sec 
Ave. Hysteresis .104 deg/sec 
Linearity .5% 

U> n = 29 Hz 

(T " - 7 
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Three-Axis Linear Accelerations 


The three axis linear accelerometer package is mounted directly 
to the airframe at the TIFS C.G., These sensors measure acceleration in the 
aircraft x, y and z axes. The DC output signals are scaled prior to low pass 
and notch filtering to remove high frequency noise. The performance specifica- 
tions for the linear accelerometer sensors are presented in Table 20. The 
channel outputs are n x , n y , n^ , and (n^-1) at the TIFS c.g.. 

TABLE 20 

3-AXIS LINEAR ACCELEROMETER 
C.G. Location 

Unico Controls, Inc. 

Model MCA-5-6X3 


n. 


X 


± .5g = ± 5V DC 
Linearity .1% 
Hysteresis < .2% 
Output Noise 1 mv 
= 40 Hz 
= .8 


*>n 


Threshold 


10" 4 g 


n. 


n. 


± lg = ± SV DC 
Linearity = .1% 
Hysteresis < .002V 
- 38 Hz 
'f = .8 


± 2.5g = ± 5V DC 
Linearity .1% 
Hysteresis < . 003V 
Output Noise — 


a) 


n 


= 4S Hz 

= .8 


Threshold 


Threshold 


10“ 4 g 


10" 4 g 


These acceleration units mounted mutually orthogonal on a common 


base. 


Axis alignment to base ± 


1 


o 
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Three Axis Angular Accelerations 


The pitch, roll and yaw angular accelerometers are mounted mutu- 
ally orthogonal in a package located at the TIFS C.G. These body axis angular 
acceleration signals are both low passed and notch filtered. These accelera- 
tion signals are available for recording but the hysteresis has increased to ' 
a level that negates their use for stability and control applications. The 
original performance characteristics of the angular accelerometers are pre- 
sented in Table 21. The outputs of these channels are pitch, roll, and yaw 
angular acceleration. 


TABLE 21 

3-AXIS ANGULAR ACCELEROMETER 

Unico Controls, Inc. 
Model MCA-6-8G 


Full Scale Output i 500 deg/sec 2 = i. 8.28V DC 

2 

Threshold .01 deg/sec 
Dynamic Linearity . 1% 

Hysteresis < 1% 

iO n = 19.4 Hz 

r = * 7 

. 2 

• Full Scale Output ± 200 deg/sec = ± 3.48V DC 

Threshold .01 deg/sec 2 
Dynamic Linearity . 1% 

Hysteresis '< 1% 

= 18.3 Hz 
r = .7 

r Full Scale Output ± 200 deg/sec 2 = 2..9V DC 

Threshold .01 deg/sec 2 
Dynamic Linearity . 1% 

Hysteresis < 1% 

uJ n = 17 Hz 

r = - 7 

These angular accelerometers have been found to be unacceptable for 
stability and control feedback usage. 
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Two-Axis Linear Acceleration Package 


This two-axis package is mounted between the pilot and copilot 
seats in the simulation cockpit. These body axis acceleration signals are 
scaled and low pass .filtered. Performance characteristics of these sensors are 
presented in Table 22. The outputs of these channels are simulation pilot loca- 


tion n. 


y-p 






and 


- 1 ’ • 


TABLE 22 

2 -AX IS LINEAR ACCELEROMETER 
(Simulation Pilot Station) 

Uni co Controls, Inc. 
Model MCA- 5 -6X2 


rt ± lg - + 5V DC 

^ Threshold 10 g 
Linearity t .1% 

Hysteresis @ zero g = .003V 
<4, = 35 Hz 

f = .7 

± 2.5g = t 5V DC 

Linearity + . 1% 

Hysteresis @ lg = .002V 
= 35 Hz 

f = -7 

Both sensors mounted on common base. Axis alignment to base i 1°. 


Radar Altitude 


The radar altimeter DC output measures distance from the TIFS to 
the closest object within an angular cone from the antennas. The transmit- 
ting and receiving antennas are flush mounted units located in a horizontal 
plane on the bottom of the TIFS fuselage just aft of the cargo door. The out- 
put signal is scaled and differentiated to obtain instantaneous vertical speed 
information. The performance specifications for the radar altimeter are pre- 
sented in Table 23 . The channel outputs are radar altitude and rate of change 
of radar altitude. 
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TABLE 23 


RADAR ALTIMETER 


Honeywell Aerospace Division 
AN/APN-171 ( V) 

HG 9000 Series 


Range 0 to 1000 feet 0 to -40V DC 
RF Frequency 4.3 GHz, 100 Watts peak power 
Pulse repetition frequency 10 KHz 
Pulse Width 35 nanoseconds low altitude 
125 nanoseconds high altitude 

System transient response: First order with .1 second time 

constant 

Altitude Accuracy — Standard Conditions 

Assume present altitude. = x (feet} 
and rate of change of altitude = x feet/sec 
Accuracy = - [ 1.5 + .01 (x) +..05 (x)] feet 

Altitude Accuracy — Service Conditions 

(worst case with component value drift not compensated for) 

Accuracy = t [ 5 + .03 (x),+ .05 (x) ] feet' 

This statement applies to dynamic range subject to the 
limitations of a first-order system with a transient 
response time constant of .1 seconds. 


Vibration and Flutter Accelerometers^ 

These accelerometers are mounted oh the wing tips, horizontal and 
vertical stabilizers and the right side force surface. Low noise coaxial wire 
is used to trunk the signals to high gain amplifiers and low pass filters. 

The accelerometer characteristics arfe; given in- Table 24. 
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TABLE 24 


GROUND VIBRATION AND FLIGHT FLUTTER 
TEST ACCELEROMETERS 

Metrix Instrument Company 
Model 502 


Voltage Sensitivity 
Charge Sensitivity 
Capacitance 

Cross Axis Sensitivity 
First Resonance 
Linearity 


25 mv/g* 

150 PCMB/g* 
6000 pf* 

3% 

3000 Hz 
± 1 % 


*Nominal values. 


Angle of Attack and Sideslip Channels 

The vanes are Calspan Corporation designed and built units. Each 
vane consists of a symmetrical airfoil surface with' a mean aerodynamic chord 
located 2.0 inches aft of the vane axis for stability. The distance from the 
inboard edge of the vane airfoil to the aircraft skin is 7 inches. A balance 
weight is installed to eliminate error due to acceleration terms acting upon 
the vane mass. The synchro output of vane angle is AC amplified, scaled and 
demodulated. The 800 Hz ripple is removed in the third order filter. Correc- 
tion terms are added to the DC vane signals to compensate for location of the 
vane from the TIFS C.G.. These terms are ^-4/V and y^/^for the angle of attack 
vane and and y»-^/V for the angle of sideslip vane. 

= longitudinal distance from C.G. to angle of attack vane. 

= lateral distance from C.G. to angle of attack vane. 

= longitudinal distance from C.G. to sideslip vane. 

= vertical distance from C.G. to sideslip vane. 

Each vane signal is corrected for position error due to the air flow distribu- 
tion at their location on the fuselage. The vane deflections are larger than 
the actual aircraft angle of attack and sideslip changes. The angle of 
attack vane rotates 1,45 degrees for a one degree change in true angle of 
attack. The sideslip vane position calibration factor is 1.25 degre es . Th e 
corrected vane signals are notch filtered. The vane outputs are adjusted to 
obtain zero volts when the angle of attack vane is parallel to the fuselage 
reference line and the sideslip vane i_s parallel with the aircraft buttock line 
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Sensor Signal Conditioning 


The basic sensor signals contain conditioning networks which consist 
of AC to DC conversion, filtering and scaling. The format of the sensor chan- 
nels (sensor dynpiics and electronic signal conditioning) except for angle of 
attack and sideslip is as outlined in Figure 57. The angle of attack and side- 
slip channels are outlined in Figure 58. Sensors with AC outputs require a car- 
rier amplifier and phase sensitive demodulator. The AC to DC conversion is fol- 
lowed by a low pass filter to remove the 800 Hz ripple frequency. Notch filters 
or additional low pass filters are installed where, needed. Table 25 is a compi- 
lation of the dynamic transfer functions for each sensor channel except for angle 
of attack and sideslip. The transfer functions for the. angle of attack and 
sideslip channels are given in Table 26. In these tables, second order transfer 
functions have the form: ^ _ 


where v*# 
and £ 


£ o/£. 


/Af 


2(?)s 


u> * 
n 


aJ, 


+ 1 


r> 


- undamped- natural frequency in radians/second 
= damping ratio. 


First order transfer functions have the form: 

1 




//v - 


— - 1-1 
oo ^ ' 


where m - corner frequency (radians/sec). 


No attempt is made to indicate the gains of each step of the signal conditioning 
process. During calibration, the gain of each sensor channel is adjusted to 
produce the desired scale factor at the final output. 


ELECTRONIC SIGNAL CONDITIONING 



CONDITIONED 

SENSOR 

O UTPUT 


Figure 57 SENSOR CHANNEL FORMAT 


Scale factor information is added to the signal identification sym- 
bols. As an example, the symbol 0.5q identifies the rate of change of pitch 
angle signal with a scaling of 0.5 volts per degree per second of pitch rate. 
Since the TIFS electronics is based on a 10 volt signal level, the maximum us- 
able pitch rate is 20 degrees per second. For trigonometric functions, simi- 
lar scaling is used. The symbol 20 sin 9 identifies a pitch angle signal which 
can have a maximum value of 30° [20 (sin 30°) = 20 (0.5) = 10 volts ] . 
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| TABLE 25 

SENSOR CHANNEL CHARACTERISTICS 

The AC amplifier and phase sensitive demodulator response 
lag is negligible and is not included in the table. 


Aircraft 

Parameter 

Sensor 

Dynamics 

800 Hz 
Carrier 
Freq. 

Filter (1) 

Additional 

Filter 

Notch 

Filter 

1(2) 

Output 

Scale 

Factor 

Output 

Sensitivity 

Pitch Angle 

>25 Hz 

Yes 

No 

No 

20 sin 0 

10V = sin 30° 



Yes, 

. No 

No 

20 cos 0 

10V = cos 0° 






0.5 A9 (3) 

. 5V/deg 

Roll Angle 

>25 Hz 

Yes 

No 

No 

10 sin 0 

10V = sin 90° 



Yes 

i No 

No 

10 cos 0 

10V = cos 0° 






0.1 0.(4) 

,0.1 volt/deg 

Heading 

- >25 Hz 

Yes 

No • 

No 

10 sin i/f 

1 I0V = sin 90° 



! Yes 

' No 

No 

10 cos 

10V = cos 0° 

Pitch Rate 

1 

Yes 

1 

No 

0.5q 

1 .5V/deg/sec 

- ztVs , 

(nbf- !7& 

: d 7 

Roll Rate 

1 

Yes 

i 

No 

: | 0 . 2p 

! .2V/deg/sec 

sZ 

(30af 308 

iA ■ t1 

Yaw Rate 

i 

Yes 

i 

No 

, |0.5r 

: .5V/deg/sec 

! S z , 2f<7)s 
— IT, 1 

(i8z) z fez 

' *h + 1 
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ITABLE 25 (CONT'D) 


Aircraft 

Parameter 


Pitch Accel- 
eration 

Roll Accel- 
eration 

Yaw Accel- 
eration 

Longitudinal 

Acceleration 
§ C.G, 

Lateral 

Acceleration 

e c.g. 

Vertical 

Acceleration 
§ C.G. 


Lateral Accel- 
eration (sim- 
ulation pilot 
station) 

Vertical 
Acceleration 
(simulation 
pilot station) 


Sensor 

Dynamics 


800 Hz 
Carrier 
Freq. 

Filter (1) 


1 


. *<•?)* H \ 

{U5) z US 

1 

a 2 ^zc-7)s (1 

(t2Z)*' fZZ 
1 

S* ,er.7)s 

(107) z 107 ‘ ' 

1 

S* Z(.6)S fJ 

(Z51)*' r 251 

1 

S z 2(.8)s^, 

(238)* ' 23 6 
1 

S z 

1 -ff 

(238) z ^238 
f 


, z(- 7)s t . 
izzof zzo 


2(-7)s 
, + tf 


(zzoy 22o 


No 


No 


No 


No 


No 


No 


No 


No 


Additional 

Filter 


1 

$ z . Z(- 9)s 

-A / 

( 85 ) 2 85 

1 

s 2 2(.9)s 

~h 1 

(8S) Z ' <95- 

1 

s z zf-9)s 

+ 1 

(, 85 ) 2 1 85 

1 


S z 2(.9)s 

h i 

(6S) Z ‘ 85 

1 


s 2 . z(-V) s 

hL / 

' 85 


1 



s‘ z(-8)s 


+ 


{esr 85 


-hf 


(esr 


, „ 

-h -hi 

85 


( asy 


2(-8)s 

85 


hi 


>tch 

Liter 

Output 

Scale 

Factor 

r— 

! 

Output 

Sensitivity 

Yes 

O.lq 

. lV/deg/sec 2 

Yes 

0.2p 

.2V/deg/sec 2 

Yes 

0. lr 

. lV/deg/sec 2 

Yes 

20 n 

X 

2V/0 . lg 

Yes 

10 n 

y 

lV/O.lg 

Yes 

4 

4V/g 


4 An^ 
I 1 (5) 

4V/g 

No 

10 n v 
y P 

IV/. lg 

No 

4 n. 

4V/g 

1 

44no 

h 

1(6) 

4V/g 



TABLE 25 (CONT’D) 




800 Hz 





Aircraft 

Sensor 

Carrier 

Additional 

Notch 

Output 

Output 

Parameter 

Dynamics 

Freq. 

Filter 

Filter 

Scale 

Sensitivity 



Filter Cl) 



Factor 


Radar 

1 

No 

j No 

No 

0.01 h 

O.OlV/ft 

Altimeter 

£-+ 1 


20 s 

(***X£") 





.10 


No 

0.2 h r 
>(7) 

.2V/ft/sec 

Right Wing 

>1000 Hz 

No 

fO 

1S4- 1 

No 


2.56| V/g 

- Tip Acceler- 

Bandpass 





ometer 

>1000 Hz 
Bandpass 

No 

1 

No 


1.73 jV/g 

Left Wing 
Tip Acceler- • 
ometer 

z(-lU + , 
(6 ?.a) z <,z.s 






Right Horiz. 

>1000 Hz 

No 

lO 

ls+ 1 

No 


2.45 IV/g 

Stab. Accel. 

, Bandpass 


7 



2.60 jV/g 

Left Horiz. 

>1000 Hz 

No 

s 2 , *f'7 > S y . , 

No 


Stab. Accel. 

Bandpass 


1Z(o 




Vertical 

>1000 Hz 

No 

70 

No 


2.78 jV/g 

Stab. Tip 

Bandpass 


fs+ 1 


. 


Right Side 

>1000 Hz 

No 

ro 

rs+-i 

No 


2.23 | V/g 

Force Sur- 
face Top 

Bandpass 






Right Side 

>1000 Hz 

No 

70' 
ts+ 1 

No 


2.52 jV/g 

Force Sur- 
face Bottom 

Bandpass 
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TABLE 25 (CONT'D) 


NOTE: 


(1) 800 Hz carrier frequency filter has transfer function g 


IN 


(— — + 1 V 

\>000 )\ 


Z _ 2 (. 8 )$ 

(3 &0) zim 380 


•b 1 


j (2) jThe notch filter has a transfer function of' — — 

in 


o _ 


( 1 08 ) 2 


+ 1 


( i 08 Y 


?(.7)$ 

f 08 . 


■b i 


i (3) O.5A0 is a small angle linear approximation of the conditioned 20 sin 9 signal, 

identifies this signal as change of pitch angle after system engagement. 


The V 


j(4) 0.1 0 is the output of a nonlinear function- generator. The function generator input is the 

conditioned 10 sin 0 signal. 

| (5) 4 A n^ is the same signal as 4 n^ . with the Ig gravity term removed. 


< j(6) 4 A n^ is the same signal as 4 n^ p with the lg gravity term removed. 
(7) The 0.2 h^ is an electronically differentiated 0.01 h^ signal. 






p,q = filtered rate gyro signals 
~V = — t rue" Airspeed’ 

, (Eq ) / a 

-4 = 31.8 ft. K&x + Q: q} = * 

* 3.8 ft. 7 


Position 

Error 

Correction 


+ (Equation 1) 



Position. rl. 

x Error . ' h ~~v~ + ~~v~ (Equation 2) 

Correction v 

- 31.8 ft. 

Figure 58 ANGLE OF ATTACK AND SIDESLIP VANE CHANNELS 


P » r z _f i 1 tered rate gyro _si g nals 
a true ai r speed + fi ) » 

3 - 3.0 ft. 
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|. TABLE 26 

Cf AND & VANE CHANNEL CHARACTERISTICS 


Aircraft 

Parameter 

Vane 

Characteristics 

j 

or 

Additional 

Filters 

800 Hz 
Carrier 
Filter 
CD 

Position 

Error 

Correction 

m 

Notch 

Filter 

C3) 

Sensor 

Signal 

Scale 

Factor 

Angle of 
Attack 


Yes 

0.690 

+ 1**2 
V V 

Yes 

•Step*?) 

. 5V/ deg 

Angle of 
Attack Rate 

. Zs 

s z as L v 

(C,7) 2 07 





. 1 Cq^. > & 9 ) 

. lV/deg/sec 

Angle of 
Sideslip 


Yes 

0.800 

_ r \ 
Y V 

Yes 

. 5 (P x +/8y) 

,5V/ deg 

Angle of 
Sideslip 
Rate 

. Zs 

32 +t 

(on) 07 





j) 

.IV/deg/sec 


CD 

GO 


800 Hz carrier Frequency Filter = 


[lira + ')( 


(380) 

q and r terms are filtered gyro signals from ‘Table 9. 

= 31.8 ft -B s = 3.0 ft 

= 3.8 ft Z+ = 31.8 ft 


2 + 


- true airspeed from | Equation 13. 2 

( tow 


'h 1 


(3) Notch filter transfer function is 




( tosY 


to& 


2(.8)s 
, 380 


+ 1 









amplitude 

Figure 59 
Figure 60 
Figure 61 
Figure 62 
Figure 63 
Figure 64 

Figure 65 

Figure 66 
Figure 67 
Figure 68 
Figure 69 
Figure 70 
Figure 71 
Figure 72 
Figure 73 


The following plots are presented for the function generator and 
and phase responses of the .sensor f i-l-ters . 

Sin 0 to 0 function generator. 

800 Hz carrier freq, filter amplitude vs. frequency. 

800 Hz carrier freq. filter phase vs. frequency. 

Notch filter amplitude vs. frequency. 

Notch filter phase vs. frequency. 

Second order low pass vilter, = 85 rad/sec amplitude vs. 
frequency. 

Second order low pass filter, co n - 85 rad/sec phase vs_. 

frequency . 

Radar altitude differentiator amplitude vs. frequency. 

Radar altimeter differentiator phase vs. frequency. 

Low pass filter = 62.8 rad/sec amplitude vs. frequency. 

Low pass filter co n = 62.8 rad/sec phase vs. frequency. 

Low pass filter u> n = 126 rad/sec amplitude vs. frequency. 

Low pass filter oO n = 126 rad/sec phase vs. frequency. 

Vane differentiator amplitude vs. frequency. 

Vane differentiator phase vs. frequency. 
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-0./1 0 , volts 










AMPLITUDE, E./E. 



Figure 60 AMPLITUDE VS FREQUENCY, 800 Hz CARRIER FREQUENCY FILTER 
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AMPLITUDE RATIO, E/E, 
















PHASE LAG, E/E. 



Figure 63 PHASE VS FREQUENCY, NOTCH FILTER 



AMPLITUDE RATIO’, E /E. 


















PHASE ANGLE 


10 

!- 20 , 
i 
i 
i 
i 
i 

|-40 

t 

I 

]-60 
1-80 
-100 
;-i2o 

| -HO 
1-160 













AMPLITUDE RATIO, E /E. 



Figure 66 AMPLITUDE VS FREQUENCY, |RAD> 
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PHASE ANGLE 
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jAMPLITUDE RATIO, E /E. 
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PHASE ANGLE 














! AMPLITUDE RATIO E /E, 
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PHASE ANGLE 










AMPLITUDE RATIO, E /E. 



! FREQUENCY, Hz 

|| Figure 72 AMPLITUDE VS FREQUENCY, VANE DIFFERENTIATOR 



PHASE ANGLE 

















3.4.2 Group Two. Air Data Signals 

The Air Data System, shown in Figure 74 provides output signals 
of altitude, altitude rate, indicated airspeed, true airspeed, true airspeed 
rate and Mach number. A block diagram of the signal flow through the Air 
Data Computer is illustrated in Figure 75- Diode function generators, mech- 
anization of inertial equations and complementary filters used to generate 
these signals are compiled in Table 27. Additional sensors in combination 
with the group one signals are required. These air data sensors are: 


p s 

- 

static pressure 




9c 

- 

dynamic pres. 

— 


T c 

- 

air temp. 

► 

INERTIALLY 





COMPENSATED 

n z 

- 

longitudinal accel. — 


AIR DATA 





COMPUTER 

** 

- 

lateral accel. 

— »** 


"i 

- 

normal accel. 



e 

- 

pitch angle 




- 

roll angle 

1*- 


<x T 

- 

inertial angle of 






attack 



fix 

- 

angle of sideslip - — 

— 


i> 

- 

roll rate 

— - 


9 . 


pitch rate 

— : 


r 

- 

yaw rate 

— 

' 


Altitude h 

* 

Altitude Rate In 

True Airspeed V 

• 

Airspeed Rate V 

Indicated Airspeed 
Mach No. M 


Figure 74 AIR DATA SYSTEM INPUT-OUTPUT SIGNALS 
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j PRESSURE 

I ALTITUDE h. 




SIN / 
FLIGHT 
PATH 
:angle 


j Figure 75 INERTIALLY COMPENSATED AIR-DATA BLOCK DIAGRAM 














TABLE 27 


AIR DATA EQUATIONS 


AIR DATA FUNCTION GENERATORS 

hp = 4 (f» § ) Figure 76, Table 30_ (Eq. 3) 

M - fg c/'Ps ) Figure 77, Table 31 _ (Eq. 4) 

{*# - fj C^c ) Figure 78, Table 32 (Eq. 5} 

y'i - *4 (% c ) Figur e 79, Table 33 (Eq. 6) 

AIR DATA EQUATION 

Vp “ <p£*77 M-ft'x (Eq. 7) 

M from Eq. 4, ~jt/^ from Eq. S 
INERTIAL EQUATIONS 

h, x = q sin 6 - sin 4> - n^ sos (Eq. 8) 

All variables are filtered sensor 
signals from GROUP ONE. 

v z “ 3 [Px 7- n y fa +- n^ot x - sin y] (Eq. 9) 


Az 

sin y 


is obtained from Eq.. 17 
is obtained from Eq. 18, 
is obtained from Eq.. 14 


Acceleration. terms are 
obtained from GROU P ONE 


h = h 


COMPLEMENTARY FILTER EQUATIONS, 


s. 


T f r 


(. Zj $ + 1'Y 


(•zqs-*- f)‘ 


Zf = Zsec. (Eq. 10) 


/? from DFG. Eq. 3 b: T from Eq. 8 
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TAB LE 27 (cONT ' D) 


h * h 


Vp 5 +■ 1 


1 ; Zz 

+ h 


71? S +■ 1 


T z - 1 sec. (Eq. .11) 


K= /< 


hp from DFG Eq. 3 h from Eq. 10 

* ') 


r (r 3 s ±l) z + Vx (r 3 s -m)* 

Vp from Eq. 7 V x from Eq. 9 


2 ^ S ( 2 


r z = 


10 sec (Eq. 12) 


V = V. 


1° TT+s-h i 


•h V 


'-4 


r, s + 1 


= IO sec. (Eq. 13) 


Vpj from Eq. 7 K from Eq. 12 

AUXILIARY EQUATION . • 

h 

sc/ 7 7 = 7 

h ‘ from Eq. 10 

^ from Eq. 13 


(Eq. 14) 
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Static Air Pressure 


This absolute pressure transducer is .pneumatical-ly- connected to 
the basic aircraft copilot static source. The sensor output is DC and pro- 
perly scaled for direct use in the system. The characteristics of this sensor 
are listed in Table 28 . 

TABLE 28 

STATIC AIR PRESSURE 

Rosemount Engineering Company 
Model 830BC 


31" Hg = 10V DC 
10" Hg = 0V DC 

Linearity ± .07% of full scale 
Repeatability < ± .02% of full scale 
Hysteresis + .015% of full scale 
Resolution < ± .02% of full scale 

For step input change, T < .015 to reach 63% of final value. 
Dynamic Air Pressure 

This differential pressure transducer is connected to the copilot 
pitot and static pressure sources. The DC output of this sensor is used 
directly in the system. The applicable characteristics of this -sensor are 
listed in Table 29. 


TABLE 29 

DYNAMIC AIR PRESSURE 

Rosemount Engineering Company 
Model 831BC 

q 0" Hg = 0V DC 

5" Hg = 10V DC 
Linearity t .07% of full scale 
Repeatability < ± .02% of full scale 
Hysteresis ± .015% of full scale 
Resolution < i .02% of full scale 

For step input change, Z < .015 to reach 63% of final value. 
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Outside Air Temperature 


A standard aircraft type EDISON temperature bulb is used for this 
measurement. A temperature change causes a proportional change in bulb re- 
sistance, The output of the temperature bulb circuit is scaled to produce 
the required temperature signal. 

Air Data Function Generators 


Diode function generators are used where nonlinear transfer func- 
tions are required. In the air data computer, function generators are used to 
obtain: 


a) Pressure altitude from static pressure, f 1 , Figure 76 , Table 
30- 


b) Mach number from a ratio of dynamic and static pressures, f_. 

Figure 77 , Table 31 . 2 

c) Square root of temperature (Kelvin) from temperature (Centi- 
grade) , fg. Figure 78 , Table 32. 

d) Indicated airspeed from dynamic pressure, f , Figure 79 , 

Table 33 . 4 

Actual function generator measurements are presented in the tables. 
A more complete description of the function generators can be obtained from 
TIFS Memo No. 53 of 21 March 1967 entitled "Diode Function Generators". 

Air Data Computations 

The air data signals required for the TIFS operation must have 
high bandwidth properties, high resolution and accuracy, yet be insensitive 
to air turbulence. To accomplish this, a complementary filter technique is 
used to blend air derived signals that are accurate at low frequencies with 
inertially computed signals which are free of turbulence at higher frequencies. 
For example, altitude derived from static pressure measurement is accurate at 
very low frequencies and inaccurate at high frequencies due to air turbulence 
contamination. Altitude derived from the double integration of vertical accel- 
eration is accurate at high frequencies, however,' d.c. offsets in the vertical 
acceleration measurement result in large errors in the computed altitude sig- 
nal due to the double integration process. The resultant composite signals 
are accurate in amplitude and phase throughout a wide bandwidth. 

The inertially computed signals h z and Vj. are formed from the 
solutions of equations 8 and 9 in Table .2.7 - The blended signals h , h , V 
and V are computed from equations 16 through 13. A procedure for determining 
the breakpoint frequency for the complementary filters is detailed in Appendix 
A. Flight path angle (sin y ) is computed in equation 14. 
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Figure 76 STATIC PRESSURE TO PRESSURE ALTITUDE FUNCTION GENERATOR 
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[TABLE 30 

TO ALTITUDE FUNCTION GENERATOR 


■fs 

In, Hg 

•Ps/* 

V. DC 

INPUT 

+ A - 20 

V.DC 

DESIRED 
OUTPUT 
+ hp/2500 
V.DC 

DESIRED 
+ hp/2500 
Feet 

ACTUAL 
OUTPUT 
+ hp/2500' 
V.DC 

31.354 

-10.451 

11.354 

-0.5200 

-1300.00 

-0.5160 

. 30.00 

-10.000 

10.00 

-0.00289 

-72.36 

-0.0020 

29.00 

-9.667 

9.00 

+0.3451 

862.84 

+0.3461 

28.00 

-9.333 

8.00 

+0.7300 

1824.90 ' 

+0.7293 

27.00 

-9.000 

7.00 

+1.1261 

2815.29 

1.1255 

26.00 

-8.667 

6.00 

+1.5343 

3835.75 

1.5329 

25.00 

-8.333 

5.00 

1.9554 

4888.46 

1.9537 

24.00 

-8.000 

4.00 

2.3904 

5975.91 

2.3888 

23.00 

-7.667 

3.00 

2.8403 

7100.69 

2.8381 

22.00 

-7.333 

2.00 

3.3063 

8265.84 

3.3030 

s 21.00 

-7.000 

1.00 

3.7909 

9477.16 

3.7877 

20.00 

-6.667 

0.00 

4.2924 

10730.9 

4.2904 

19.00 

-6.333 

-1.00 

4.8157 

12039.3 

4.8121 

18.00 

-6.000 

-2.00 

5.3618 

13404.5 

5.3601 

17.00 

-5.667 

-3.00 

‘ 5.9331 

14832.7 

5.9299 

16.00 

-5.333 

-4.00 

6.5322 

16330.6 

6.5306 

15.00 

-5.000 

-5.00 

7.1625 

17906.2 

7.1624 

14.00 

-4.667 

-6.00 

7.8278 

19569.4 

7.8277 

13.00 

-4.333 

-7.00 

8.5327 

21331.8 

8.5353 

12.00 

-4.000 

-8.00 

9.2830 

23207.6 

9.2858 

11 . 00 - 

-3.667 

-9.00 

10.0858 

25214.6 

10.089 

10.00 

-3.333 

-10.00 

10.9500 

27375.1 

10.952 




OUTPUT, volts 
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TABLE 31 


%/f s TO MACH NO. FUNCTION GENERATOR 




DESIRED 

ACTUAL 


rc 


OUTPUT 

OUTPUT 

M 

Ps 

V. DC 

V. DC 

MACH NO. 

0.00 

+10.00 

0.914 

0.910 

0.0457 

0.01 

+9.00 

2.388 

2.382 

0.1194 

0.02 

+8.00 

3.368 

3.363 

0.1684 

0.03 

+7.00 

4.118 

4.114 

0.2059 

0.04 

+6.00 

4.748 

4.752 

0.2374 

0.05 

+5.00 

5.298 

5.302 

0.2649 

0.06 

+4.00 

■ 5.794 

5.793 

0.2897 

0.07 

+3.00 

6.248 

6.246 

0.3124 

0.08 

+2.00 

6.668 

6.667 

0.3334 

0.09 

+ 1.00 

7.060 

7.060 

0.3530 

0.10 

.0:00 

' 7.430 

7.430 

0.3715 

0.11 

-1.00 

7.780 

7.775 

0.3890 

0.12 

-2.00 

8.112 

8.110 

0.4056 

0.13 

-3.00 

8.430 

8.426 

0.4215 

0.14 

-4.00 

8.734 

8.734 

0.4367 

0.15 

-5.00 

9.026. 

9.024 

0.4513 

0.16 

-6.00 

9.308, 

9.308 

0.4654 

0.17 

-7.00 

• 9.580 

9.577 

0.4790 

0.18 

-8.00 

• 9.842 

9.842 

0.4921 

0.19 

-9.00 

10.096 

10.091 • 

0.5048 

0.20 

-10.00 

10.342 

10.341 

0.5171 







volts DC 
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j'TABLE 32 

± c TO -ftK FUNCTION GENERATOR 





INDICATED AIRSPEED, V.j/30 
— volts 
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j TABLE 33 

TO INDICATED AIRSPEED FUNCTION GENERATOR 


9c 

In. HG 


lb/ft' 


-2 ? c 
V. DC 


INPUT 
4 + 10 

V. DC 


DESIRED 
OUTPUT 
Vi /30 
V. DC 


ACTUAL 
OUTPUT 
Vi/ 30 
V. DC 


V i 

Knots 


0.00 

0.25 

0.50 

0.75 

1.00 

1.25 

1.50 

1.75 

2.00 

2.25 

2.50 

2.75 

3.00 

3.25 

3.50 

3.75 

4.00 

4.25 

4.50 

4.75 

5.00 

5.25 


0.000 

17.6825 

35.365 

53.0475 

70.73 

88.4125 

106.095 

123.7775 

141.46 

159.1425 

176.825 

194.5075 

212.190 

229.8725 

247.555 

265.2375 

282.92 

300.6025 

318.2850 

335.9675 

353.65 

371.3325 


0.000 

0.500 

1.000 

1.500 

2.000 

2.500 

3.000 

3.500 

4.000 

4.500 

5.000 

5.500 

6.000 

6.500 

7.000 

7.500 

8.000 

8.500 
9.000 

9.500 
10.000 
10.500 


- 10.00 

-9.00 

- 8.00 

-7.00 

- 6.00 

-5.00 

-4.00 

-3.00 

- 2.00 

- 1.00 

0.00 

1.00 

2.00 

3.00 
4.0Q 
5.00. 

6.00 

7.00 

8.00 

9.00 

10.00 


1.000 

2.4035 

3.3936 

4.1501 

4.7852 

5.3423 

5.8443 

6.3042 

6.7300 

7.1276 

7.5025 

7.8570 

8.1959 

8.5188 

8.8259 

9.1244 

9.4114 

9.6879 

9.9549 

10.2140 

10.4646 

10.7102 


1.0000 

2.4033 

3.3929 

4.1517 

4.7834 

5.3423 

5.8469 

6.3036 

6.7300 

7.1282 

7.5023 

7.8562 

8.1959 

8.5132 

8.8253 

9.1206 

9.4106 

9.6847 

9.9539 

10.2098 

10.4658 


30.00 

72.10 

101.81 

124.50 

143.55 
160.27 

175.33 
189.13 
201.90 
213.83 
225.07 
235.71 
245.88 

255.56 
264.78 
273.73 

282.34 

290.64 

298.65 
306.42 
313.94 
321.30 




3.4.3 Group Three. Angle of Attack, Sideslip and Gust Computations 

High feedback gains of angle of attack, sideslip and their rates 
are desirable for the TIFS control system, To achieve these feedback gains, 
inertial signals must be used. The vane signals contain gust information 
which make them unusable for feedback. 


Inertial signals are computed for angle of attack oc x , angle of 
sideslip fi x , angle of attack rate cc x and angle of sideslip rate fix • 
Gust signals are computed from the difference between the vane signals and 
these inertial computations. The equations for these computations are listed 
in Table 34. 


The inertial computations are the following: 


® Inertial angle of attack rate is computed from the following equa- 
tion mechanized in the TIFS sensor electronics. 


generated from . 57.30 r ,i P ft t 

% force equation *** ~ — 'p [ ^ + cos $ J + ^ ~ 5 j $ 


(Eq. 15) 


where: , p, q, cos 0 are conditioned basic sensor signals 

* from Table 9. 

V is obtained from equation 13. . ■ 

fi x is obtained from equation 18. 


• Inertial angle of sideslip rate is computed from: 


generated from _ 57.3 g 

y force e quation Px y [ 5 <j 


+ sir) 


fl- 


■P &JL 
57. 5 


where: , p, r, sin 0’ are conditioned basic sensor signals^ 

V is obtained from equation 13. 
ccj is obtained from equation 17. 


(Eq. 16) 


The inertial angle of attack is computed from the following equa- 


tion. 


cx-x - -57. 3 


j y + sin & 

cos (p 


Ax 

57.3 


S'z.n 


in (j) 


(Eq. 17) 


where: sin y is obtained from equation 14. 

sin 0, sin 0, cos 0 are co nditi oned basic sensor signals 
from Tabl e 9. 

fijr is obtained from equation 18. 
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Inertial angle of sideslip cannot be obtained in the same manner 
as inertial angle of attack. In TIFS, a complementary filter Is used to gen- 
erate inertial sideslip angle. In this filter, low frequency information from 
the corrected angle of sideslip vane is blended with the inertial angle of 
sideslip rate term. The complementary filter equation mechanized is: 

VAME /MirRT/AL ,f s . \ 

' A — -> / ' — 

- (^^ocTTf r + /** JzTTly*- ’ rwsec <&• 18 ^ 


fix 


( j3 r + fig) is 'the angle of sideslip vane signal corrected for position error 
and for location from the CG (Table 26). 

fij- is obtained from Equation 16. 

Angle of attack and sideslip gust terms are computed from the difference be- 
tween corrected vane and inertial ly computed terms. The gust equations mech- 
anized in the TIFS sensor system are: 


tegusi = (ot x + &g) (Eq. 19) 

where: ( cc z + oc.^) is corrected angle of attack vane from Table 26. 

ccx is computed from Equation 17. 

A small DC term may be present on the <X 3 ust signal because of 
improper vane alignment in the local airstream at all airspeeds. A DC washout 

(high pa ss) fi lter with higher frequency rolloff chara cter istics is available 

on the computer patch panel to eliminate this DC term. The filter transfer 
function is - ' $ ’ s A(fd + 1 )(A + i )) ' 


= (/5l ) - A CEq. 20) 

where: {fix* fig) is corrected angle of attack vane from Table 26,. 
fi L is computed from Equation 18. 


c& 


gust 


= ( ce x + a 3 ) - a. z 


(Eq. 21) 


where: ( o& x is the derivative of the corrected ang-le of attack vane 

signal from Table 26.. 


cc z is computed from Equation 15. 


A 


gust 


“ C fix + fig 5 “ fix 


(Eq. 22) 


where: (/S z + /S„) is the derivative of the corrected sideslip vane signal 

y from Table 26 . 


and 


fix 


is computed from Equation 16. 
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TABLE 34 


CX, X : 

Ax ' 
Ax = 


ANGLE OF ATTACK AND SIDESLIP EQUATIONS 


57.3 


ANGLE OF ATTACK COMPUTATIONS 
— sin -h Sin & ~ Ax / 3 sin <fi 


COS <f> 


-J 


(S7.3 )(32.2 ) 
V 


n„ + c.os <f> 


+ ? " 


Mr 

57-3 


CEq. 17) 
(Eq. 15) 


ANGLE OF SIDESLIP COMPUTATIONS 


( 57 . 3 )( 32 . 2 ) 


V 

VANS 

Ax * Ag 
(Zs + i) ' 



(TTS + 1) 


X = f Sec. 


CEq. 16) 
(Eq. IS) 


GUST EQUATIONS 
Corrected 



Vane 

i Signals 

Inertial 



cc - 

gust 

(o:j 

"1 

* *3) 

CC X 

CEq. 

19) 

A 

gust 

(Ai 

+■ fig ) 

Ax 

CEq. 

20) 

as = 

gust 


+ *e) 

- os x 

CEq. 

21) 

A gust 

(/ 

+ Aq) 

1 

H 

CEq. 

22) 
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Excellent correlation of vane and inertially computed angle of 
attack and rate terms for an elevator doublet can be seen in the flight record. 
Figure 80. The inertial angle of attack and inertial angle of attack rate 
signals at the TIFS C.Q'. are derived from independent computations (Equations 
17 and 15) . These signals can be directly compared with the vane angle of 

a ttack s ign al and its rate term corrected for position err or and location f rom 

the C.G. The an gle of sidesli p v ane rate signal_and_ the inertial ly computed 
sideslip rate signal for a rudder doublet can be directly compared in Figure 
81- . The sideslip vane is corrected for position error and for location from 

the C.G. This correc ted vane signal can be directly compa red with .the ou tput 

of the sideslip complementary filter '{fix) • Mechanization of the d x fix 
computations were included to provide noise free signals that can be used in 
high gain feedback loops. A comparison of the noise content between the dif- 
ferentiated vane signals, and the computed signals can be made from Figures 80 
and 81. In addition, the inertial computations allow the- generation of gust 
terms (Figures 86 and 87) which are the difference between the corrected vane 
signals and the inertial computations. 

If the vane and inertially computed angle of attack signals are 
compared for lateral aircraft motions, a difference between the vane and 
inertial signals is apparent (Figure82). The independent computations of 
OC x (Equation 17) and 'a x (Equation 15) maintain correlation during these 
lateral maneuvers . Therefore, it is suspected that the angle of attack vane 
is affected by aircraft sideslip. To check this effect, a- percentage of the 
corrected angle of sideslip .vane signal was subtracted' from the corrected 
angle of attack vane signal. The equation used was the following: 

^“compensated = +O.S (ac x +,<x g ) - K On? ( + fy) . 

a dramatic improvement in correlation between the compensated angle of attack 
vane and the inertially computed signals was obtained for an aileron doublet 
(Figure 83) and a rudder doublet (Figure 84) with a compensation factor of 
K = 0.125. 

To check the correlation of angle of attack vane deflection with 
sideslip, records should be obtained at various airspeeds and angles of attack 
for lateral-directional inputs. 

The angle of sideslip vane response for an elevator doublet input 
is presented in Figure 85. No difference between angle of sideslip vane and 
computed sideslip is apparent from this record. The above record was taken 
at a zero sideslip angle. To determine if crosscoupling exists from angle 
of attack to sideslip at steady state sideslip angles, additional flight 
checks should be made. 

The aircraft configuration for the records in Figures 80 to 85 was 
165 knots indicated airspeed, 10,000 feet MSL, gear up, flaps at zero, and the 
angle of attack and sideslip vanes mounted on the nose faring,. For the records 
of light to moderate turbulence, the airspeed was 150 knots indicated, 1800 
feet MSL, gear up, and flaps at zero.. The oc. and a. gust records are presented 
in Figure , 86 and the fi and /3 gust signals are presented in Figure 87'. From 
these -figures,' ht;ls evident that the inertially computed signals remain clean 
in the presence o f l ight to moderate turbulence. 
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Figure 80 ANGLE OF ATTACK RESPONSE TO ELEVATOR DOUBLET AT 165 kt IAS 
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Figure 81 AIRCRAFT SIDESLIP RESPONSE TO A RUDDER DOUBLET AT 165 kt IAS 
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Figure 82 ANGLE OF ATTACK RESPONSE TO A RUDDER DOUBLET AT 165 kt IAS 
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Figure 83 ANGLE OF ATTACK RESPONSE WITH SIDESLIP CORRECTION 
' — FOR AN AILERON INPUT 
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Figure 84 ANGLE OF ATTACK RESPONSE WITH SIDESLIP CORRECTION 
FOR A RUDDER INPUT ■ 
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Figure 85 ANGLE OF SIDESLIP RESPONSE TO AN ELEVATOR INPUT AT 165 kt IAS 
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Figure 86 ANGLE OF ATTACK RESPONSE IN LIGHT TO MODERATE TURBULENCE AT 150 kt IAS 
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Figure 87 ANGLE OF SIDESLIP RESPONSE IN LIGHT TO MODERATE 
- TURBULENCE. AT 150 kt IAS 
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Complementary Filter Fidelity 


To verify the overall fidelity of the TIFS . complementary filters, 
flight, records of the input and output signals of the V, V, h and fi £ com- 
plementary filters were evaluated. The derivative input terms V z , *fi x and 
fix were integrated electronically during ground playback to aJllow direct com- 
parison with the filter output signals. . If close correlation is maintained 
between signals and their derivatives oh both a short term and a long term 
basis, proper operation of the complementary filter is verified. 

During the flight record ground playback, the inertial ly computed 
V x input signal was integrated electronically for approximately 20 seconds. 

The V complementary filter output can be verified by direct comparison to the 
V> input in Figure 88. 

For the V complementary filter, a similar. comparison of the V and 
V p input signals can be made by int'egration of the V input signal. Proper 
V complementary filter operation. is indicated by comparison of the V output 
signal with both the integrated V and the V input signals in Figure 88. 

In the h complementary filter, accurate correlation of the iner- 
tial input signal with the filter output can be seen from Figure 89. The hj 
input signal was electronically integrated during playback of .the flight 
record and can be directly compared with the h output signal. The sin y com- 
puted signal waveform can be directly compared with the h signal in Figure 89. 

A magnitude comparison can also be made, with proper scaling, since the record 
was taken at a relatively constant indicated airspeed of 155 knots. 

• 

The input 'signals f$ x and {fi x + fi g ) to the fi x complementary 
filter can be compared by integration of the fix signal during ground play- 
back. The integrated fi x signal and the ( fit + fig') signal have excellent cor- 
relation (Figure 90) with the fi x complementary filter output.. Also shown 
in Figure 9Q is the differentiated sideslip vane signal {fix + fig) which can 
be directly compared with the inertially computed fi z signal. 
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DIGITAL TAPE RECORDING SYSTEM 


. TIie primary flight data recording device on TIFS is a system 

built around an Ampex 9-track, 800-byte-per-inch digital recorder. Any pre- 
selected 58 analog data signals available in the TIFS electronics can be read- 
j.y. Patched to the recording system. Each signal is sampled, converted to a 
digital word, bookkeeping digits added, formed into an IBM 370/195 compatible 
format, parity bits added for error detection and bit location specified. This 
is accomplished m an. interface/control unit built by Interland Electronics to 
Gal span Corporation specifications. The input data channels can be sampled at 
selectable rates of 50 or 100 updates of each channel per second. The 10V max- 
imum analog signals are converted to a 10 bit-plus-sign digital word which re- 
sults m discrete bit steps of 10 MV. Each analog input channel is filtered 
with a second order low pass filter prior to sampling to avoid aliasing errors. 
The cutoff frequency is automatically selected to 7.5 Hz and damping of .8 for 
50 sample-per-second rate and 18 Hz with damping of .8 for the 100 sample-per- 
second rate. - F F 


A complete ground playback 
poses. Any 8 of the 58 channels can be 
strip chart recorder in analog form. 

Accuracy: 


system is available for quick-look pur- 
selected and played back on an 8-channel 


Total system accuracy from input signal level to digital codes 
is ± 0.2% of full scale for steady state signals. This includes multiplexing, 
digitizing linearity, comparator and reference drift, and quantizing errors. 

No aperture or aliasing errors are included in this figure since these values 
depend upon signal rate of change and frequency content. 


Recorder: 


Input Level: 

0 VDC ± IV for "0", 4 VDC ± IV for "l". 

Input Impedance: 

10K ohm minimum. 

Input Format: 

NRZ (nonreturn to zero) mark (NRZ-IBM) . 

Bit Packing: 

800 bytes per inch. 

Interchannel Timing: 

6 microseconds max at 60 in/sec. 

Dropout Rate: 

Approximately one bit in 10 6 bits. 

Temperature: 

0°C to 60°C. 

Vibration: 

5 g random per procedure XIII of MIL-E-5272C. 

Shock: 

25 g for 11 milliseconds . 

Altitude : 

12,000 feet. 

Tape Speed: 

15 in/sec (100 samples/sec/channel) 
7.5 in/sec (50 samples/sec/ channel) 
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Performance Specification 


... . The di g ital t a P e recording system operates to the following 

specifications over the entire system environment specified. 

Input : 


Signal Capacity: 

58 channels of analog information and up to 
6 true-false inputs (0 and +28 VDC) . 

Frame Size: 

60 channels normal, or 30 channels when one 
multiplexer is not installed. 

Frame Rate: 

100 frames per second and 50 frames per 
second for 60 channels. 

Clock Stability: 

System timing shall derive from a crystal 
oscillator stable within ± .05%. 

Impedance : 

1 megohm minimum shunted by less than 100 
picofarads . 

Level: 

0 to 5 volts bipolar. 

Encoding: 

• 

Code: 

Each sample is encoded into a 10 bit-plus-sign 
code of sign-plus.- absolute-value configuration. 

Linearity: 

Conversion linearity will be within t .05% 
- 1/2 LSB (least significant bit). 

Aperture: 

Encoding time is approximately 12 micro- 
seconds permitting data frequencies up to 
10 Hz without increase in error. 
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APPENDIX A 

INPUTS TO. ISIC PROGRAM OF FLEXSTAB 


INTRODUCTION 


This appendix contains a reproduction of the inputs to the Internal- 
Structural Influence Coefficient Program which generated the parameters used 
in obtaining the normal mode results presented in this report. The structural 
model made provision for representing vertical and lateral translation of the 
engine on its mounts and rocking, bending, and torsional degrees of freedom of 
the side force surfaces. These motions were essentially eliminated in the case 
considered in this appendix by assuming very high stiffnesses for the beams 
used to represent these motions. 


- .V;’-- Figure A-l is a reproduction of the FLEXSTAB elastic axis plot and 

' .j, V- . - * 

shows. 'the 'locations of the concentrated masses used to model the inertia prop- 

V' . , , >5 .*■ • 

erties of TI-FS . - 


The FLEXSTAB thin body representation was used for the wing, hori-^ 
zontal tail, vertical tail, and side force surfaces while the FLEXSTAB slender 
body representation was used for the fuselage and nacelle. It had originally 
been intended to use an additional slender body at Y = 9.182 m because the side 
force, surface (SFS) elastic axis intersects the wing forward of the wing • elastic 
axis.-; A strategem introduced to avoid this involved assuming that the SFS, 
elastic axis intersected the wing elastic axis and then was swept abruptly for- 
ward to its true position. The swept portion of the SFS elastic axis was made 
very stiff so that it would have no effect on the computed vibration modes. 


FLEXSTAB does not allow concentrated masses to be located at. junc- 
tion points. Consequently, the masses which were to have been located in a 
slender body at Y = 9.187 m were divided and placed slightly inboard* and out- 
board of this station. 
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A few details should be mentioned which cannot be seen on Figure 
A-l . FLEXSTAB requires the wing elastic axis to be perpendicular to the X axis 
at the junctions with the -fuselage-, -nacelle and SFS- elastic axis . This was 
accomplished with short jogs in the wing elastic axis. 

A discussion of the cards used in the I SIC input is given in 

Reference 7, 
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33rX ™ 

308 

0.124 2.643 

0.213 

3.361 

0.302 

4.078 

* 33F2 ' v - 

— : — : — -3 39"'‘ — 

— 07391 47797' 

0 .480 

5 #513 

7. S69 - 

6« 232 

33F3^ T 

f. • .. . 910 ; 

0.603 6.503 

0.613 

6.589 


- 

33F4 


“'SSTIFFNESS'TJTrnr 





, . 7J-5F ' 

1 ’31 2 , 

0.0 3.65 

E082.1C 

EOS 



- 35 Fl 

■"’* • • -.'..r-313~: J T 

0.0“ — “ : 3769 

E 062.10' 

“E“tT6 



0 •„ r 2 $3 

r : - : v" 3i4- .* 

;0'.0 " 3.28 

E 061.85 

E 06 
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/■- ' 35F3 1:^ 

>■; '.t-tis- r' ; ; 

'.'07 0: 2.41 

E'"C6l7'43 

'E 08 
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;- : :3i6 v- 

**Q t /9 i«42 

E 061.01 

E 06 

, 


>' . . •■ 35-5 

. ' '3 !?• .7~ 

’> . ii * • ■ * f » 2 3 

t osb. isli 

E' 05“ 


• 

” Si 1 

■ : 318 ' 

- 3 , 62 . 

E 054.0 

E 05 . 

» * t* 


_ _ . .£ J 


~'Ti7~ . 1752 

c, • 0 * 0 

w J? » 

- 

• T 


! ' 3 23' '• 

1 a D ** ** * 6*3 x 

c 341.1 _ 

E 05 

* * f 



■ j 
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APPENDIX B 

FLEXIBLE AIRCRAFT EQUATIONS OF MOTION 


The flexible equations of motion of the TIPS aircraft have been com- 
puted at two flight conditons, climb and cruise, for symmetric motion only. 

Time prohibited computation at the third flight condition. 


The form of the equations of motion are as follows: 

Ax +■ Bk + Cbi- = O 

where 

x' *[/*,**-> <f, e , q, , , q 3 , , r? 7 ] 

A, B, C are square matrices of dimension 11 x 11 

U = component of inertial velocity in X direction, m/sec 

-us = component of inertial velocity in Z direction, m/sec 

^ = pitch rate, rad/sec 

6 - pitch attitude, rad 

r?- are structural mode deflections of the aircraft 

The seven structural modes are those modes discussed in Section 2.3 
of this report in order of increasing vibration frequency. Sensor equations 
have not yet been formulated. The equations of motion were computed using the 
Residual Elastic formulation in which the differential equations are corrected 
for the static effects of higher order structural modes above the seventh mode 
The corrected equations are included in Case 2, Climb and Case 3, Cruise. 


The set of stability and control derivatives of the rigid body de- 
grees of freedom corresponding to the above equations of motion are labeled^as^ 
Revised Modes, Cases 1, 2 and 3. The column labeled "elastic increment" repre- 
sents the static elastic correction for the structural modes not modeled 
dynamically in the equations of motion. 
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TIFS CASE 2, CLIMB 


no 

o 

•js. 




ppint cut of the sysiin "matrix *»»** 

NOTF - DENOMINATORS tfAVfc F-EM CLEARFO 

numbY°’ OF Tohf" IM SVFT£M= 11 
. NUMBEF. OF COLUMNS (TOTAL). = _ 11 





. .THERE A c £ 11 VARIfi D L£5 IK THE SYSTEM 



— 




UL(Trp u icc nilTD11T UAOTABI F5 







. AND _ 0 ARE INPUT. VARIABLES . . 



_ 



. _ -JHE 

ORDERING OF. THE VARIABLES WITHIN THE SYSTEM IS 

AS FOLLOWS 






♦ ¥ m 9 * nilTPUT WFPTAPlF<i * * * * * 




OflCtTTnil NIMF 

_ .__T_¥PE_ , . 

. .. POSITION - NAME . TYPE 

POSITION NAME 

_ _ TYGE 

1 

If _ _ 

. JURE LANE. H0TT0N.. 

2. . H. AIRPLANE hdtiok 

3 

~ Q . 

AIRPLANE MOTION- 




F IIF1 ..AIRPLANE., MOT TON 

. 

IIEP. 

ATFPLANF MOTION 

t 

7 

1 «ri 

... AIRPLANE. MOUO N 

8 UE<t - AIRPLANE .MOTION. 

- . 9 

.UE5. 

A I = °LA HE— MOTION — 

i 10 

UEfi.— 

AI RELANE .JSOTiatL. 

11 ..UE7. . AIRPLANE .MQIION . 




--- 


V 




h 



; TIFS CASE 2, CLIMB 


M ATP T^ir TTtLF OnFFFTCTEHt. HilB 

*. 'STORAGE" — IN CORE .• .• 

^ paRT r . r*,--- ■ ttr . 2 . . 7 - 

2 ■ ■ • - 5 1 i i i ft « : ii 10 ? 

'-. 15 g 4 F 4 ftfc - . 163 gE 416 > 2 MAEtfl 6 _ 7 x 760 . 46+02 . L ix 3 afiflE 3 :fl 3 - . 73 S 4 E . 4 . C 3 x 4 . 3 A. 7 E 4 03 . - . 3 22 fi£ l fl 3 » ^ 5 E t a ^ .... 

r-2 '. 22576*85 -. 1949 E 4 U 7 .1960E405 .3fl39E*03 .54926404 .4S87E404 -.7142E+04 .2306E*04 --.1681E405 4^44156+05-7 

.MWM8 '' -657I.F+BA 8. 7x4fi85£+04 .16246+05 -. .3985EJ Q5 -- .29 84E+.05 -..7J541EAJ.4 -. 75626*05 -'- '•.‘jWJPWi 


— itt? jyrH*?-: k » - 


*'4 • 0.. ~ 0. -.10C0E+01 0. 

|l5 19611 TF4Bg . 3B33E4Q3 .93296+01 fix- 

1.-6 « 2659E+5Z .12676+03 .9156E+83 B. 

r T — T47flF+<n -.jgnfiE.+na — -.1S33F+B4 — (U_ 
• -.3899E+M — #83472402 -.1928E404 0. 

■ 9 #19816482 -.711BE+12 .22386+1)3 D. 

18 .162BE+01 .23382402 -.2163E403 0. 


0 « 0 » 0* fl* , ‘0* 0* 4 

.7456F+n3 .79402402 -28736*03 -.64296402 .UB4E403 - . 2373Fttl3 / ' ' . fl iT.BF+ tfc- 

« 1307E401 .17536*04 -.17796*02 -.7747E402 .35526*02 -. 2113E*03 v , t l bJi , 4’383E40T 

.6993EtH3_ T.S418E4A2 2303iBE*fl4 .J.3ZZ6EQ3L _t.J. 1Q 864.02 .Zgfl9£*_rU— ?■ ' T . -? 1 ^9F . tn 3-- 

.Z045E402 -.3867E402 .1996E+03 .3930E404 .3618E40Z . 1385E403 , ir £yl2119E*Q3. > . 

.67 496*01 .26386*0? .1864F*n7 -. 75246*0? tfi!t&ZE±a4 -.,92 Q4.E+.0.2 .. jffl>y39 ) 7F*Q3._ 

.27246401 .4502E+01 .14256*02 -.1934E40Z .1H2SE4B2 .17986+05 S&«l*HE4Bj - 


TIFS CASE 2 , CLIMB 


STORAGE 


I.A ■ 0.-.- 
I D 

KIO * - C. ■' 


— IN COKE 

— 11.. It — 1 

— 1 OF 


S_»* 01TER11._ 


g 3 4 5 6 7 .* 9 10 11 

467 26+111 -- 233 IF..H? 0 . 3 fi 24 EJ-.fi 2 . 17 j 94 £.tAZ . 14 .filE. 4 fiJ. xlOJZE+OIL. _=x& 13 . 7 E.tJU .JtSiaEULl - , r . 2 B 97 F+ l) g 

. 24506*05 -.17796404 0 . . 2963 E 4 D 3 -. 12486+03 -. 23856*03 -. 2934 E +02 -. 3464 E 402 . 4113 E 403 . -. 1230 E 404 

-?B-j 7 F+n 4 . 7923 F +116 fi . _ 7 L. 19 . 41 E 4 .fll- . 1425 E+J 14 . _.=.« J. 37 . 7 E 4 H 4 — - xl 7 . 41£+ 04 =. 4 JL 26 E±J 13 - « 7 Qfi . 8 E . t .fl . 3 r-« gfl 796+04 .. 

a. o. . 1000 E 401 0 . o. o. 0 . o. o. ®« ’ 

-. 4 ? 31 F.B 7 -.? 3 f( 3 F*n? fi. . 12 B 1 F 4112 -? 677 F+(lt -. 2797 F .-01 . 53 . 9 JE.t-lU :, » 4663 E + .flJ . 37 1 5 E+IU ;r . 9 ft 31 F 4 tl l 

I' 2 Q 21 E 4 fil — . 7105 E+ 01 0 . . 38226 * 01 . 69866+01 -. 3795 E 401 -,Z 291 £* 01 . 2968 E 401 . 7605 E 400 348 6 E 401 

TB.tF.fl 4 -A? 9 ? F *ni r . 7 L. 1212 E 4 HI .--. 4 . 918 E +91 . 10 S 4 Etfi 2 . 9014 .EJJ 11 .-^ 2242 EJ 0 J =j t 832 EtC l. ~ .., v: . fi 4 fl flE t 0 1 — 

-. 6681 E 4 Q 1 . 93 SflE 401 0 . . 3624 E 401 -. 3936 E 401 . 571 SE 401 . 2 fi 8 lE 402 . 6948 E 400 -. 165 fiE 401 479 O 2 E 40 B 

— imoF.in — 4 & 3 X E 4 J 1 n. . 4257 EtDi . 298 BEJ 2 J — 2 S 81 £tOJ . lUlfi E lfi l. . 1312 £i 0.2 . 5633 E tl l II Aj '~ . ~fi 747 E 4 nri . 

-.16936401 — . 2554 E+ 01 0 . . 8014 E 400 127426400 -.93166400 -. 8421 E , 400 -.24886400 .. 1111 E 402 'l'-» 1734 E 401 ■ 


• .'‘STORAGE - 
jTr'ylATRIX.SI 

:• part — — 


or..— a.— 
■ i: o. - 

,3~;; o.~_ 
4" o. 


IH CORE 


TIFS CASE 2, CLIMB 


l^’.l v» 


4 

.. 0 . 

0. 

“o. 


• 4 'j'aS'S-V 

• 

' 5 6 7 8 9 10 M 

. 10 75E4D 0 - x812fl.E=fi J ..1197-E-01 ..4.372 EtlQJ — zs3Bb 3E-I12 _.„. 360 . 3 E- g l_ i 

. 1343E402 -• 8637E40 1 .4139E401 .31S5E401 -.2656E400 .4656E401 ' 

. 2695E4, 02„ -.3259E4.02 ..5.31.9E *-02—^.22426111 «.53&8E*fll... .. a1 6 . 86 E 10 2..; 

, 0. 0. 0. 0. 0. ; i 


• 6357E-01 .19506*01 ^6413E-01 -.1275E-02 -.13276-01 .28006-01' i 

.5662E-Q1 .687ZEtHJ .1823Eiai ..393?Er 01 . 4605670 1 — j . 21 . 25 , E-r Ol-j 

•1022E-01 .SSO7E-01 -.1336E400 .2083E*01 , .8183E-01 -.1164E-01 i 

,2033Erll — 13 15F-01 . 1131E-02 .J.1-4ZE-J2 iJOliE+ll ^SitaiErllZj 

.4029E-02 -.25996-02 -.5905E-02 -.51B5E-02 -.23S8E-01 .2018E401 , 

a hcoCinA . BfH7F-(H .Q7t*^F-n? -a?U7^F-H1 .?242E-Q1 1 





206 


TIPS CASE 3, CRUISE 


. print out of the system matrix »»»*■ 

_ H0TE.--L-DE1I0HINAT0RS HAVE BEEN CLEAREO „ 






NUMBER 0^ 

ROWS IN SYSTEM^ 

11 









NUMBER OF COLUMNS (TOTAL) = 

11 









THERE ARE 11 

VARIABLES IN THE 

SYSTEM 









OF WHICH li 

ARE OUTPUT VARIABLES 









ANO 0 

APE INPUT VARIABLES 








THE 

ORDERING OF THE VAPIABIES WITHIN THE SYSTEM IS 

AS FOLLOWS 







* * * * * OUTPUT VARIABLES * 

* * * * 





POSITION 

NAME 

TYPE 

POSITION NAME 

TYPE 


POSITION 

NAME 

TYPE 

i 

U 

AIRPLANE 

MOTION 

2 W 

AIRPLANE MOTION 


3 

Q 

AIRPLANE 

MOTION 

4 

theta 

AIRPLANE 

MOTION 

5 UE1 

AIRPLANE MOTION 


6 

UE2 

AIRPLANE 

MOTION 

■V 7 

UE3 

AIRPLANE 

MOTION 

8 UE4 

AIRPLANE MOTION 


q 

UE5 

AIRPLANE 

MOTION 

10 

UE6 

AIRPLANE 

MOTION 

11 UE7' 

, A IRPLANE MOTION 







t'if : MATRIX;TZJi.E'~- 

r*!' STORAGE-'"— — 0 

<•'•■>, MATRIX -SIZE -- 

[;.A.>panT -- — 

i;Cl'.;S.H-;i074E+03 .. 

.'38S6EVn» .-•■ 
'.'■i3'V-.U32E4B4 

1-fS ;t..'-«5S39E+H2 ' 
ft- G^_. .9G20E±gl_: 
'J-7' .•;• '.«1435E*fl2 ■ 

i :<t ■ .7590E4Q1 

-lH' ~ -8422E4tll 
iv 11.'. - ..954624«7- 


i TIFS CASE 3, CRUISE 

•'.-it* U' 




«' g 2 <r '. 

••2362E403 

.1(632485-, 

Uj- Zl 

,.514SE4D3,..v 
-.2229240 3.; 

-i2E27E+03 . • 


’’ 3 ' '■ 

^•3212E405 
- ,3 g 4 9E*0X- 
V-V 11652487 

.-.55S6E402 

. -.28862404 

364 22* 34, 

• 40042+03 
-.33842*03 
-.2794E+03 


.24222+06 -.49002+02 

_Zl 2A.33£*84 t lA9nE±H4_ 

0. -.15342+05 

. Pt Oj.-.— 

0. .79812+03 

JL, .2Q86E*01 

0. .3184E+02 

Jl. ..694.8£4°2- 

0. .2117E+02 

JL. .75.172+01. 

0. . -.2820E+02 


-.1248E+D3 -.4456E+03 .7505E+02 -.1744E+Q3 .3807E4Q3 --.2038E+04 

.1708E+8 5 .1372 E+0 5 -.227 6E+C5 . 7 193E4-0A- -« g£95E4flg. — :; -.17S8 E 4C6 -. 

.5283E+05 -.1356E+08 -.1024E+06 -.2517E+35 -.2464E+06 ,^;..6B33£405 


.5283E+05 

0. . — 

•Z439E4Q3 

_ «1.959EtBA.. 

-.1823E403 

-jlSAfcEAM 

•8585E+02 

. _»JL273E±52_ 
-.1459E403 


•6343E+03 
-.7395E+02 
.3314E4C4 
• 6777E+03 . 
.62Z2E+C2 


•.1894E403 

-.2562E+03 

.45832+03 

_t.4.315£444. 

-.8459E+02 


.3447E+0J 

.108GE4Q3 

-.3022E+02 

..i2L9E+J3__ 

.6S88E+04 


-.26712+03 .145QE+03 -.12B2E403 


-•7216E+03 .;-.,'j436 386+04 
-.6789E403 .1359E404 

.773DE403 tT«69V5E483 

-jlSPZ. 6 E+03_ ,S722Etg3.. 

3P2SE4B 3 - -f .'il83,0E*J,4 
1 4257E+0 1" l\'“ i2l73E409 


TIFS CASE 3, CRUISE 


■'""'"MATRIX TITLE — COEFFICIEHT MATRIX FOR S ** 01TERH 

•l" STORAGE IH CORE . 

MATRIX SIK — 11 X 11 


PART 


_2£_ _J5 


' -9 -.0. 


1 2 3 4 5 6 

.2470E+05 -,5438E40i -.2629E+02 0. ^7lG55E+02 -.3984E401 

. — 2132E+04 D. ._,5j)i6£+03. -J.1S05E+Q3 

, ’.29D6E404 . 7500E+06 0. -.3455E+03 .2S39E404 

, 0. P. .lOfiQJE+Ol .0, ... ,_J. .... 

, -.134SE+02 -. 2503E+02 0. .2040E+02 .3661E+ul 

, . iftLSF+tn -.9876E+01 o. ifiiZSki.ttl . . 9 . 8592+ 0.1. 

. -.3954E+Q1 Tl315E402 0. -.2092E401 -.8732E401 

-.75P0P+01 . 1505E+02__._C.«_ .62.97E+01 _-;*.7367£_+ai 

. ^14^E401 6054E401 0, .8S80E401 .4911E401 

- ' / - .17.55E +01 — r«.3j23EEt01 — 5. 41216E401 .JJ1JE1U. 

.8063E+01 • 1727E402 0. -.7440E+01 .6877E+00 


-.1758E+01 -.37S6E401 

-.A659E±fi,3 — =jl8Ji3£±flZ. 
-•2129E+04 -.330 3E + 04 

. 0 . .. « . . . 0 . . - - - - . 

-.6318E+C1 .10482402 

-.6998E4 Q1 -. 4 1 4 1E.4jQ.l- 

•1511E402 .15972402 

_ .9024E4O1 JL333EE+02. 

-.5177E401 .2091E401 

-.1 651E.4C1. *J^2.a£+Dl- 

.7936E4C1 -.53S6E4S1 


-.13812401 
.. S.672£402__ 
-.5429E403 

’ T75S9E401~ 
.5Q78E4Q1 — 
-.40082401 
. ...8.13&E+Q0_. 
•1792E402 


.5160E+01 

...87-76£*0-3. 

•1771E404 


-;i249E402 

..9RP8F4n4 

-•3319E404 


.8288E401 • ..-.2562E402 
■28B8E4P1 £ >E r 8330F4Dl 
-.5345E401 . - . 1260E4O2 

- . 4W39F4D1 . .44412400 

.1869E4D1 ^>-.‘43202401 


..1392E401 -,1219E402 |“j. .-.5293E4CI2 


TIFS CASE 3,' CRUISE 


” MATRIX TITLE - -- COEFFICIEMT MATRIX TOR S ** 02TERM 

i ,.v-L: STO RAGE IM COSE I : 

, ■- '.MATRIX SIZE — 11 X- 11 


r- 1 - o». 

ZJt Bj 

. .a o. - 

FTs. ■ o. 


■9S o. 
■ ■ 0 .- 


ay9.'ta« 

S*.# 1 ' n. 



-.1151E400 

_r.l-412E+02 

-.2372E402 

. 0 .__ 

•1841E401 

-*68fi2!--01 

.6150E-01 

«4.483E-D2_, 

-•3115E-01 

_TLM.63Er.Q2. 

.i611E4C0 


6 Z a 

-. 8754E-01 . 1060E-Q1 .4920E-B1 

. -l 9655E481 „4752£t01 l3.865E.*1L1_ 

-.35722482 .6233E402 .1163E401 

.o. a. a.. 

-.8637E-01 -.6815E-01 .53372-01 

■ 1938 E 4Q1 «7971E.r0 1^ .5767ET8e._ 

♦8578E-01 .18612401 .2268E-01 

75162 -UJL__ 17.5 t D 0 jZa59E.HU-. 

-« 2111E-Q1 .5237E-02 .3080E-02 

r«.5845E.r02 -.3753E -I1 7 -.31122-82 

.86412-01 .1077E-01 -.2733E-01 


-.50312-02 .3788E-01 

.^.MUZEiOD l5IL48£4JU-. 

•2901E401 .1932E402 

0 . 8 ■ 

-.36572-01 .4159E-01 

-.1941 2-01 . 34.79 Er01_ 

•5Q26E-01 -.3136E-01 

jZ366£r.0i =l2Z10.E=01_ 

•2007E401 .8249E-02 

. 2206E-0 S .2017E4P1 

. 26012-01 . -.3157 E-Ill 


-.1702E4C0 

-.16202402 

-•7309E+C1 


*'— • 9C07E-02 
..f-;3437e-Bl" 
, j 2183E4 01 " 
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TIPS CASE 1 , 


,-REtfISED-HODES ' 


- ' 



■ 

3L. . .. * 




* STATIC STABILITY DERIVATIVES 

* 

* 

¥ 



+**+*+*♦+*«««**+#**##***+*»««*## 


RIGID 

ELASTIC 

runsEMF wr 


TOTAL 

UNITS 

• 

CL ( l ) 

CD.JUJ 

. 326523 
xa.5?3fi 

-.007895 
-.nnnn a 


.318629 
. nnqi ?u 

• 

CM < u ) 

• i 69167 

-.008027 


.061160 


CLUI) 

... . 1.99*973- 

. . _ -_^cna623 


1 . 982160 

M / P AO 1 

C D < U J 

— -CM.1U1 

.263675 

-.00130 5 
-.ni.n?3A, 


.282370 
fl Qflft qfl 

(i/RAD) 
M /Q Am 







CLULPHA) 
CD (ALPHA) 

.10 J 5 86 
« v Zj.4x.fl2 

.003229 

.. 000022- _ 


. 100615 
-.-02010.9 . 

(1/DEG) 
(i /nPR» 

CM(ALFHA) 

- • - 31751 

-.300485 


-.032236 

(1/DEG) 

CL (Q) 
CD (Q ) 
CM IQ) 

15.3677 35. 
1. 351531 
-39 . 52*6 9Q 

--^L850X8. 

-.0 22846 
^718.435.. 

-_L5. 2022X7 .. 

1.328685 
-38^8 05656 

(i/Rfln> 

(1/RA0) 

(1 /RAD) 







CL(DE) 
CD IDE ) 

. -14656 
.1111550 

- .0 00 80 5 

. - -.0 0009.8 


.014351 
--.0 0i45Z. 

(1/DEG) 
tlZQEGi. . _ 

CM ( DE ) 

“.-57256 

.003716 


-.053541 

(1/DEG) 

- .CL (IHF.LLS.T 1. .. 

— ll 1 

“.n nnnnn 


_ nnnniin 

M/NPKTflN) 

CO(THPUST) 
CM’( THRUST ) 

-.100 l'C4 

• L L ii t U 1. 

-.000000 

..oaoooo .... 


-.003004 
. .OQoaoi. ... 

(1/ NEKTON) 
-(X/.NEHTQN) 
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/TIFS CAS E I , LANDING-CONFIGURATION. -REl/.ISED-HflOES. 


****** ** ************ ************* 


- ’ 


* 

DYNAMIC STABILITY 

. i s 

DERIVATIVES ♦ 

. . * 





** **************** *************** 




RIGID 

nnMPUTFo 


■ RIGID 
OUTPUT 

ELASTIC 

INCREMENT 

TOTAL 

UNITS 


cl <im 

1 .990973 
. 783875 


1.990973 

. 783675 

.009897 

.noi ?40 

2.000970 

.784916 

(1/RAO) 

(1/RAO) 

CM (U) 

-.058454 


-.058454 

-.021495 

-.079949 

(1/RAD) 

CL (ALPHA) 

CO (ALPHA) - . 

CM (ALPHA) 

. 100566 

.oaacMi2 

- • 0 31751 


• 10 05 86 

♦CC0884’ 

.(tr.pii? 

.101470 

,0711194 

(1/DEG) 

(1/flFG) 


-.031751 

-.001832 

-.033582 

(1/OEG) 

CL (Q) 

CO IQ). - 
CM IQ) 

15.387735 
1 * 351531 - 


15.387735 
1 .351531 

.322058 
_ .Q39898 

15.709793 

1.391479 

(1/RAD) 

(1/RAD) 

-39.524090 


-39.524090 

-1.250171 

-40.774261 

(1/RAD) 

CY (D£T A) 

-C-LIBEIA) — 

CH (OLT A) 

. oooaoo 

• OQO flO Q-. . 


•000 J CO 

- ft .1 a n r. i 

• 0 30 QCi 0 

»{) Oflfjfl ^ 

•ooaooo 

(1/DEG) 

(i/nFG) 

.000003 


• 0 3 0 j CO 

.000000 

.oooooo 

(1/OEG) 

^Y(P) 
C1(P) - 

CN (P) 

.000 BOO 
. OS 003 0 . 
.000000 


• 03 C 0 CO 
•asocial*. 

•CU0303 

+ C v 0 n n n .. 

.003000 

.nnnnno 

(l/RAO) 

(l/RADI 


.000000 

.000000 

.000000 

(1/RAD) 

CY (P.) 
Cl(R) 
CN (R) 

,OQ0(>0O 
. 0 0(1 300 


.DOuOOO 

.oaoyoo . 

» C 0 G u 0 0 

. . .iunflii a 

.000300 

.oooooo 

(l/RAO) 

.000030 


.003000 

.000000 

.oooooo 

(1/RAD) 

CL(A-DOT) 
CD lA-firtT) 

-2 .734893 
- .404575 


-2.734693 

• .012566 

.0 LI 20 6 

-2.722327 

-.403322 

(l/RAO) 

11/RAD) 

CM (A-DOT) 

-14 , 0*. 075 1 


-1** .0407 51 

-.261637 

-14.302387 

(l/RAO) 

CL (Q-OOT) 
CQ(Q-QQT) 
CM (Q-DOT) 

-17.115211 
-2 » 3*216-8 
-£1 .069772 


-17.115211 
t£. 3421.68 

-.387189 

-*U485E2 

-17.502400 

-7.399730 

(l/RAO) 

(l/RAD) 


-21.069772 

1.203445 

-19.866327 

(1/RAD) 

CY (B-OOT) 
CK9-00T) 

.000000 

OOCOQO . . 


. 00QCC3 
Q Q 

.oooooo 

.0 00 3 0 0 

.000000 

.OOOOOO 

(l/RAD) 
(l/RAO) 

CN (B-DQT) 

.000030 


.000003 

.000000 

.000000 

(l/RAD) 

C Y (P-0 OT) 

.000003 


.000000 

.Doaaoj. 

.00000 U 
.000000 

.OOOOOO 
.000000 

(l/RAD) 

(l/RAD) 

CN(P-OOT) 

.oooooo 


.000000 

.oooooo 

.003000 

(l/RAO) 

• CY (R-OOT) 
Cl (R-QOT) 

.000030 
, .UQ00D& 


.000000 

. . jlMiOO 

.oooooo 

.oooooo 

.000000 
.oooooo. 

(l/RAD) 

(l/RAD) 

uN(R-OOT) 

.008000 


.030003 

.oooooo 

.000300 

(l/RAO) 


'1 




ORIGINAL PAGE IS 



_\T I F S CASE 1 , _U.NDING.XQNFI(iliRATIQN,. REV IS ED ..NODES 


*** ******* 

* _ _ ... *_ 

* ACTIVE CONTROL DERIVATIVES * 

*** **************** ***** ****** 


* SYMMETRICALLY DEFLECTED CONTROLS * 


RIGID RIGID 

.-C-QM.RUI.E_0. XUTEUX .. 


ELASTIC 

.INCREMENT 


TOTAL 


UNITS 


CL (OS) 

-ca.(QS). 

CM (OS) 

-uL.lDS~.DOXI- 

CD(DS-OOT) 


.012626 

..U-D XT 6 8 — 

-.005329 
--X.6695.83 . - 
-.221476 


.012626 
..^0017 59 
-.005329 
rl. 6635 83 
-.221476 


■ . 0 00 06 9 
• • U-x0 ill 9 
’ . 0 vl) 0 36 
=*022295 
■ . u C2 82u 


.012557 
.000750 . 
-.005365 
-1. 631.6 7 6 
-.224298 


CL (DS) 


CM (OS) 


CONIRQJ 2_- DLE 

.023603 .023603 -.0G0197 

... 0 0178 2 .0 0X7X2 r -.fl.SQj] 26— 

-.003323 -.C03323 .000376 


.023407 
- . M1Z5Z. 
-.G02947 


CD (OS-OOT) 


CL (OS) 


•.360244 
=..3197.62 


.014856 


-.380244 

-.313732. 

CONTROL 3„..r_.ELV_I 

• ii!4856 


.003571 
1.15662 8 t 


-.000144 


-.376673 

•lxCiftkia 


.014712 


(1/OEG) 


CM (OS) 


■ C 57258 


-.057258 


.UGQ58G 


-.056678 


(1/OEG) 


CO (OS-DOT) 


-.151547 


-.150547 


-.004290 


-.154837 

2.202565 


(1/RAD) 

tl/RAD) 


CL (OS) 

CD (DS) 

CM (OS) 

_CJJUCLSx.n.DIL. 
: CD (DS-DOT) 


.020792 

Jl0.2.1iL9L 

.077093 

.260639 

,169031 


CO NTROL 4_ 

.020792 

.002 1 09 

-.077093 

^U-?6M.33- 

-.169331 

1.549343 


.000716 
.iP.00063. 
-.0 02720 
•.9 48898. 
-.006093 


.021508 

_. .M.219.8- 
-.079813 
-1.309537 
-.175924 


(1/DEG) 

(1/RAD) 

(1/RAD) 
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T-rce pacc 9 n tmr. REVISED MODES 


********************************* 

* * 

* STATIC STABILITY DERIVATIVES * 

* * 

******************************** 



RIGID 

ELASTIC 

TOTAL 

UNITS 



INCREMENT 




cl ( ; ) 

. 32942 8 

-.010397 

.319331 


onto) 

. .01 5T?3 - 

. - . Of Q.iit8 

_ . .... 

- 

CM ( n ) 

.070229 

-.010361 

.059869 


r.niji 

1. 554462 

-. J15283 

1.539179 . 

(1/RAD) 

coej) 

. '.532=3 

-. 1G15A9 

.156705 

(1/RAD) 

CM(U) 

-,05?73U 

. -.IlBlSu .. 

__t.Q7.1124. 

..(1/AA01 

OL< ALPHA) 

. 102325 

. 0*’3304 

.102630 

tl/DEG) 

Tl ( ALPH 6 ) 

.014888 

.036024 

.014912 

(1/DEG) 

C M ( AL C H A ) 

-.0 318=4 

-. 3 v. i 647 - 

-.032501 

( 1/OEG) 

CUD) 

15.524966 

241726 

15.283240 

(1/RAD) 

cn ro) 

. 938506 

-.019707 

. 918799 

(1/RAD) 

Cm(0> 

-39.797429 

.929954 

-38,867475 

(1/RAO) 

r L ("OF ) 

• 015 004 

“ ":.oTiosr “ 

.013953 

(l/DLG) 

cntoi) 

. 001 058 

-.000084 

.006975 

(1/DEG) 

CM (OE) 

-.057701 

. OQ 4 823 

-.052876 

(1/DEG) 

CL t THPl)?T) 

WLOlC 

0030 Oil 

. .__..QjJ<liDLG 

( 1/NEHT ON) 

Cn(THRUST) 

-.O0COQ3 

-.300300 

-.000003 

(1/NEWTON) 

CM (THRUST) 

• 0 COOQ 1 

. 00 0000 

.008001 

(1/NEWTON) 



quality 



212 




j TIFS CASE 2, CLIMB, 

REVISED MODES 


. • * -*.■ ’ 1 


' ■ ?**«*.**** 4* ** 4 ****** 5g 

' '• - 

> ■ •*» 

OY N A MIp'- STA BIUT Yi 0 ERI V AT IV ES ' ^>.-4 fe 

’V -■ "• V. ' *•*}.- s 5,?/,, a a; 

J * , f 

,, 1 , 1 ' 

V 5 ,, -.JSi}bA^W~X'n 


' RIGID 
COMPUTED 

. RIGID ’• 

output’ . 

: ELASTIC., 




* / 

j* v 

. — * 

x*;* „,;i . v, ; *■.*&.. 

CL(U) 

coun . 

1. 554462 
.158253 

1.554462 

.158253 

.0070 26 . 
.008579 ' 

•1.561488 vV(l/RAO)' 
- .^58832 " 4 ' -<l/RftD> -: 

CH(U) 

-.032734 

-.032734 

-.015076 •' 

;-i 047610:',; ' • (.l/RAD) /; 

CLtALPHA) 

CO(ALPHA) 

.102325 

.014888 

.102325. 

.814888 

.001173 

.880049 LJJ. A. 

\il0349$’ '(l/OEG) 
.014988 : 11/DEG1... 

CM(ALPHA) 

-.031854 

-.031654 

-.002426 ■ 

-.034280 " (1/DEG) 

CL (Q) 
CQ< (?) 

15.524566 

» 9.12.5.16. 

-39.757429 

15.524966 

^3385.06. 

.429217 

j.0351.37_ 

15.954183 tl/RAD) 

CH(Q> 

-39.797429 

-1.659416 

-41.456845 (1/RAD) 

C Y ( BETA) 
Ci(BETA) 

” .000000 
. onoonn 

.000000 

.000000 

.000000 

.008000 

.000000 11/OEG) 

.-6MM0 (1/OEG) . 

CN (BETA) 

.000000 

.000000 

.000000 

.000000 (1/DEG) 

_______ — 

Ci ( P) 

'.do'ob’ijft 
,0 00.000 

.000000 

.000000 

.000000 

.JLMMO 

.000000 H/rM) 

.000000.. _(l/RAO) — 

CN(P) 

.000000 

.000000 

.000000 

.003000 (1/RAD) 

CY<**) 

Cl( 

•000000 
, ooqooo 

.000000 

.cnnnoo 

.000000 

.000000 

.000008 (1/RAD) 

.oooonn ( 1 /rao) 

CN(R1 

.000000 

.000000 

• .000000 * 

.000000 (1/RAp) 

CL(A-OOT) 
CD( A-DOT» 

-2.962224 

-.312355 

-2.982224 
-.312355 ■ 

,016536 

.000864 

-2.965688 (1/RADI 

-...311472 U/RAP). 

■ CM ( A-DOT) 

-14.290997 

-14.290997 

-.350494 

-14.641491 (1/RAO) 

CUQ-OOT) 
CO(O-DOT) . 

-17.890330 

-1.694267 

-17.890330 

-1.69.4262 

-.525633 
-.843827 . 

-18.415964 (1/RAO) 

CM(Q-DOT) 

-20.940744 

-20.940744 

1.629167 

-19.311577 (1/RAD) - 

3 

CY(B-OOT) 
■ CUB-DOT) 

.000300 

.nnnono 

.000000 

.000000 

.000000 

,-JHLOMQ : 

.000000 ■ tl/RAD) 

_,_0.Q O.G 0 0 MZRjS D) 

, CN(B-DOT) 

.000000 

.000000 

.000000 

.000000. ■ ll/RAD) 

CY(P-DOT) 
ct (P-nnn 

.000000 

•oonooo 

.000000 

•nnnnno 

.003000 

.009000 

.000800 ’ tl/RAD):, 
.Ttinnnn M/Pim 

I CN (P-OOT) 

.000000 

.000000 

.000000 

«>>'c:j3 (1/rad) 

CY(R-OOT) 
• ni (R-nnn 

. .000000 
.000000 

• .000000 
.000000 

.000000 

♦ ononoo ... 

.030530 (1/RAO)’ 

-nn 3.0 i usaqla. 

'CN(R-OOT) ■ 

.000000 

.000000 

.000000 

.uoaasa ci/radi 
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’ TIFS CASE 2, CLIMB, REVISED MODES 


l 


*********************** ******* 



— - — ■ 

"" 

¥ 

ACTIVE CONTROL DERIVATIVES * 





¥*¥¥*¥***¥♦*¥¥*¥*¥¥¥¥¥¥¥¥*¥¥¥¥ 




♦ 

SYMMETRICALLY DEFLECTED 

CONTROLS * 



’ NAME 

' RIGIO " 


RIGIO 

ELASTIC 

TOTAL 

UNITS 

— 

COMPilTEO . 


OUTPUT 

INCREMENT 






. . CONI R0.L 1 -_&IL3_._ 




CL(ns> 

.012741) 


.012740 

-.CC0C92 

.012649 

(1/OEG) 

CO(QS) 

* 0 C Q 56? 


.000565 

-.000008 

.000557 

(i/OEG) 

CM(OS) 

-.005358 


-.005358 

-. 00 i) J 49 

-. J054C8 

(1/OEG) 

CL(OSrl)Or) 

• -1.725481 


-1.725481 

-.030091 

-1.755571 

(l/RAD) 

co(os-dot> 

-.156273 


-.156273 

-.0)2541 

-.158815 

(1/RAO) 

-C.H(J.5-.DQT.l 

. . _1 ,.322956 . . 

. .. 

... 1.322656 . 

. _ .091488 

1,304144... 

Il/RAD) 




CONTROL 2 -■ OLF .. 




CL(ns) 

.023871 


.023831 

-.000261 

. 0 2357C 

(1/DEG) 

CD(OS) 

.001322 


.001322 

-.000023 

.001299 

(1/DEG) 

CM(OS) 

-.003354 


-.003354 

. 0035Q0 

-.002854 

(1/DEG) 

CL(OS-OOT) 

.-^2. 8.76501 


. ~2a 8711501. 

. ...J139240. . . . 

t2. 6.31261 

_(1/RAD)_ 

CD ( OS-DOT ) 

-.272916 


-.272916 

.003140 

-.269777 

(l/RAD) 

_CM<QS-DOTJ. 

. .-.967404 


-.967404. . 

-.209804 

-1.177209 

(l/RAD) 




CONTROL 3 - ELVT 




CL(OS) 

.015004 


.015004 

-.000192 

.014812 

(1/DEG) 

co<ns> 

.001058 


aOOi658_ . . . 

.. _-,MD316 _ . . 

_ ..0.0104.3 . 

.. .(.1/DEG1 . 

CM(DS) 

-.057701 


-.057701 

.003773 

-.056928 

(1/DEG) 

’.CLdlS-DOI) 

, -1.18348.9 


-1.183489 

-.046871 

-1.230360 

(l/RAD) 

C D ( OS-DOT ) 

-.106788 


-.106788 

-.003870 

-.110658 

(1/RAO) 

.CH(DS-QOT). . . . 

2.153215. . 

- 

2.153215 

. .168493 . 

2.321708 

(l/RAD) 




. CONTROL. ..4 - HT 




CL ( OS) 

.020983 


,020983 

.000954 

.021937 

(1/DEG) 

,_£LQJ[ DS.l 

.flJD 1 46.4 


_ .001464 

. 080378 ... 

.901543. 

(1/DEG). 

CMtOS) 

-.077642 


-.077642 

-.003611 

-.081253 

(1/DEG) 

_CL.LD.SrQnU 

-1 ..3 14917 ... 


.._ si. 314917. 

... s.,1) 66445 

-1.381362 . 

(i/RAD) 

, CDC OS-DOT) 

*.121715 


-.121715 

-.005486 

-.127202 

(l/RAD) 

' 0M(DS-.D.QL1 

1.61 8584 


1^618584.- . 

.,J239812_ . _ 

1.858396 

(l/RAD) 



Qg eoor 



j. 


[ TIFS CASE 3, -CRUISE ,-£NG4NE-I-N ~y-REVI560- MODES-- 


* STATIC STABILITY DERIVATIVES * 
********¥*«********?**¥*¥******* , 


- 

RIGID 

ELASTIC 
-- INCREMENT-- 

TOTAL 

UNITS 

CUO) 

CD (3) •- 

CH ( 0 ) 

.153030 

- -it 55 581 

.578810 

-.034849 
■ — -.486-474 
-.029248 

.318181 

.6055G-7-r — 

.049563 



CL ( U) 

CD (U) 

CM(U> - • - 

a 55838Q 
.• 13762 
— -.402574 

■ ' -.081902 
-.000181 
-.0-86716 

.476478- - 

..13580 
- -.009685 

(1/RAO) 
(1/RAD) 
•(•1/RAD) ■ 

CL < ALPHA) 
C"! ALPHA) 
CM( ALPHA) 

.115 243 
.003238 
-»C 32883 

.001063 

-.000012 

-.302350 

.111306 

.003226 

-.035233 

(1/DEG) 

(1/OtG) 

(1/DEG) 

- CL(Q>- • 
CD(D) 
CMtC) 

16.452161 
..21674 
-42.526824- - 

- -.738504 

. 0 060 8u 

2. 710491 - 

15.711656 • - 

.027754 
* *—35.316334 

(1/RAO) - 
(1/RAD) 
(1/RAD) - - - 

CL (DE) 
CD(DE) 
CM( DE) 

. .15891 
.,,;r.'i,T46-- - 
-,o 61 319 

-.003193 
. 0000 32 
.014343 

.012697 

-.800013 

-.046976 

(1/DEG) 

(1/DEG) 

(1/DEG) 

-CL ( THRUST) 
CD (THRUST) 
CM (THRUST)- 

-.ooo&oo 

-. -JCOO^l 
• 0PC0C0 

-.000800 * 
. 0 0 0 0 00 
.000000 

-.otoooc — . 
-.008001 
.000000 

-(i-Z-NEWTON)- - 
(1/NEWTON) 
(1/NEKTON) 



| TIFS 

CASE 3, -CRUISE ,ENGII 


t v'l V”’ - T 

::vi. 




* ?V* **• 

• 


«*»«»****##»**««**«**«*«#»«»««#«* ‘ 

. . , ,<i » > ,■ , J. 

' ‘ ./ ? ..** s . ' 





* 

* “« * 

' 1 j * ' * ‘ 



* DYNAMIC STABILITY 

DERIVATIVES * 












• • -A-, 



RIGID 

RIGID 

ELASTIC 

TOTAL 

Units- • 




TMfiHFMFNT 






t t 



CUD) 

.558380 

.558300 

-.025468 

.532832 

(l/RAD) 

■G6(J)) - - - 

- .3-13763--- 

.413663 

.663108-- 

-.0-1370 i n/RAOV— 

CM(U> 

-.002974 

.008767 

.058520 

.067287 

11/RAD) 

CL ( ALPHA) 

.110243 

.103323 

.004106 

.107429 

(1/DEG) 

CG( ALPHA) 

,003238 

.003211 

- -.000024- - - 

— - • - .00318-7-- 

- C-t/UEG) - 

CM( ALPHA) 

-.032883 

-.002379 

-.008656 

-.011035 

( 1/DEG) 

CL ( 0) 

16.450161 

16.450707 

1.562276 

• 18.012983 

(l/RAD) 

GUI 9> 

.021674 

.015748 

— -.013132- - 

-.002616 - 

(-1/-RAD)- 

CM(0> 

-42.026824 

-41.683927 

-6.G31389 

-47.715317 

(l/RAD) 

CYI95TA) 

. • 3 c 0 0 0 j 

.000000 

♦ Du 0000 

.OOOOOO 

(1/DEG) 

C1<4~T-A} 

. J .--.ongeoo-- 

.609069 - — 

* - 

.OOOOOO — 

— <4/©eG»- 

CN(95TA) 

. 030003 

« L DC C Ci} 

. oaoooi 

.000000 

(l/OEG) 

0 Y ( p ) 

.00GC00 

.OOOGCO 

.000000 

.001)000 

(l/RAD) 

Ci( D ) - 

,300003 

- — - - " — ** 

.OOOOOO 

- -- ---.OOOOOO 

— M/-RAQ-) — 

CN(P) 

, 0 CO QOS' 

.000000 

.000000 

.000000 

(1/RAO) 

CY(R) 

.oocooa 

.000000 

.COOOOO 

.003000 

(i/RAO) 

GMR)- - - 

-.030003 — 

- ,600600-- • 

.OOOOOO-- 

-.000080—- 

. -«y.RAD».- 

CN( °) 

.080000 

•COOOOO 

, .000000 

.000000 

(l/RAD) 

CL(A-DOT) 

-5 . 270920 

-5.070920 

.060618 

-5.010-301 

(l/RAD) 

GO-< A-DOT) 

- -.049114 - 

C69116 - - 

-.-.000404- 

-.061518 - 

- -(-1/RAO)- - 

CMt A-OOT) 

-16.460042 

-16.460042 

-1.378287 

-17.838329 

' (1/RAO) 

CL(G-DOT) 

-23.981482 

-23.981482 

-2.193535 

-26.175017 

(l/RAD) 

CD-tG-GOt) 

- •« • 1C1GB3 - 

. _ —-.101-0 83 --. 

. . . .015327 

- — - - -.085756- 

- (1/RAO) - 

CM(D-OOT) 

-19.811921 

-19.811921 

6.873396 

-12.941525 

(l/RAD) 


; C Y( 9-DOT) 
--GM 9-DOT) 
' CNCR-DOT) 


.occooo 
. 000003 — 
. oooooo 


.uococo .oooooe .oaoooo (i/rad> 

.000008 --.003000 - .OOOOOO tlVRAD)— 

.COOOOO ■ .DOJOOG .OOOOOO (1/RAO) 


. CY(P-OOT) .000000 ,003000 .000300 ,000000 (1/RAD) 

i„Gl-<n~eoT) .000008 .OOOOOO- .00-0000 .003000— (-l/RADl- 

CN( p -DOT) .000000 . 000000 . 030000. .GOOOQO (I'/RAD) 


'CY(R-DOT) .0 00000 .COOOOO . 000000 .OOOOOO - (1/RAD) 

-G-HR-DOT)- -• ,000000 • - • - .00000-3 -.003000 ,00(1000 -(1/RAD)- 

' CN(R-DOT) .OOOOOO .COOOOO .OOOOOO' .OOO'OOO (1/RAO) 
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— -jTIFS 


NAME 


CUDS) 

GD< DS-) 
CM(OS) 

6K OS-DOT) - 
CDf OS-nOT) 

c mt os-dot) 


CL(OS) 

CD(OS) 

CM(QS) 

-OL4 OS-GOT) 
CO ( OS-DOT) 
OH(-0S-GGT) 


CKOS) 

-GD4 OS-) 

CHIOS) 
-CL4-9S-O0T)- 
C 0 (OS-DOT ) 
--GM4-D-S-DOT) 


• CKOS) 
-C-G40S) 

< c m ( ns) 

-GL f DS-DQT4 
C Df OS-DOT) 
-C-MT-OS-DOT) 


CASE 3 , —cruise, eng i ne i-n-t- r e v ^se o -h od eg 


_j» « 

* ACTIVE CONTROL DERIVATIVES * 


* SYMMETRICALLY DEFLECTED CONTROLS * 





R IGIO RIGID ELASTIC TOTAL UNITS 

COMPUTED - - — OUTPUT- INGR£HE+fT 


v . 013731 
04) 5163*-* 
-,005797 
- - —2..1 57571 

'0 0 0 670 

■■ - 4,469189-- 


—CONTROL-— ! - A - I - L - S 1 

..013731 -.000318 ..01341-3 41/, DEG.) 

* - - O 0 163 — - -- — -—.-000000 — .3 00103 — ll/OEG ) — 

-.CO-5797 -.000204 -.006001 (i/OEG.) 

- - -—2,457571 -.1-19514 2-. 2 7-7085 (4/RAD4- 

-.000670 .0CJ682 .000012 (1/RAD) 

1-.-4 59 4 39- — - • .238863— 4,765654 44 /RAB-) - 


.025615 
.000 352 
-.003513 

— - -3.643225 

-.017502 

- -• -1.414440 


CONTROL- 2- - DL-F — 
.025615 

- --.005352 • — — 

-.003513 

— 3.543325 

-.017502 

- --1. 414440 


.000901 
.030002 - 

.001844 


, 1522 52 - 

.001598 
•-824275 


.024713 (1/DEG) 

• -.-305-354 (4/-0E64— 

-.001668 (1/DEG) 

-.019100 (1/RAO) 

2.238715 ( 1/RAD- T- 


.015 891 

.... -.800046- 

-..0 61319 

— — -1,536891- 

.000309 

--- 2 , 8043 - 17 - 


.022141 

— ,0 000-15- 

-.082134, 

—-1,732485 

-.003353 
2,2-15583- 


- ■ CONTROL -3 ELVT 

.015891 -.003720 

. .000046- ,0.;8M6 

-,061319 - .0,32889 

-1,536891 —.-200-245 

.000309 .001558 

2, 8040 17 ,7-24096— - 


.015170 (1/DEG ) 

,000 040 — (-1-/BE-64— 


-.058430 (1/DEG) 

-1.7-37135 (-1-/RAQ4— 

.001867 (1/RAD) 

3 ,52011 3 44 /RAO— 


GGN-TROL 4—--HT 

.022141 .003456 

,000015 --,06-00 2-8 

-.082134 -.013122 

-1,732485 ,283140- 

-.003353 .002180 

2,-215583- 1.-045986 


.025596 (1/DEG) 

,00 0013 4 I/O EG)— 

-.095256 (1/DEG) 

—2,01 262-5 (4-/RAG4- 

-.0lill73 (l/RAO) 

— 3,2-31669 (1VRAO )- 



origin® BBS* IS 

ca poor. 

APPENDIX C 

PRESSURE DISTRIBUTION DATA LANDING CONFIGURATION 

This appendix contains pressure distribution data computed by FLEXSTAB 
for the landing configuration. The coordinates given are of the aerodynamic 
panel centroids in the reference axis system. Panel numbers start on the 
leading edge and progress rearward on the thin bodies. The pressure coefficient 
Cp is defined as follows. 

^ _ -Piowe* ~ upper. 

» - J 

where and are lower surface and upper surface pressures 

§ is dynamic pressure. 
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_TIFS CASE 1,-j 


-MDOLES- 


************************************* 

J_._ *- 

♦ STEADY AIC PRESSURE DISTRIBUTION * 

Jt... — — 

************************************** 


* NUMBER OF SLENDER BODIES — — - 2 

♦ NUMBER OF THIN BODIES 7 

* MACH NUMBER — — — .23 

♦ ANGLE OF ATTACK 7.0227 

- (DEG) _ 

* SIuESLlP ANGLE “ 

(OEG) 

♦ DYNAMIC PRESSURE 2819.68 


. (NEW TONS 


DATA FOR SLENDER 300Y NUMBER 1 


* 

SLENDER BOu Y‘ 

NAM E 

FUSELAGE 


¥ 

Y-COORDINATE 

--------------- 

, — — . ^ ^ 

0000 

l METERS) -- 

♦ 

Z-JUORCI NATE 


, — — — * 

CC00 

( METERS) 



nf CVMMPTWY 

DM 



* 

NUM3ER OF DOUBLETS 


20 


NO. 

X 

CP(X) 

CP <Y ) 


CP(Z) 


l METERS 1 






4 


. ir. ir.o 

.ti i> b u ij j 


. 35240.8 

2 

7172a 

. 0 J 0 0 G 0 

.0 L000 u 


, 150461 

3 

3 .4£2d75 

. GQQQ 0 0 

• 0 00 0 0 Q 

_ _ . 

• b 17 4_4i 

A 

A . 334 l25 

.CJCOOC 

« C C C C 0 u 


.022176 

5 

b.215 175 

.cuoca 

a^DDuil 


.-^012.919 

6 

7.596325 

• u u 0 U 30 

• 0 000 jO 


.012294 

J 

A. 377475 

• fi n o i ? 3 

. o on no o 


.055169 _ 

' 8 

10.358625 

.006110 

.000000 


. 302103 

9 

11 .73^775 

^uJOOflQ 

. C i: fl il 3 l 


.098100 

1m 

13.120 925- 

.-J 0000 

.ooooco 


.113105 

ii 

14.5G2G75 

.000000 

.000000 


.116748 

1 ? 


".■•Mg' W i 1 1 IJ tfiTt 




13 

17.264375 

.000000 

.000000 


.062232 

14 

18. .645 525 

.0000110 

*0000011 


051194 . 

15 

20.026675 

.000000 

.000003 


.014893 

lb 

?1 .407325 . 

.nnoooo 

.0 0000(1 


r. 014801 ...— 

' 17 

22.788975 

.000000 

.0 30000 


-.053082 

i fl 

2A.17., 125 

-0.10000. 

.000000. 


1 1 1 II U L 

19 

25.551275 

.000000 

.0 00000 


-.152734 

2D 

.9321*25 

.000000 

. 0 0QGQ 0 . 


r« 1940 66 





- 

.. - * 
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TOGINAL PAGE IS 
HHE££>R QUALITY 


TIFS CASE 1 j . L A N CI NG CONFI GU R A T I ON , R E VISED ho pes 



DATA FOR SLENDER BODY 

NUMBER 2 






« 



f UCTCDC1 





i— 

ON/OFF PI ANF nr cvmmctov 











NO. 

X CP<X> 

- - .IMET&ftS) 

CP (Y ) 

CP(Z) 


4 -9.4 9 4.1.0 — — . 00 8470 ^.005354 .454995 


2 

9.924-00 

• 00 03 C 0 

-.005486 

.288405 

3- 

lu .354003 

.03-0885 

- , g 21Q39 

.-258457 

4 

iw »7 34 Q J 

.003000 

-.036605 

.235744 

• 5. 

- 

. . .853868 

- * i j 3661 9 

.577912 

6 

ii « 644 l. 0 v 

-.255242 

- ♦ fi 55411 

« 444494 

7 

1? tO 0 3 



.RQflPFn 

8 

*.<2 « 5 14 u 

-.187054 

- .063615 

.555473 

9 

13«934vJ'^- 

- --.11579a 

-.022563 

.563264 

I'- 

15 » 3 64 4 g u 

"• j8 >)6 33 

-.074722 

.561508 

ll 

13*79fu0u 

-.0h3455 

-.056098 

.294158 

_ 

-»^34 wv 3 - — 

- i V 1 1 v QR 

— - D ^7 *?') 

.?*? m 

13 

14.#654vC « 

-. -J7296 

-.071154 

.269764 

14 

ip • 3 u 

-^uj713u -- 

- a .551521 

.213768 

15 

lb • 51h U00 

-.014200 

-.239102 

1.540847 


TIFS CASE 3, LANDING CONFIGURATION, REVISED MODES 


DATA FOR THIN 8O0Y NUMBER 1 


* THIN 800 Y NAME 

* DIHEDRAL ANGLE -- 

* on/off plane of symmetry 

.. - NUMdER- OF. -PANELS.. 

* NUMBER OF Y-ST ATIONS 


VER TAIL 

9G.0GuQ„lDEG^l 

ON 

5 

3 


PANEL 

NO. NO. 


X 

IMETlRS) 


X/C 


CP CP . CP 

SYMMETRIC, .. . RIGKI_SID£ LEFT SIDE 


Y-ST ATION i <Y = 2.1803 (HETERS)) 




23.825982 .096621 

25.735J.Z2 ..599597 .... 



.000000 .000000 


Y- STATIQN 2._LY_=_ 3.7276 < MET ERS.11 


... 1 3. .. .2*»^4a!iilBfl ._.593Z0£ .UMU .00001111 .000000 

2 4 25.803215 .599249 .000000 .000000 .800000 

Y~-7miON~3 {V = T.2604 (HETERSU 


1 5 24.964816 .084403 .000000 .000000 .000000 

2 E 25 .2.70618 ..598118 .0000 00 ^UIMM : .000000 
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TIF$ CASE 1 » i L^WH W -CO HF IGURA W M-r -R E tfl SeO -HO O ES 


‘ n T£’"* M 5 Y “f 



DATA FOR THIN .BODY NUMBER 

2 






♦ 

1 nlN D UU I WSnt 

4. ft??A IDPG* 1 






S- 

NUHBER OF - PANELS -y — t — — 

18 



* NUMBER OF Y -STATIONS — 3 . „ 

PANEL X X/C CP CP CP 

NXU (METERS) = - - — SYMMETRIC R IG HT - S IDE L EFT STDE 


Y-STATION 1 (Y = 1.1290 (METERS) > 


1 7 11.939303 .049954 3.470485 

2 8 12.327974 . 14996-2 1.225009 

3 9 12.910960 .299974 1.005248 

4 10 13.688321 .499990 .670036 

5 11 14..465662 .700006 ,412171 

6 12 15 -. 24309 . 3 - .900222 . 2821(14 


3.470485 3.470485 

1.225C09 1.225009 

1.005248 1.005248 

.670036 .6700 36. 

.412171 .412171 

_.28Zifl-4.-. .287104 





-13- 


12.0225-33- 


>049900- 


14 

-15- - 
16 

.17. . 
18 


12.396591 
-12.S57679— . 
13.705795 

-1.4, .452911 

15.202027 


,149917 

.299944 

.499979 

—.7110 0.1. 4 

.900049 


3.8126 39 3.812639 3,81263-9. 

1.310160 1.310160 1.310160 


.669765 .669765 .669765 

.4168 63 .416663 .416863 

.309292 .309292 .309292 


Y-STATION 3 (Y = 2.6670 (METERS)) 


1 

19 

12.113334 

.049921 

4.191875 

4.191875 

4,191875 

z 

22 

12.47145U 

.14393 5___. 

JU24727X 

. 1.247277 

1.247277 

3 

21 

13.008625 

.299955 

. 988019 

.988019 

.966019 

.4 

. 2 Z 

13. 724898 

.49998 3 

.625773 

.625773 

.625773 


23 

I**. 441098 

.700011 

.384563 

.384563 

.384563 

6 

24 

15.157323 . 

.900039 

. .30J6956- . . . 

. -.346956 

.30 6956 
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drigxnau page is 

OS £OOR QUALITY 


TIFS CASE 1,L 


„REVISEQ. MODES. 




DATA FOR THIN BOOT 

NUMBER 

3 




» tutu ar 




““ ' ‘ 



- * nTHFngai thir.i f 


ft. snnn f npfc- 1 




♦ ON/OFF 

♦ NIIHRFP 

PLANE OF SYMMETRY 

OFF 

ZP 




* NUHBER 









4 


_NC.. . 

PANEL 
. NQ. 

X 

_ (HFTFRS1 

X/C 

CP CP 

SYMMFTRin RIGHT 5TnF 

CP 









. 

Y-STATION 1 (Y = 

4.9075 

(METERS)) 


1 

^ 2 

25 

26 

■ 12.349734 

1 2. £36784 

.045424 

3.876601 

3,876601 

3.676601 

3 

4 

27 

-2.8 

1,3.067254 

13.6412H3 

.272692 

.45448 6 

,1.14060 6 
.804061 

1.140606 

.804061 

1.140606 
. 8040 61 

5 

6 

9 

9 

14.071751 

i4.35.flk.aa- 

.590858 
.681679 

.586679 

.516295 

.580679 

.516295 

.580679 

.516295 

7 

a 

31 

32 . 

14.717271 

15.148276 

.795321 
.93183 8 

.5171 91 

. 231261 

.517191 
. 2 3126 1 

.517191 

.231261 



Y-STATTON 2 (V = 

5.9056 

(HFT FRSn 


i 

33 

12.461199 

. D 45 87 8 


3.633942 


z 

_3 

34 

as 

12.737235 
13.1611 kO 

.138698 

.276381 

1.582404 

1.225154 

1.582404 

1.225154 

1.582404 

1.225154 

k 

S 

36 

37 

13.703055 

14.117S till- 

.460765 
.599074 

.864290 

.676378 

.884290 

.676378 

.884290 

.676378 

6 

2 

38 

30 

14.392774 

14.727249 

.691192 

.802935 

.596626 

.593626 

.598626 

.593626 

.598626 

.593626 

8 

40 

15.121034 

.934493 

.277121 

.277121 

.277121 



Y-STATION 3 (Y =* 

7.2223 

(METERS)) 


1 

2 

41 

42 — 

12 . 60 7 7 31 
12.668745 

.046714 
.140 65 8- 

3.576674 

1.596329 

3.576674 

1.596329 

3.576674 

1.596329 

3 

4 ' 

43 

44 

13.260268 

13.782279 

.281575 

.469456 

.1.255788 

.924902 

1,255788 

.924902 

1.255788 

.924902 

5 

16 

4S 

---'46 

14*173799 

14.434743 

.610 371 
.704290 

.712820 

.630955 

.712820 

,630955 

.712820 

.630955 

- • ':i,7 .. 

?- -• • • 

... /• -HJ< 

: -48 

14.738016 

-15.084022 

.813444 

.937978 

.622751 

.298002 

.622751 

.298002 

.622751 

.298002 

* . ' r ' 

• ■tfia©.:' 

y-STAT-IOM -4. - t v a 

8.5201 

(METERS)) 


* i * * - *», 


_ig'.7j?1932 

i> 

' .047523 

3.629523 

• 3.629623 

3.629623 







1 . CIO. 37 

3 

5i 
$ 2 

13* 365797 

-1 T . flKft A7A 

• 286715 

IIi93954 

1.293954 

,95194? 

1.293954 

,951542 


53 

5h 

!<**22515<f 

4 tx - 1.7 AA ftC 

.621562 

_ 7A 79XX 

.725361 

.725361 

,638074 

.725361 

,638074. _ 

7 

-8-- 

55 

S6._. 

14.744285 .823842 

— 1^^044208 .-941486 

• 62 7961 

- .627961 

.402423 

.627961 
— .302423 — 
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TIFS. CASE. It, 


-U . 


_R£V_ISED._HQDES_ 


DATA FOR THIN BODY NUMBER 4 


* THIN BODY NAME TIP-HING 

* ON/OFF PLANE OF SYMMETRY OFF 

..■■NUMBER OF P ANELS - t -— -- -- -- -- r 411 

* NUMBER OF Y-STATIONS 5 


PANEL X X/C CP CP CP 

JiO. NO. (HFTFRS1 SY MM FT R I C RIGHT SIDE LJEJEX JSIflE 




Y-STATION 1 (Y 

= 

9.8687 ' 

(METERS) > 


1 

57 

12.901415 

' .048304 


3.326991 

3.326991 

3.326991 

7 

5 ft 


■1 45768 


1.3ft 5255 

} * 3B5255 

!• 3B5255 

3 

59 

13.472978 

.291944 


1.133156 

1.133156 

1.133156 

4 

60 

■FVTnnT^B 



. . B473-.Z 

.847347 


5 

61 

14.273148 

.633033 


.658125 

.658125 

.658125 

ft 

42 


WJIU11 


mi iiiiii 

TN i i i i 


7 

63 

14.745695 

.834466 


.591697 

.591697 

.591697 

fl 

64 

15. (10519a 

.945865 


.284817 

„ .284817 

.284817 . 









66 

1 3.049673 

. ±i 4fl 52 3 


3. 371692 

. 3.371692 

3.371692 

2 

66 

13.258289 

.146736 


1.414688 

1.414688 

1.414688 

3 

ft 7 

1 3.571310 

iRX L'l ilcl 



1.158699 

. 1.158699— 

4 

68 

13.988682 

.490429 


,860510 

.860513 

.860510 

6 

60 

1 6. 38 1 7 56 

.637749 


mm mm 

Mil — 


6 

70 

14.510416 

.735935 

■ 

.582305 

.582305 

.582305 

7 

71 


. B3883 A 


.588938 

.588938 

— i n n i m 

8 

72 

14.957633 

.946472 


.262833 

.262833 

.282833 

- 


Y-STATION 3 (Y 

X 

12.5394 

(METERS)) 


7 1 

73 

• 13.193857 

,048639 


3.324212 

3.324212 

3.324212 

2 

74 


,147.399 





3 

75 

13.665264 

.295489 


1.147312 

1.147312 

1.147312 

4 

76 

14.042410 

.492979 


.851411 

— JU1W 


5 . 

77 

14.325360 

.641144 


.649983 

.649983 

.649983 

t, 

?a 


bSVvJl t T.tT. 


A trt: Lrl 

i n i mm 

— mill i 

7 

79 

14.709035 

.842053 


.582296 

.582296 

.582296 

a 



Hl'l .* :TFTW 


■ 1 IIIIII 

.276879 


-f 

7* ‘ 

Y-STATION 4 (Y 

X 

13.8572 

(METERS)) 


.•.V' .'iw'".- 

fti 

13.338050 

.048730 


3.176815 

3.176815 





1 1 1 8 0 2 5 


4 . 

HKKUWMI 

^nr?:i ht— 

3 

■"O c " 

63 

13. 758630 

.296970 


1.099952 

1.099952 

1.099952 




, ^g «7£ 




- T >1168 39 

5 

85 

14,3475 82 

.644587 


t 

.620825 

.620825 

.620825 



5JL-5d3^ 



, 54 At) ft it 

,64aor4 

.5481164 

7 

87 

14.687618 

* 84528 6 


,562789 

.562789 

.562789 

.6- — ■ 

-82- - 

-14.86 2 895 

— *948740 


.266846. 

.266646 

*2.66846 — 


X=ST AI-IOU -5- -I Y._ = 15.2704 ..{HETERSl-l- 


1 . 8-9- 13.. 4925 83 „ .046301 2*76202-7 2 . 762(177 2* 7 62 077 . 

2 9 C 13.636635 .148237 1.154320 1.154320 1.15432D 

3 . ai- ...13*857624 .298133__ *955560... .956580 .956560 



* 92 

Q3 

14.150107 

i a 3 lift 

.498017 
. 6i»79? fl 

.712907 
. 55 2433 

. .712907 

.552433 

.712907 

.652433 

b 

94 

14,515425 

.747847 

.522601 

.522601 

.522601 

7 

q c 


- 848569 

. 451678 

.451 67ft 

.451678 

8 

96 

14.811067 

.953C41 

.222237 

.222237 

.222237 
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TIFS CASE 1 » -LANDING- CONFIG URAT ION* .REVISED MODES 


DATA FOR THIN 30DY NUMBER 5 


* THIN BODY NAME 

* -DIHEDRAL-ANGLE -r — 

* ON/OFF PLANE OF SYMMETRY 


TOP- SFS 

96.5C.ia »OEG. 1. - 
OFF 


* NUMBER OF Y-STaTIONS 


PANEL X 

— NH»— tHETERS) ... 


CP 

. SYMMETRIC 


CP 

. RIGHT SIDE 


CP 

LE£T_SJCQE_. 
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TIFS CASE 1 , lANQiNG_£DN£HaiRA.T10iU_R£VIS£Q..M0QE3. 


OttTA FOR THIN BODY NUMBER 

7 



.. * nxwpnPAi ANGLE ---------------- 

-nnnn (rifr,.) 



* ON/Grr rLANfc. Ur jTnntIKT m -- -- - 
5 MIIMRCg DP PANELS 

±2 


* NUMBER OF Y-ST AT IONS 7 

PANEL X X/C CP CP CP 

.__NQ» NO. (METERS) SYMMETRIC RICMT_S THE LEFJL_SI_QE. 




Y-STATION 1 (Y = 

1.8599 

(METERS)) 


X 

P 

• 1 J 5 

24. l&w 630 
Pu.51 PPl-3 

.049296 
,211464 . 

.274402 
_ «. 0816-92. 

.274402 

.08169? 

.274402 
.081692 

3 

u 

10 7 
in a 

25.208321 

?5.73i68n 

,*♦8575 0 
.691 964 

.000754 
— . ?2 82 68 

.000754 

-.228268 

.000754 
-.228268 

5 

(y 

10 9 
11 Q 

25.955251 

Ph .PQ .hf. 7 

• 780 C5 6 
.912193 

-.066623 

-.036582 

-.066623 
-.036582. . 

-.066623 

-.036582 



Y-ST AT TQM ? ( Y = 

P.6630 

(HETERS) ) 


1 

1 1 1 

?U. P78597 

.049365 

.596473 

,.5.9647.3. 

.596473 

2 

3 

112 
1 1 3 

24.647538 
?5-?7i 5 5n 

.211517 

•U8577L 

.191823 

.045508. 

.191823 
.OU550 n 

.191823 
.045508 

b 

114 

. 11 <5 

25.740697 
pr t 9ui i „9 

.691968 

.780050 

-.218478 
.. -.054499 

-.218478 

. , -.054499 _ 

-.218478 

-.054499 

• *;•* * 6 ■ 

116 

26.241728 

.912175 

-.028189" 

-.028189 

-.028189 

. *+* ' ■ 


Y-STATION 3 (Y = 

3.4163 

(METERS) ) 


1 

? 

117 - 

■1 1 A 

24.445SC8 

pu.77t.uf.'; 

.049123 

.211333 

1.210552 

.358009 

1.210552 

.358009 

1.210552 

.358009 

3 

119 

i?n 

25.33L850 

P5.7U9153 

.485688 

.691956 

.10 1203 
-.199463 

.101203 — 
-.199463 

.101203 

-.199463 

5 

6 

121. 

ip? 

25.927846 

26.195835 

.780070 
.912241 

-.034434 

-.014823 

-.034434 
-.014823 

-.034434 
=UP14823 . 



Y-STATION 4 1Y = 

4.1863 

(METERS)) 


1 . 

1 P 3 

24.616144 

.048861 

1.454053 

1.454053 

1.454053 

2 

3 

124 

1P5 

24.904225 

25.391473 

.211134 

.485595 

.434977 

.131876 

.434977 

.laiazfe 

.434977 

.131876 

4 

5 

12 6 
127 

25.757798 

25.914287 

.691942 

.780090 

-.186354 

-.021296 

-.186354 

-.021296 

-.186354 
: -,124296. , 

6 . 

12 8 

26.149019 

.912313 

-.005700 

-.005700 

-.005700 

' 4- ~ i s 

i ' ; 

Y-STATION 5 CY * 

4.9017 

(METERS)) 



■P- 



PANEL 

X 

X/C 



CP- . 

LEFT -SIDE- 


1 

12 9 

24.774660 

.048898 

1.027734 

1.027734 

1.027734 




1 *11 1 * 9 

T 337U60 

. 337UC9 

mm 

3 

131 

25.447790 

.485609 

.112249 

.112249 

.112249 


1 7 P 



— T 1 87U07 

-.187407 


5 . 

133 

25.901690 

.780088 

-.024021 

-.024021 . 

-.024021 


L3-4 



-.006762 

-.006762 
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-j* m 


_11S T q4fl q Qft 

136 25.135465 .21Q864 

—13 Z £5 *-4.935-03 ,435.46.9 

138 25.773204 .691924 


.278050 


-.1970 93 


.278050 


-.197093 


.276050 

.n*6na«; 

-.197093 


26.065502 


.912410 


-.013098 


-.013098 


Y-STATION 7 ( Y 


6.2311 (HETERSH 


-.013090 


25.069245 


.047624 


25 . 552450 .465156 

25.78^.764 . ... .6.91830- 
25.678262 .780189 

26. i, 245 Z3 .9126 53 


.783135 

. ?3f,1 S4 

.057298 
209 8-33— 
-.045924 
021131.. 


.783135 


.057298 

-.045924 


.783135 

.£ 36154 . 

.057298 


-.045924 
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APPENDIX D 

COMPLEMENTARY FILTER CORNER FREQUENCY DETERMINATION 
COMPLEMENTARY FILTER CONCEPT 


A complementary filter, block diagram is presented in Figure D-l. 



T(s)j 

T(s ) 2 

e c 

(S) 


Signal obtained from air pressure measurements 
Signal obtained from inertial measurements’ 
Transfer function of filter (1) 

Transfer function of filter (2) 

Composite output signal 
Laplace variable 


Figure D-] BLOCK DIAGRAM FOR ARBITRARY INPUTS e AND &■ 

T & 

The signal flow depicted in block diagram form in Figure D-l is described by 
Equation D-l. 

e c = e T(s ) 1 + e- L T(s) z (D-l) 
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Also since e c , e p and e L are analogs of the same variable. Equation D-l can 
be rewritten as Equation D-2. 


1 = T(s) 1 -h T (s) 2 . (P-2) 

The transfer functions of the complementary filters are forced to yield a sol- 
ution to Equation 2. That is, one transfer function is assumed and the other 
is computed from Equation 2. The resultant composite output signal has no 
"bumps" at or around the frequency that output data transitions from one input 
signal to another. A simplified example is presented to further illustrate 
this point. For the block diagram illustrated in Figure D-2, <5 and 

the transfer functions of filter (1) and (2) are ~7'(s) 1 and T(s) t respectively. 



Figure D-2 COMPLEMENTARY FILTER BLOCK DIAGRAM FOR e^= e; 

Assume a transfer function for T(s) j and compute T(s)g from Equa- 
tion D-2. 

From Equation 

» e p T(s) j -f e t T(s) z 

also, =■ - e. 

Therefore, = & \_J r (s) 1 i-T(s)^ (0-3) 

From Equation D-2. 

T(s\ = 1 - T(s) 1 

Substitute T(s) 2 from Equation D-2 to Equation D-3 

= e ^T(s) t + 1 -T(s ) f J = e J 
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Thus, it is clear that the composite output signal precisely equals the input 
signal in amplitude and phase for all frequencies within the bandwidth of the 
computing electronics. 

CORNER FREQUENCY DETERMINATION 

The following procedure can be used to determine optimum corner 
frequencies for the complementary air data filters. Power spectral density 
functions describe the errors in each input signal and the square of the stan- 
dard deviation (energy) for each signal is computed from the integral (over all 
frequencies) of the density function. The density function for the output sig- 
nal is computed from the product of the input density function and the squared 
absolute value of the filter transfer function. The integral (over all fre- 
quencies) of the output density function results in a relationship between the 
square of the output error standard deviation and the filter corner frequency. 
Differentiating this equation and solving for the filter corner frequency for 
the derivative equal to zero yields an optimum corner frequency. To further 
clarify this concept, the air data vertical speed complementary filter corner 
frequency computation is outlined. The. system input signals, filter transfer 
functions and output signal are shown in Figure D-3. 



Figure D-3 BLOCK DIAGRAM OF AIR DATA ALTITUDE RATE 
COMPLEMENTARY 'FILTER 

where 

hp = altitude computed from pressure measurements 

hz ~ vertical acceleration computed from body axis 
accelerations and attitude angles 

h = composite altitude rate signal 

h results from the body to earth axis transformation of the body axis 
accelerations , /?,, , and and by pitch and roll anele. 

* 7 J • - • 
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In general form: 


hr = -f(”x s n y , O, 4) 

An incremental change in h T is equal to the total differential, that is: 


•• 'df . Sf . 

Ah x = o An x. + sTT Ar) n 

x 9n x * drty y 


7 dO 9<p 


•CD-4) 


The partial derivatives of Equation D-4 

transformation matrix, and the differential variations in each signal repre- 
sent sensor errors. The significant sensor errors used for this computation 
are static offsets. The magnitude of the offsets used are: 


An z - i S.O x 10 ^ g's A 9 - t 1 deg 
Afiy = + 1.0 x 10 3 g's A 4 - t 1 deg 
An a = t 2.5 x 10~ 3 g's 

This defines the errors in h ^ and results in an impulse when transformed to 
a power spectral density function. 


- The errors in h ^ result from air turbulence, and power spectral 
density functions describe the cl r if and uf velocity components of the tur- 
bulence. The assumption was made that the major component of error is intro- 
duced by, and subsequent analysis was performed on, the Y axis gust term. The 
Dryden form of the spectrum for the turbulence velocity is used in this analy- 
sis. The velocity spectrum described by the power spectral density function 

= °u Jr 7 [ i-i-(£uj) z ] is converted to pressure variations at the pitot 


£ 






Y M 

*2 

1 

A H ) 

7 T 

. ** (■?«>)* . 


where: 


L = 1000 ft altitude 

1/ = 200 ft/sec aircraft velocity 

Vy ~ 20 ft/sec mean value of wind velocity in Y axis 

& er . ~ 3 ft/sec standard deviation of wind velocity in 
Y axis 


The altitude rate complementary filter is redrawn in Figure D-4 with input and 
output power spectral density functions labeled. 
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Figure D-4 ALTITUDE RATE COMPLEMENTARY FILTER WITH ERROR 
DENSITY FUNCTIONS AS INPUTS AND OUTPUTS 

(co) 

^ = power spectral density function of the errors in the 

50 pressure derived altitude signal 

- power spectral density function of the errors in verti- 

h r cal acceleration 

$(co) - power spectral density function of the errors in the 

^ composite output signal 

The standard deviations of the input and output signals are: 

% ’Z, s< C dM CD " 5 

- r 00 $<*>) m_< 


h ttco 


$(co) related to ^ and £(co) ^y. 

h h f> h T 
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Equation D-9 is the standard deviation squared of the error in h as obtained 
from the operations specified by Equations D-7 and D-8. 

To obtain an optimum Z for the filter,-^, j is obtained. 
This expression is set to zero and solved for T . 9 


9 
3 Z 



= 0 


erf? is plotted in Figure D-5 to gain additional insight in the shape of the 
error function and how sensitive the error function is to variations in filter 
comer frequency. As can be seen, the function has a rather broad minimum and 
an acceptable result is obtained for filter time constants that range from 
0.5 to 5 seconds. A value of 5 seconds is used for the air data altitude rate 
complementary filter. The development of the altitude, airspeed, and flight 
path acceleration signal filters is similar to the .procedure riutlined for the 
altitude rate system. . . . 
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APPENDIX E 

COLLECTIVE AILERON SYSTEM PRELIMINARY DESIGN 


E. 1.0 GENERAL 

This appendix describes the preliminary design of a servo driven 
collective aileron system to be used for maneuver load and structural mode 
control on the AF/TIFS airplane 

E. 2.0 MECHANICAL DESIGN. 

In this section, the existing TIPS aileron control system is de- 
scribed briefly to provide a better understanding of the mechanical arrange- 
ment of the proposed Class II modification. 

E.2.1 Existing TIFS Aileron Control System 

The existing TIFS aileron control system is shown in Figure E-l. 
In this configuration, a centrally mounted electrohydraulic position servo 
actuator, when the variable stability system is engaged, moves both ailerons 
through the standard C-131H cable and drum system. The multiple drum assem- 
bly shown in this figure is mounted behind the rear spar at the center of the 
aircraft (buttock line 0) and is a standard C-131H component. The mechani- 
cal system from control wheels to stops at the aileron quadrants was stressed 
for a condition in which two pilots apply couples at the control wheels in 
conjunction. Pilot applied wheel loads are given in this reference as two 
tangential limit forces of 80 lbs each applied in conjunction at the wheel 
rim. The TIFS position servo has relief valves set at 535 psi differential 
pressure across the piston which limits servo output forces to 525 lbs. This 
is equivalent to one pilot applying two tangential forces of 63.28 lbs in 
conjunction at the wheel rim. The maximum servo output therefore, represents 
79.1% ox the maximum single pilot applied limit force and 39.5% of the 
maximum two-pilot applied limit force. The control system from the stops 
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Figure E-l TIFS AILERON CONTROL SYSTEM 



to the aileron is designed fox an aileron down limit airload of 410 ft lbs . 

It should be noted that this same reference assumes that the pilot applied 
forces at the center drum are distributed unequally 2/3 to one aileron and 
1/3 to the other. Using the same assumption, the maximum limit airload which 
may be applied to one aileron by the position servo with relief valves set at 
535 psi is 525 x 2/3 x 475/1.64 x 1.92/12 = 162.19 ft lbs. 

There is a geared tab on each aileron. The gearing provides hinge 
moment relief (i.e. up tab motion for the down moving aileron and vice versa). 
For each degree of aileron movement, the tab will move approximately the same 
amount in the opposite direction. The tabs are adjusted for trim purposes by 
operating a knob on the center pedestal of the pilot's cockpit. Motion is 
transmitted via gears, pulleys and cables to jackscrews in each wing. Each 
jackscrew actuates a tab via two push-pull rods and a rocking lever. The max- 
imum aerodynamic hinge moment produced on a tab is 18 ft. lbs. (Reference E-l, 
page 167) and trim movement on the tab is limited to ± 7°. 

The TIFS airplane does not have the standard C-131 aileron -rudder 
interconnect springs. Consequently, manual operation of the TIFS aileron con- 
trol feel is different from that of a standard C-131. There is no spring 
resistance in the TIFS system, but there is a higher friction level. This 
is due partly to the drag of the aileron servo which is in series in the sys- 
tem at all times and partly to the re-routing of the aileron cable controls 
in the wings through the direct lift flap areas. The effect of the proposed 
modification on the system friction is analyzed in Section E. 2.2.4. 

I 

E.2.2 Preliminary General Design of Collective/Differential Aileron System 

It is proposed to modify the aileron system as shown in the 
schematic diagram, Figure E-2. The major changes will be: 
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1. To remove the existing central servo actuator and install a 
servo actuator in each wing driving each' aileron from an 
adjacent position. 

2. To remove the standard C-131H multiple drum assembly at the 
center of the fuselage and install in its place a new double 
drum assembly. 

3. To install a coupler/ differential unit in the fuselage (sta- 
tions 243 to 281) connected to the cockpit control wheels 
through the existing cable system and connected by divided 
cable systems to the new double drum assembly at the center 
of the fuselage. 

4. To install hydraulic connections between the existing position 
servo subsystem and the new wing servos and coupler/differen- 
tial unit. 

Brief descriptions of the above changes follow: 

E. 2.2.1 Wing Mounted Aileron Servos 

The proposed installation is shown schematically in Figure E-3. 

The actuator bodies will be mounted on the aft side of, the rear spars. A new 
double quadrant will be mounted inboard of each actuator. Motion of an 
actuator piston rod will be transmitted to the aileron quadrant via a short 
run-around cable system. This mechanical arrangement has been selected to 
give good system dynamics in either direction with a minimum of structural 
rework in this area. The new double-quadrant support ribs will be braced to 
the existing quadrant support ribs. The trailing edge structure will be 
modified to provide clearance for the servo installations inside and to provide- 
large access doors on the underside. The necessary servo control valves, hinge 
moment limiting valves, maniford blocks, etc. will be installed in the general 
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Figure E-3 AILERON SERVO ACTUATOR IN WING 


area o£ the actuators. All hydraulic piping will be solid type, no flexible 
hose being required. 

E. 2.2.2 Double Drum Assembly, Rear Spar, Buttock Line 0 

To provide separate cable drives between 'the coupler/ differential 
units and the ailerons, it is necessary to replace the existing solid drum 
assembly by a new double drum assembly. This new assembly will be inter- 
changeable with the existing C-131H assembly, but separates the cable drives - 
to the two wings to permit independent aileron operation. Existing pulley 
assemblies mounted on the top aft corner of the rear spar will be used for 
the two cable systems. - ’ 

E.2.2.3 Coupler/Differential Unit 

This assembly. See Figure E-4, has two • major functions. When the 
airplane is being flown manually using the normal aileron system, a rigid 
mechanical connection is provided between the left, and right ailerons. When 
the airplane is being flown on the ACT system, the left and right ailerons are 
mechanically disconnected at the coupler/differential unit and are then inde- 
pendently controlled by the wing-mounted actuators. Through the action of the 
differential linkage, the cockpit control wheel motions will reflect only the 
rolling deflection of the ailerons. The unit consists basically of three para- 
llel shafts of which the forward one is connected to the pilot control wheels, 
the center one to the right aileron and the rear one to the left aileron. The 
coupler units operate on the- center and rear shafts only and the differential 
linkage connects to all three shafts. Before describing the purpose and action 
-of the differential linkage the coupler unit operation must be explained. 

E. 2. 2. 3.1 Coupler Units 

The internal arrangement of the coupler units is shown in Figure E-5'. 
The coupler bodies are trunnion-mounted on levers attached to the rear shaft. - 
The coupler rods have slotted ends working on rollers between levers mounted 


241 




PRESSURIZED: FULL COLLECTIVE AILERON CONTROL 



PRESSURIZED: COMBINED COLLECTIVE AND DIFFERENTIAL 


Figure E-5. COUPLER UNITS 


2 43 



to the center shaft. Free single-acting hydraulic pistons inside the bodies 
are spring loaded against the hydraulic action. If the hydraulic chambers are 
depressurized,. the., spring forces- on- -t-he- pistons drive them to the "locked" end 
of their travel. The enlarged diameter portions of the rods are then 
trapped between piston and body. Linear movement of the rods relative to the 
coupler bodies is then prevented. Since the rollers on the center shaft are 
at the forward ends of the coupler rod slotted ends, the center shaft must 
then rotate in unison with the rear shaft. When- hydraulic pressure is applied 
to the pistons, they are driven to .the forward end of their travel compressing 
the springs more and permitting the coupler rods to slide freely between 
pistons and ends of bodies. The combination of coupler rod movement and slots' 
then permits the center and rear shafts to rotate independently in the same 
or opposite directions as required by the movements of the servo actuators 
mounted in the wings near the ailerons. The springs in the coupler units are pre- 
loaded to exert 1289 lbs on each piston in the unpressurized "locked" position. 
Spring force increases to 2160 lbs when hydraulic pressure moves the piston to 

. 2 

the opposite end of its travel. -The net piston area is 1.96 in so that a 
pressure of 2600 psi will produce 5096 lbs of force on the piston. The lowest 
pressure which will balance the spring in the "unlocked" position is. 1102 psi. 

The supply of hydraulic power to- the uncoupling mechanism will be controlled 

v ^ # 

by a solenoid operated three-way valve. Energized',' this valve will pressurize 
the mechanism. De-energized, pressure will be shut off and fluid in the mech- 
anism will be bled through the flow restrictor to the hydraulic return line. 

See Figure E-14. The rate at which the free- pistons can move is cpntrolled by 
a restrictor. Total piston displacement for the two couplers is 2 x 1.96 
x 1.6 = 6.27 in^. If the time is 0.5 second, then the average flow rate through 

3 

the restrictor will be 12.54 in /sec or 5.25 gpm. 

E.2.2.3.2 Differential Linkage 

The action of the differential linkage is illustrated in Figure E-6 
for the "coupler locked" configuration and in Figure E-7 for the "coupler un- 
locked" configuration. The design of this linkage is such that in the coupler 
unlocked condition, when the ailerons are being independently driven by servos. 
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CONTROLS AND AILERONS IN NEUTRAL POSITION 




Figure E-6 NORMAL AILERON OPERATION - COUPLERS UNPRESSURIZED AND LOCKED 

' i 
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COLLECTIVE OPERATION, AILERONS DOWN 
NO PILOT CONTROL FEEDBACK 



COLLECTIVE OPERATION, AILERONS UP 
NO PILOT CONTROL FEEDBACK 



COLLECTIVE AND DIFFERENTIAL OPERATION 
PILOT CONTROL FEEDBACK DIFFERENTIAL ONLY 


igure E-7 COLLECTIVE/DIFFERENTIAL AILERON OPERATION 
COUPLERS PRESSURIZED AND UNLOCKED 
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the cockpit control wheels will display only rolling motions of the ailerons to 
the pilots at all- times. The advantage of this is that while the ailerons are 
being moved by the servos, the normal use of aileron wheel motion as an indica- 
tor of roll input will be preserved. Total aileron deflection will be displayed 
to the pilot by meter movements similar to those used in the TIFS direct lift 
flap and side force servo systems . 

E.2.2.4 Control System Friction 

- Measurements made during TIFS Flight #486 indicated that the force 
required at the rim of the control wheel to overcome breakout friction is about 
14 lbs. This represents a force in the control cable of about 58 lbs and was 
measured -with the aileron servo actuator in an unpressurized condition. The 
actuator is a Miller floating seal type modified at Calspan. This type was 
chosen because of the relatively low friction level which can be achieved with 
it. Actual bench tests of the aileron servo gave breakout force values of about 
8 lbs for the unpressurized condition and 35 lbs 'pressurized. It may be con- 
cluded then . that the friction in the aileron system due to cable bending, 
pulleys, bearings, pressure seals and rubbing strips must be in the order of 
50 lbs when the cables are tensioned to standard C-131H requirements. Compared 
with the existing system, the new collective aileron control system will have 
increased friction in it, estimated as follows: 


existing TIFS aileron system 58 lbs 

additional aileron servo 8 lbs 

\ 

coupling/differential unit (aileron mode) 5 lbs 

Total 69 lbs 


This represents about 16.5 lbs force at the control wheel rim, an 
increase of about 18% over the forces presently experienced in flying the 
TIFS airplane manually. 

From the above estimate and discussions with the regular TIFS pilots, 
it is predicted that the handling qualities rating for ILS tracking or cross- 
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wind landing tasks will not be degraded from the current "acceptable but "unsat- 
isfactory" level. In the event that handling qualities are, after preliminary 
tests", j'udged~" T, unacceptable" , a power boost will be added to the aileron system. 

E.3.0 ELECTRICAL AND HYDRAULIC CONTROL DESIGN ’ 

The existing TIFS VSS electrical and hydraulic control systems 
must be understood in order to study the impact of the proposed collective * 
aileron modification. 

E . 3 . 1 Normal • Engage Sequence 

The engagement procedure applies electrical and hydraulic power- to 
the Variable Stability System. Four time-sequential operations are required to 
engage the VSS. These steps are: 

1. Depress Master Power Switch. (SAFETY PILOT) 

2. Depress Feel System Engage Switch. (TEST PILOT) 

3. Depress Pressurize Switch. (SAFETY PILOT) 

4. Depress Engage Switch. (SAFETY PILOT) 

These four engage steps are' electrically interlocked and the system can only 
be engaged if the switches are activated in the order shown. Indicator lights 
display the engagement state in the simulation cockpit, the safety cockpit and at 
the test engineer's consoles., 

The four-step engage sequence is illustrated in Figure E-8 and is 
described in more detail below. 

1. Master Rower Switch — applies electrical (AC and"DC) power to 
the variable stability system electronics. Input power is also 
applied to the VSS sensors. Power is connected to the ready bus 
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Figure E-8 
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through an 80 second time delay. The time delay allows the 
attitude gyro to erect before the system can be engaged^. Power . 
from the ready bus is applied to the’ feel system activate switch 
through the safety trip system.' 

2. Feel System Engage Switch — applies hydraulic power to the feel 
system, i.e., pressurizes the hydraulic line to the feel system 
actuators and simultaneously closes the bypass valves across 

the feel system actuators. In addition, the feel system hydraulic 
flow control valves are connected to the servo amplifiers. At 
this point the feel system is engaged. Power is applied to the 
pressurize switch. 

3. Pressurize Switch — applies hydraulic pressure to the surface 
servo hydraulic lines. However, bypass valves across the servo 
actuators remain open. (See Section E.3.2). Power is applied 
to the contacts of the engage switch. 

4. Engage Switch — closes the bypass valves across the servo actu- 
ators of the servos that have been manually preselected before 
system engagement or selected after system engagement. Simul- 
taneously, the surface servo hydraulic flow control valves 

are connected to the respective servo amplifiers. 

E.3.2 Existing Hydraulic Servos and Safety Features 

The aileron position servo presently installed in the TIPS airplane 
is powered from a special variable stability system (VSS) hydraulic supply. 

The three surface servos, elevator, rudder and aileron, form a subsystem con- 
trolled by a shutoff valve in the pressure line. This valve, normally closed, 
energized to open, works in conjunction with a bypass valve which connects the 
subsystem pressure line to the return line. This bypass valve (Bypass I on 
Figures E-9, E-10, and E-ll) is normally open and energized to close at the same 
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time that the shutoff valve when energized, opens the pressure supply to the sur- 
face servos. The following description of the various hydraulic servo states 
applies, therefore, to all three surface servos... figure .E-9 shows the state of. 
the hydraulic section of- the' servo when the VSS is OFF. The safety pilot, by 
moving the control, forces oil from one side of the actuator or the other through 
the normally open solenoid valve. Upon pressurizing the system, the solenoid 
valve between the pump and the Moog valve opens and leakage flow through the Moog 
valve raises the pressure inside the actuator as shown in Figure E-10. However, 
the solenoid valve across the struts (Bypass II) has remained open. So the --safe- 
ty pilot can still move the controls forcing oil through the open solenoid valve. 
When the system is engaged, the solenoid valve across the actuator closes (Figure 
E-ll). - The computer can now control the actuator position- through control of the 
Moog valve. Disengaging the variable stability system returns, the hydraulic 
section of the servo to the state shown in Figure E-9. 

Two groups of safeguards are provided in the existing system to en- 
sure that the flight control surfaces can be easily- operated by a safety pilot 
when the VSS system is "off". The first group consists of the solenoid-operated 
shutoff valve (Bypass II) fitted as a. bypass control in each surface position 
servo piping system (see Figure E-9). When de-energized, these valves are normal- 
ly open, and there are neither differential pressures across the servo actuator 
pistons nor any impediment to the free flow of fluid from one side of each piston 
to the other side. The second safeguard is provided in case any one of the 
shutoff valves just' mentioned fails to open. when de-energized. This consists of 
a single solenoid -operated shutoff valve (Bypass I) connecting the pressure 
supply for all surface position servos to the common return line, and a pair of 
check valves fitted in each position servo manifold block. The Bypass I shutoff 
valve, also normally open when de-energized, dumps the supply pressure, thus 
permitting flow through the check valves. Should any servo Bypass II valve 
fail to< operate, fluid can escape freely' through one of the check valves when 
movement of the relevant surface by the safety pilot causes the servo piston 
to move in the cylinder (see Figure E-12) . 
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In addition, the safety pilot is able to overpower the servo at 
any time by exerting forces higher than those determined by the relief valve 
settings in the existing aileron control system, see Section E.2.1. 

E.3.3 Safety System 

The safety system consists of both manual (dump, buttons) and auto- 
matic disengagement capability. The function of the automatic safety system. 
Figure E-13, is to supplement the action of the safety pilots in situations 
where rapid transients can occur or where parameters cannot be monitored due 
to high workload. The safety trip system protects the airplane by disengaging 
the variable stability system thus restoring normal aircraft control to the 
safety pilots . A safety trip or system "dump" is indicated by blinking red 
indicator lamps at all aircraft control stations and a modulated tone in the 
intercom system. The system dump can be initiated in four ways: 

1. By safety pilots: 

The safety pilots can monitor control surface positions 
and rates of motion; aircraft attitudes, rates and accelera- 
tions and the aircraft environment. If either pilot feels 
and unsafe situation is approaching, he can disengage the 
system by depressing a dump button located on his control wheel. 

2. By simulation pilots: 

The simulation pilots have similar dump buttons which they 
may depress if they cannot maintain control of the simulation. 

3. By test engineers: 

The two test consoles located in the cabin are equipped 
with panel-mounted dump buttons to be depressed if a computer 
malfunction is detected by a test engineer. 
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4. By automatic safety trip: 


The automatic safety trips consist of electronic circuits 
that continuously monitor selected critical aircraft para- 
meters. If any of these signals exceed a predetermined safe 
level, a system dump results. Parameters monitored are: 

a) control surface hinge moment. 

b) control surface rate of motion. 

c) aircraft lateral and vertical acceleration. 

d) variable feel system controller rate of motion. 

e) fuselage, horizontal stabilizer, wing, and vertical fin 
structural loads. 

When any of these monitored variables cause a system dump, 
the safety trip indicators can be deactivated (reset) by either 
safety pilot by depressing the safety trip reset button located 
on the safety cockpit center isle console. 

E.3.4 Collective/Differential Aileron Hydraulic System 

It is proposed to modify the VSS flight control hydraulic subsystem 
as shown schematically in Figure E-14 by removing the present centrally mounted 
aileron actuator, adding two wing-mounted actuators and adding the hydraulically 
controlled coupler/differential unit. Hydraulic power to the two coupler assem- 
blies in the latter unit is controlled by a solenoid-operated three-way valve so 
that the unit may either be pressurized to unlock or depressurized to lock. The 
coupler units will be pressurized only when the aileron servos are selected and 
the engage button is pressed. It should be noted that the two coupler assemblies 
provide fail safe action, for, if hydraulic power is shut off, the chambers 
filled with fluid are then opened to the return line and the pistons being 
spring loaded will return to the locked positions. 
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Figure E-14 COLLECTIVE/DIFFERENTIAL AILERON HYDRAULIC SYSTEM 







E.3.5 Aileron Electrical Control System 

Existing circuitry in the Control Logic Engagement chassis (PNL 11) 
will be used for the engage and disengage functions of the two aileron servos 
and for control of the coupler/ differential unit. Position servo electronics 
for both aileron servos will be contained in the existing surface servo card 
enclosure (Rack 05) . 

A block diagram of the collective aileron control logic is pre- 
sented in Figure E-15. The sequence of engagement and the function of each 
valve will be as outlined in Section E.3.1 with the exception of the coupler 
three-way valve. This valve is energized only if the system is engag.ed and 
the aileron servos are selected. When the three-way valve is activated, the 
couplers function as described in Section E. 2. 2. 3.1. 

The method of applying collective and differential inputs to the 
modified aileron system is shown in the simplified'block diagram. Figure E-16.J 
A more complete block diagram of the mechanization of each aileron position 
servo is presented in Figure E-17. This mechanization is identical t o t hat used 
on the direct lift flap and side force surface servo loops presently installed 
in the TIFS. Position feedback is supplied by an actuator-mounted, infinite 
resolution potentiometer. Differentiated position feedback and high-pass 
filtered actuator- differential-pressure signals will be used for servo stabili- 
zation. The balance and hold circuit nulls the command signals prior to system 
engagement, thereby preventing undesirable transients. 

Hydraulic components such as servo valves and actuators will be 
similar to the type used in the present TIFS VSS systems. The electronic com- 
ponents will be similar to the types presently used on TIFS. 

E.3.6 Modification of Safety System 

The high speed structural limits on direct lift flap deflection as 
a function of speed and normal acceleration will be recomputed as a function of 
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Figure E-15. BLOCK DIAGRAM, COLLECTIVE AILERON CONTROL LOGIC 
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Figure E-16. SIMPLIFIED BLOCK-DIAGRAM OF COLLECTIVE 
AND DIFFERENTIAL CONTROL SYSTEM 
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the collective aileron deflection. These limits will then be implemented, 
as they are at the present time, by safety trip functions. These functions 
will be dependent on collective aileron deflection and single aileron deflec- 
tion as well as the present dependent variables. 

In addition, individual aileron servo force output and surface rates 
will be monitored. 

E.4.0 LOAD ANALYSES 


Load analyses for the hydraulic system and the electrical systems 

follow. 

E . 4 . 1 Hydraulic 

The TIFS VSS hydraulic system is supplied from an engine driven 

3 

pump which can deliver 42 gpm. There are also two SO in accumulators, three 
3 3 

100 in accumulators and one 400 in accumulator in the system to supply peak 
flow demands above the pump capacity. 


The collective/ differential aileron modification will require 
hydraulics to actuate the uncouplers and to drive an additional servo. The 
couplers will require approximately 3.3 gpm for less than 0.5 seconds at each 
VSS engagement. This demand will occur when minimum flow rate is required for 
the servo actuators. ..Thus, the coupler flow rate demand has no impact on 
system capacity. The additional aileron servo will have a peak no-load 
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requirement of 2.2 gpm. This will change the maximum instantaneous peak 
demand from 77.1 gpm to 79.3 gpm. This increase of 2.85% should not signifi- 
cantly affect the operation of the VSS hydraulic system. 

If a power boost should be installed to augment manual operation 
of the aileron control, (a possibility indicated in Section 2.2.4) its power 
requirements will depend on the boost configuration chosen. If one of the 
servos which will be mounted in the wings for this modification also functions 
as a boost cylinder, there will be no change in demand. If a separate boost 
actuator is installed between the pilot control wheels and the forward torque 
shaft of the coupler/ differential unit, it would have a 2.2 gpm instantaneous 
peak demand but it would only operate as a boost cylinder when the two wing 
mounted servos were unpowered so that the power required would be less than 
when the VSS system is operating. 

E.4.2 Electrical Load Analysis 

The increase in electrical power requirements for the aileron mod- 
ification is estimated to be less than 5 amperes of MAIN DC BUS power and 0.1 
KVA of 115V, 400 Hz power. Should an aileron boost system be installed, an 
additional 4.0 amperes of MAIN DC BUS will be required. The estimated break- 
down of power requirements follows, but in all cases, the increase in loading 
is well within the TIES excess capability. 

COLLECTIVE AILERON MODIFICATION 


System Condition 

VSS power off 

VSS power on (not engaged) 

VSS engaged with aileron servos 


Power Required 
None 

2.0 amps MAIN DC BUS 
0.1 KVA L.H. Alternator 
selected 5.0- amps MAIN DC BUS 

0.1 KVA L.H. alternator 
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POSSIBLE AILERON BOOST SERVO 


System Condition 



Power Required 

VSS power 

off 



4.0 amps. MAIN DC BUS 

VSS 

power 

on (not engaged) 


4.0 amps MAIN DC BUS 

VSS engaged with aileron servos 

selected 

l'.O amps MAIN DC BUS 

The impact 

of the additional loads to Table III of the current 

load analysis will be 

as follows: 




1 . 

VSS FLIGHT ELECTRICAL LOADS 

(ADDITIONAL) 

LOAD 


L.H. ALT. 40 KVA . 

DC 

GEN. 375 AMPS 


Taxi 

Term. Area 

Cruise 

Taxi 

Term. Area Cruise 

Basic 

— 

— 

— 

4.0 

1.0 1.0 

VSS 

0.1 

0.1 

0.1 

2.0 

5.0 5.0 

Former Total 

11.73 

11.73 

11.73 

333.3 

329.6 312.9 

Load 






Former % of 

29 

29 

29 - 

89 

88 83 

Capacity 






New Total Load 

11.83 

11.83 

11.83 

339.3 

335.6 318.9 

New % of 

30 

30 

30 

90 

89 85 

Capacity 






2. 

FERRY : 

FLIGHT ELECTRICAL. LOADS (ADDITIONAL) 

LOAD 



DC GEN. 375 

AMPS 




Taxi 

Term. Area Cruise 

Basic 



4.0 

4.0 

4.0 

Former Total Load 


312.1 

302.0 

270.3 

% of Capacity 



83 

81 

72 

New Total Load 



316.1 

306.0 

274.3 

% of Capacity 



84 

82 

73 



3. ELECTRICAL LOADS WITH LEFT ENGINE FAILED 


TWO 200 AMP TRU'S 




Taxi 

Term. Area 

Cruise 

Basic 


4.0 

4.0 

4.0 

Former Total Load 


246.1 

239.0 

222.3 

Former % of Capacity 


62 

60 

56 

New Total Load 


250.1 

243.0 

226.3 

New % of Capacity 


63 

61 

57 

4. ELECTRICAL 

LOADS 

WITH MAIN DC 
No Change 

BUS DISABLED 


5 . ELECTRICAL 

LOADS 

LOADS 

WITH MAIN AC 
DC 

Taxi 

BUSES DISABLED 

GEN. 375 AMPS 
Term . Area 

Cruise 

Basic 


4.0 

4.0 

4.0 

Former Total Load 


303.4 

301.7 

270.0 

Former % of Capacity 


81 

80 

72 

New Total Load 


307.4 

305.7 

274.0 

New % of Capacity 


82 

81 

73 


E.5.0 STRUCTURAL INTEGRITY 


The proposed collective/differential aileron modification will not 
pose any severe structural problems. The servo actuators will be force limited 
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by relief valves as shown in Figures E-9, E-10, E-ll and E-12. The relief valve 
settings are yet to be determined. The quadrant assemblies in the wings are 
designed for two different conditions. The first involves two pilots exerting 
maximum design force on the control wheels with the aileron quadrants being 
held against stops. The case has a cable tension force of 996 lbs ultimate. 

The second case involves a maximum aileron download of 410 ft lbs limit which 
requires a cable tension force of 1325 lbs ultimate to react. An actuator 
with dimensions similar to that used on the existing centrally mounted servo 
would require a differential pressure of 899.5 psi to produce an output force 
of 883 lbs (1325/1.5 = 883 lbs limit). 

The wing servos would be mounted on the aft side of the rear spar. 

The servo body would be rigidly attached to the spar. Servo output forces 
would appear therefore as shear forces along the spar web and as a moment in 
the wing construction plane approximately. The moment would.be reacted as forces 
normal to the spar web and adequate local vertical structure is necessary to 
transfer these forces into the spar caps and the upper and lower wing skins . 
However, existing rib structures on the forward side of the rear spar may pro- 
vide adequate support without" reinforcement . 

The double drum assembly is interchangeable with the existing drum 
assembly and loads in this area will not exceed existing structural capacity. 


’ The coupler/differential unit will be mounted along the center of 
the aircraft between stations 243 and 281, and between the underside of the 
floor and the topside of the radome box structure. It is not anticipated that 
detail stress analysis will produce any major problems in transferring local 
loads into the fuselage structure at this point. 

Some new guide pulley brackets will be required fore and aft of 
the coupler/differential unit. All pulley brackets and other components involved 

in this modification will be designed in accordance with various pilot load 
conditions. 
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. When the ailerons will be operating collectively, they will be 
augmenting the effect of the direct lift flaps. It will be necessary, therefore, 
to change the safety trip signals on the latter in order to stay within the 
load factor limitations which result from the combination of wing bending capa- 
bility, effect of the cargo door cutout on fuselage bending strength, horizontal 
tail constraints and aerodynamic buffet boundaries. New boundary diagrams 
for the load factors will have to be prepared. 


E. 6.0 WEIGHT AND BALANCE 


The estimated net additional weight will be about 150 lbs. Removal 
of the center servo and its valving will reduce weight by about 30 lbs at 
station 390. The installation of the wing servos and new center drum will add about 
95 lbs at station 468. The coupler/differential assembly will add about 85 lbs 
at station 263. The net effect on the C.G. will therefore be: 

' lbs. ins. 

center servo removed 30 lbs at 390 = -11,700. 

wing servos added 95 lbs at 468 = +44,460 

coupler/diff .unit added 85 lbs at 263 = +22,355 ' 

+55,115 

The added weight represents an increase of about .31% of zero fuel 
weight and the moment change is about .3%. The effect on the C.G. of the TIPS 
.airplane will be negligible. 

E . 7 _ STRUCTURAL/ AERODYNAMIC TESTS 

The aileron control system will be subjected to proof loading tests 
after the collective/differential control modification is completed. Ground 
vibration tests of the modified aileron system will also be performed to deter- 
mine, among other things, the effect on wing bending modes of the added masses 
of the actuators out at the ailerons. 
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' The collective ailerons at full trailing edge up deflection are 
predicted to raise the stall speed by- about 5 knots. The actual stall speeds 
will be determined by flight 'testing.. Lateral control at the stall with full 
up and full down collective aileron deflection will also be investigated in 
flight . 


The wing portion of the ground vibration tests in support of the 
flutter flight tests which were made in 1971 will be repeated with the new 
aileron servos engaged to define any change in wing modes. The airplane will 
also be flutter checked in flight with the servos engaged out -to the limit 
•speed of 295 knots. 

The aileron hinge' moment characteristics predicted for the airplane 
' in.Reference 13 wi ll be used for de sign of new stru cture. These will be com- 
pared with the total aileron hinge moment measured in flight by the differen- 
tial- pressure sensor.in the present aileron actuator and later with the indi- 
vidua! aileron moments as determined by similar sensors in the new actuators. 
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